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FORECAST 


the physical, psy- 


Laboratories for space science at Martin are now studying and forecasting 
chological, and biological factors that will affect man in space... another tremendously 
fascinating program which attracts persons with exceptional professional abilities. If you have 
these abilities, you are invited to communicate with N. M. Pagan, Director of Technical and 


Scientific Staffing (Dept. 4-D) The Martin Company, P. O. Box 179, Denver 1, Colorado. 
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Normal binary eutecti 3 ; | A New binary eutectic 
structure (250X) structure (250X) 
+4 Opportunities for + 
ienti i - * 
We, at the Research Laboratories, are convinced that it is t scientists in the phys- ¥ 
worthwhile. For example, one of the fundamental scientific pro- t ics of solids, liquids, 
grams being undertaken is a study of the mechanism of solidifica *% gases, and plasmas. * 
tion of polyphase alloys. ¥ 
This study is still in its early stages. However, we now have a 3% Currentstudiesrange 
better understanding of the solidification process. And we have ¥ from the fundamental {= 
demonstrated that the microstructure of certain eutectics can be * 
significantly altered by controlled solidification. The contrast is % Properties of matter + 
illustrated above. % to the application of £ 
It is possible that these binary alloys with preferred orientation - scientific knowledge f 
may have unusual physical and mechanical properties. x Pha & 
Perhaps you are interested in corporate—sponsored studies into :- to promising new ‘ 
the fundamental nature of matter. If so, we can offer a research $ products. x 
environment that seems to stimulate scientists to extra achieve- x - 
ment. It may be due to the way we encourage an internal cross- ¥ owe * z 
fertilization of ideas; or because of a unique complex of comple- t * 
mentary services that free the scientist from routine analysis and x For more specific information, x 
experimental chores. % please write to Mr. W. F. & 
Walsh, at the Research Lab- 
oratories. 
* 
RESEARCH LABORATORIES ies 
400 Main Street, East Hartford 8, Conn. 
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AERO-THERMODYNAMICS AND 
FLOW IN TURBOMACHINES 


By Micuazt H. Vavra, U. S. Naval Postgraduate 
School. Designed to give an understanding of 
flow phenomena and to show rational methods 
for formulating them mathematically. This is 
done by using a theoretical approach stressing 
basic physical processes rather than empirical 
methods. In combining this approach with his 
own experience the reader will gain an idea of de- 
sign possibilities. 

Part I sets up basic flow relations applying 
to any branch of fluid mechanics. Part II covers 
types more closely related to flow in machines, 
and Part III is directly concerned with flow in 
machines. One of the book’s most valuable fea- 
tures is its use of vector- and tensor-analysis to 
establish fundamental flow relations. 1960. 609 
pages. $14.50 
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POWERED FLIGHT 


By R. L. Carrot, Bell Aircraft Corp. A simple, 
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tion is made at this point to questions of powered 
flight with aerodynamic lift. 1960. 275 pages. 
$8.50. 
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Edited by Maj. Gen. Oris O. BENnson, JR., and 
Dr. Husertus StrucHoip, both of USAF Aero- 
space Medical Center, Brooks Air Force Base. A 
symposium sponsored by the School of Aviation 
Medicine and arranged by Southwest Research 
Institute. Stresses the environment of man in 
space—extreme conditions of space flight and 
effects of gravitational, radiational, and magnetic 
fields. Also: space vehicles, weightlessness, and 
other topics. 1960. 645 pages. $12.50. 7 


SURFACE EFFECTS ON 
SPACECRAFT MATERIALS 


Edited by F. J. Crauss, Lockheed Aircraft Corp. 
Record of a first symposium, covering methods 
for calculating required characteristics, avail- 
able materials, laboratory techniques for meas- 
uring radiation characteristics, and effects of 
space environment on behavior of materials. 
1960. 404 pages. $11.50. 


THEORY OF THERMAL STRESSES 


By B. A. Botey and J. H. Weiner, Columbia Uni- 
versity. Unites thermodynamical foundations 
and practical material of interest to engineers 
and researchers. Both thermal stress and theory 
of heat conduction are covered. 
590 pages. Prob. $15.50 
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Scope of ARS JOURNAL 


This Journal is devoted to the advance- 
ment of astronautics through the dissemina- 
tion of original papers disclosing new scientific 
knowledge and basic applications of such 
knowledge. The sciences of astronautics are 
understood here to embrace selected aspects 
of jet and rocket propulsion, spaceflight 
mechanics, high speed aerodynamics, flight 
guidance, space communications, atmospheric 
and outer space physics, materials and struc- 
tures, human engineering, overall system 
analysis, and possibly certain other scientific 
areas. The selection of papers to be printed 
will be governed by the pertinence of the topic 
to the field of astronautics, by the current or 
probable future significance of the research, 
and by the importance of distributing the in- 
formation to the members of the Society and 
to the profession at large. 


Information for Authors 


\lanuscripts must be as brief as the proper 
presentation of the ideas will allow. Ex- 
clusion of dispensable material and concise- 
nes- of expression will influence the Editors’ 
acceptance of a manuscript. In terms of 
standard-size double-spaced typed pages, a 
typical maximum length is 22 pages of text 
(including equations), 1 page of references, 
1 page of abstract and 12 illustrations. 
Fewer illustrations permit more text, and vice 
versa. Greater length will be acceptable 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for full 
articles, may qualify for publication as 
Technical Notes or Technical Comments. 
They may be devoted to new developments 
requiring prompt disclosure or to comments 
on previously published papers. Such manu- 
scripts are published within a few months of 
the date of receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
publication, according to the above-stated 
requirements as to subject, scope or length, 
but which nevertheless deserves widespread 
distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as a 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be double spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ posi- 
tions and affiliations in a footnote on the 
first page. Equations and symbols may be 
handwritten or typewritten; clarity for the 
printer is essential. Greek letters and unusual 
symbols should be identified in the margin. If 
handwritten, distinguish between capital and 
lower case letters, and indicate subscripts and 
superscripts. References are to be grouped at 
the end of the manuscript and are to be given 
as follows. For journal articles: Authors 
first, then title, journal, volume, year, page 
numbers; for books: Authors first, then title, 
publisher, city, edition and page or chapter 
numbers. Line drawings must be clear and 
sharp to make clear engravings. Use black 
ink on white paper or tracing cloth. Lettering 
should be large enough to be legible after 
reduction. Photographs should be glossy 
prints, not matte or semi-matte. Each illus- 
tration must have a legend; legends should be 
listed in order on a separate sheet. 

Manuscripts must be accompanied by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
nent sponsorship. Full responsibility rests 

th the author. 

Preprints of papers presented at ARS 

eetings are automatically considered for 
ublication. 

Submit manuscripts in duplicate (origi- 
nal plus first carbon, with two sets o 
illustrations) to the Managing Editor, ARS 
JOURNAL, 500 Fifth Avenue, New York 
36, N.Y. 

ARS JOURNAL is published monthly by 
the American Rocket Society, Inc. and the 
American Interplanetary Society at 20th & 
Northampton Sts., Easton, Pa., U. S. A. 
Editorial offices: 500 Fifth Ave., New York 

6, N. Y. Price: $12.50 per year, $2.00 per 
single copy. Second-class mail privileges 

thorized at Easton, Pa. This publication 

s authorized to be mailed at the special 
rates of postage prescribed by Section 132.122. 
Notice of change of address should be sent 

the Secretary, ARS, at least 30 days prior 
to publication. Opinions expressed herein 
are the authors’ and do not necessarily reflect 
the views of the Editors or of the Society. 
© Copyright 1960 by the American Rocket 
Society, Inc. 
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Scientific objectivity characterizes the 
examination of natural forces in the 
experimental laboratories at Los Alamos. 
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PROGRAMS SCOUT 


Builders of more large, thin wall, high 
strength solid propellant rocket engine 
cases and nozzles for development pur- 
poses than any other company in America. 


A small experienced organization geared 
to handle your development and prototype 
requirements for static and flight tests in 


the shortest possible time. 
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the many advantages of 


Chlorine Trifluoride (CIF3) 
as a rocket fuel oxidizer! 


Desirability 
—Chlorine Trifluoride is one of the most 
powerful rocket fuel oxidizers known. Its 
performance compares favorably with 
other oxidizers far more difficult to 
handles 
Properties 
-—CTF has many characteristics desirable 
4 a storable propellant, with a wide liquid 
Com : temperature range. It has a low freezing 
Patibini, (—105°F) and a moderate boiling 
y point (53°F). Vapor pressure is 80 psia 
140°F. 


Storability 
—Unlike most liquid propellants, Chlorine 
rifluoride requires no refrigeration. It is 
handled and stored safely and conveni- 
ently in standard steel cylinders. No spe- 
_ cial storing precautions necessary. 


Availability 
—CTF has been produced commercially 


by General Chemical for almost 10 years. 
Readily available in ton and 150-Ib. steel 
cylinders and in smaller sizes. 
Stability 
1 —CTF has very high shock resistance; is 
i shipped nationally via common carriers. 
; It is thermally stable up to 450°F. 


Dependability 


fuels over a wide range of pressures and 
temperatures. 


deal storable 


Flexibility 

—CTF can be used under all geographic 
and climatic conditions. No need for fixed 
supply points. 


Compatibility 

—CTF can be used satisfactorily with 
properly passivated steels, nickel alloys, 
copper and other common metals. 


Density 
—CTF’s high density (1.825 gm/cc at 
68°F or approximately 15.3 Ibs./gallon) 
leads to outstanding density impulse val- 
ues with a variety of rocket fuels. 


Write today for a free copy of our techni- 
cal bulletin, ‘Chlorine Trifluoride.”” Com- 
pany letterhead, please. 


llied 
GENERAL CHEMICAL DIVISION 
hem ical 40 Rector Street, New York 6, N. Y. a 
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DIMAZINE has many characteristics that 
make it excel as a storable fuel. 


Its remarkable stability has been proven in 
tests of more than two years’ duration, even 
in mild-steel vessels. Its high thermal sta- 
bility and low freezing point permit wide 
variations of storage temperatures. It is not 
shock sensitive and is unaffected by many 
types of accidental contamination. It is high- 


Putting 


FOOD MACHINERY 
AND CHEMICAL 
CORPORATION 


ill STORABILITY 


Just one of the many outstanding attributes of 


DIMAZINE 


unsym-Dimethylhydrazine 


ldeas to 


? 


TEXTS 


ly compatible with most metals and appro- 
priate sealant materials. 

Excellent storability, high performance, ease 
of handling, ample availability . . . these 
and many other advantages are combined in 
DIMAZINE, the outstanding storable fuel. 
We will be pleased to assist in your evalu- 
ation of DIMAZINE and to supply detailed 
data on its properties and handling. 
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FOOD MACHINERY AND CHEMICAL CORPORATION 
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Physics 


HIS paper summarizes the results and analysis of seven 

rocket experiments designed and performed to system- 
atically study the physics of the generation of artificial elec- 
tron clouds in the upper atmosphere. 

‘he several experiments can be discussed in four phases. 
For example, the first involved the daytime Aerobee rocket 
release of potassium at 121 km and served to establish the 
feasibility of solar photo-ionization generation of artificial 
electron clouds. The second phase was the Nike-Cajun 
nighttime release of cesium at 101 km to establish the role 
of thermochemical generation of electron clouds. The third 
phase was designed to study the correlation of optical and 
radio-radar data as obtained from a visible, photographic 
electron cloud. This was accomplished by the morning twi- 
light releases of cesium and sodium at 128 and 116 km, re- 
spectively. The final phase of this series reported here in- 
volved three releases of cesium at 91, 82 and 69 km, in order to 
determine the altitude below which long-lived, dense artificial 
electron clouds cannot be generated by the techniques em- 
ployed here. The results of this last phase yield valuable 
physical data pertinent to several of the electron decay proc- 
esses dominant in the lower altitudes. 

The overall analysis of the extensive optical and radar data 
has yielded information on the following physical parameters: 


1 Chemical yield of contaminant. 
2 Thermal ionization efficiency. 


Received April 26, 1960. 

! The original experimental work was performed while authors 
were associated with the Geophysics Research Directorate of the 
Air Foree Cambridge Research Center. Most of the data analy- 
sis and interpretation presented in this survey was accomplished 
by the authors at their present address under the sponsorship of 
contracts: U.S. Army Signal Corps DA 36-039-SC-78971; U.S. 
Army Signal Corps DA 36-039-SC-78925. 


F. F. MARMO, 
L. M. ASCHENBRAND and 
J. PRESSMAN 


Geophysics Corp. of America 


Upper atmospheric wind velocity. 


Optical and gaussian half-width of cloud with time. 
Ambipolar diffusion. 

Neutral diffusion. 

Photo-ionization cross sections. 

Mutual neutralization processes. 

10 Rate of chemical consumption of alkali by atmospheric 
oxygen. 

Finally, the role that artificial electron clouds can assume 
for practical utilization for radio-radar propagation is em- 
phasized with some discussion indicating the extent of the 
available potential for this engineering purpose. _ 


Introduction 


Artificial electron clouds are formed by the ejection of a 
suitable chemical contaminant at pre-selected altitudes with 
the generation of the electrons due to photo-ionization and/or 
thermal ionization. The first successful experiment on the 
generation of an artificial electron cloud was performed by an 
Aerobee rocket release of about 18 lb (275 moles) of nitric 
oxide gas at an altitude of 95km (1).2— The success of this ex- 
periment, whereby an artificial electron cloud with an initial 
average density of about 4 X 10° electrons per cm* was de- 
tected for more than 10 min by radar techniques, indicated 
that it would be feasible to use other contaminants to create 
artificial electron clouds at various altitudes with a wide range 
of densities, time durations, sizes, shapes and other param- 
eters. 

Since May 1957, a series of seven additional upper at- 
mosphere rocket experiments has been performed to examine 


2 Numbers in parentheses indicate References at end of paper. 
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Dr. Frederick F. Marmo studied physical chemistry at Boston University, where he received the | 
B.S. degree, Phi Beta Kappa and graduated summa cum laude in 1949. His M.S. and Ph.D. de- aa 
grees in chemistry were received at Harvard University in 1951 and 1954, respectively. 

Mr. Leonard M. Aschenbrand received the B.S. degree in physical chemistry from CCNY in 1947. 


. - Mr. Jerome Pressman received the B.S. degree in mathematics from CCNY in 1941. He later 4 
attended the University of Chicago and the University of Puerto Rico while taking advanced © 
- courses in meteorology and physics at NYU, 1946-1949. 


From 1951 through 1958 all three scientists were involved in theoretical and experimental studies 

relative to the study of upper atmospheric physical processes while employed in the Photochemistry 

a Laboratory at the Geophysics Research Directorate of the Air Force Cambridge Research Center. 

: u Dr. Marmo was chief of the Chemical Physics Branch, Mr. Aschenbrand of the Atmospheric Energy © 

4 Section and Mr. Pressman of the Theoretical Photochemistry Section. ie 

a All three were responsible for initiating, planning, conducting and eventually publishing papers — 

; and reports on the program of the creation of artificial electron clouds and synthetic chemilumines- 
cence in the upper atmosphere, using techniques of releasing chemical constituents. 
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- In 1949 he was awarded an M.A. degree in chemistry from NYU. 
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these possibilities, and to ascertain the feasibility of employing 
these clouds as a research tool for upper atmosphere re- 
search. In addition, the possibilities of employing artificial 
electron clouds for propagation purposes was explored. These 
initial experiments established the foundation of a new tech- 
nique of study of physical and chemical processes associated 
with Earth’s upper atmosphere. A substantial amount of 
qualitative and quantitative data has been obtained from 
these experiments pertinent to the studies of upper at- 
mospheric winds, diffusion, microscopic cross-section values 
(electron attachment, ion-electron recombination, charge 
exchange, etc.), chemical consumption, expansion of high 
pressure gases into an attenuated atmosphere, electromag- 
netic propagation research and implications, ionospheric re- 
search and thermochemical and explosive effects, to name 
Ve 


Fig. 1 Ionogram obtained 30 min after chemical release as re- 
corded by C-3 vertical ionosonde located at White Sands Proving 
Ground, N. Mex., for experiment conducted in November 1957. 
This display depicts the slant range vs. maximum frequency. 
Each altitude marker represents 100 km, whereas the frequency 
scale goes from 1 to 25 mcpersec. The artificial cloud is shown 
between the normal E and F regions at about 170-km slant 


range. The maximum frequency return is about 8 mc per sec 


but several. The purpose of this survey article is to briefly 
present the basic concepts employed, a chronological summary 
of experimental events, the data observed and the results of 
analysis of these data. 

The first successful experiment designed for optimum-type 
artificial electron clouds was accomplished by the release of 
atomic K from an Aerobee rocket at 121 km. This flight oc- 
curred at about 0932 MST on Nov. 19, 1957. The electron 
cloud formed as a result of this experiment was readily de- 
tected by conventional radar techniques. Fig. 1 depicts one 
of the more interesting representative pieces of radar data. 
This result is an ionogram obtained from a C-3 vertical iono- 
sonde, located at White Sands Proving Grounds, N. Mex. 
which clearly shows the detection of the manmade electron 
cloud formed between the normal E and F layers (at about 
170-km slant range) at 30 min after contaminant release. 

Interpretation of the results pertinent to the ionized cloud 
which was formed in this experiment suggested that some ad- 
ditional process to photo-ionization was responsible for thie 
initial electron production. The chemical nature of the ex- 
perimental method indicated that the effect could be due to 
thermal ionization. In order to examine this important cf- 
feet, an experiment was devised in which atomic Cs would |e 
released during the nighttime hours to eliminate the role of 
solar photo-ionizing flux. The choice of this element arose 
from the consideration that, at a given temperature within 
the range considered, Cs should produce considerably more 
thermal electrons than other constituents of the alkali metal 
family; thus, a maximum-type effect could be anticipated. 
In addition, since this was to be a nighttime release, it was 
considered feasible to produce a cloud easily visible and, 
therefore, easily photographed. 

Some of the obvious advantages gained by optical detec- 
tion are: Accurate determination of release altitudes by tri- 
angulation optics, correlation of radar returns as a function 
of cloud size and shape, determination of yield and chemical 
rates, wind and diffusion data, and spectral characteristics of 
the emission. The method employed for rendering the cloud 
visible was to ‘‘dope” a Cs 100-km release with about 10 per 
cent Na so that the latter would react in a chemiluminescent 
manner with the ambient atomic oxygen in the upper at- 
mosphere (maximum number density at about 100 km). Of 
equal importance was the fact that results of the November 
17 Aerobee flight suggested the desirability of utilizing the 
smaller Nike-Cajun two-stage rocket. Known advantages 
of using the smaller rocket included lower cost, ease of han- 
dling at launching, solid fuel and superior trajectory charac- 
teristics. However, it was recognized that, for the production 
of maximum-type electron clouds, such utilization would incur 
the sacrifice of larger payload. (The Aerobee’s instrumented 
payload is capable of carrying aloft about three times more 
constituent.) In spite of this, the decision was made to utilize 
the Nike-Cajun vehicle whenever possible. 

The successful release of this Cs-Na mixture was accon- 
plished from a Nike-Cajun rocket at an altitude of 101 km at 
the pre-dawn hour of 0245 MST on May 20, 1958. The elec- 
tron cloud created was detected by radar; the cloud was as- 
sumed to be due to thermal ionization. It was clearly visible 
to the naked eye; the photographic and spectral data were ob- 
tained for as long as 6 min. Interpretation of the optics! 
data presents the possibility of a study of the cloud dynamic= 
vs. time, as well as information pertaining to chemical con- 
sumption, kinetics of the reactions involved and some of th: 
physical processes previously mentioned. 

The next group of experiments of the series was primarily 
concerned with the correlation of optical and radar data from 
ionized clouds as formed by solar ionization of alkali meta! 
vapor. This could conveniently be accomplished by creating 
a light emitting photographic cloud, preferably one that 
emitted in the visible. The method of rendering the 
cloud visible is based upon the detectability of Na 
resonance radiation against a morning-twilight-sky back- 
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ground (2). Again, if we merely “doped” the chemical pay- 
load of Cs compound with a small amount of Na, an easily 
observable cloud would be formed. A morning-twilight 
release was adopted, since it was essential that the cloud be 
solar illuminated for an extended period, and that a reduced 
intensity sky background exist for the early part of the experi- 
meut. Two different experiments were performed as part of 
this phase of solar induced ionization. A release of a mixture 
of Cs and Na (2 per cent) was conducted at 128 km on May 
21, 1958, and, following that, a payload consisting entirely 
of Na compound was released at 116 km on May 22, 1958. 
Both experiments were successful in that they afforded the 
opportunity of providing good photographic and spectral 
coverage for over 8 min. Fig. 2 shows a series of three pic- 
tures obtained on the May 21 experiment using a light ampli- 
fying unit (Lumicon) capable of obtaining series photographs. 
The first picture was obtained 2 sec after release time, and the 
latter two show the cloud development after 1 and 2 min, 
respectively. 

The expansion characteristics are clearly indicated, and 
the utilization for data of this type is obvious (3). 

The radar results from the former release provided positive 
results for as long as 19 min, resulting in valuable correlated 
radar and optical data. The results of the Na release were 
uniform in that all operating stations reported negative radar 
results. This negative result appears to be in agreement with 
the deduced relative photo-ionization probabilities of Cs, K 
and Na, namely, 6.5 X 10-4, 1.8 X 10-5 and 1.0 X 107 
sec~!, respectively (4). This further indicates, as expected, 
that thermal ionization does not play an important role for 
Na as far as production of a detectable artificial electron cloud 
is concerned. 

The last series of experiments reported here consisted of the 
daytime Cs releases at 70, 80 and 90 km. This series was 
conducted during September 1958 to determine the altitude 
below which the generation of a detectable, high density, 
long-lived solar induced artificial electron cloud is precluded 
due to the rate of photo-ionization electron production, re- 
combination, electron attachment processes and chemical 
consumption of the contaminant. The basic engineering 
technique concerning these releases of alkali metal vapors in- 
volved a stainless steel canister so designed and constructed 
as to be within the requirements of payload, configuration of 
nose cone, high acceleration and, finally, adaptability to the 
rocket body and instrument package. No further details 
are given here, since much of the results reported here have 
been previously reported in a more extensive manner in a 
series of papers (3-8) dealing with the experimental data in 
detail, theoretical considerations, evaluation and interpreta- 
tion, ete. Finally, the extensive radar and optical probes 
utilized have also been described previously (3-8). 


JuNE 1960 


Photographs taken from television screen of Lumicon depicting the dynamics of cloud formation for experiment conducted 
May 21, 1958. The picture (a) was obtained 2 sec after chemical release, (b) and (c) after 1 and 2 min, respectively 


Data Analysis and Interpretation 


In addition to their use as a powerful radar propagation 
tool, the scientific uses of electron clouds are varied; i.e., the 
scientist is given a valuable research tool for the investigation 
of physical and chemical processes of the upper atmosphere. 
Many of the scientific areas explored by means of electron 
cloud techniques are reported here. This work is restricted 


to discussions of the following: or 


Thermal ionization efficiencies. 
Chemical yield. 

Ambipolar and neutral diffusion coefficients. 
Horizontal wind velocities in the upper atmosphere. 
Wind shear. 

Chemical consumption of alkali constituent by atmos- 
pheric oxygen. 

7 Lower atmospheric electron decay processes. 

These do not exhaust the possible areas of investigation, but 
rather serve to illustrate the potential utility of artificial 
electron clouds as a research tool. The following sections are 
devoted to presenting the brief analysis on the experimental 
data associated with items 1 through 7. 


Analysis of Thermal Release: Chemical Yield and 
Thermal Ionization Efficiency 


A sodium-doped Cs payload was released at 101 km during 
a nighttime experiment to determine the magnitude of ther- 
mochemical effects upon the production of artificial electron 
clouds. Subsequent analysis of radio-radar and optical data 
pointed toward the fact that it was indeed possible to ascribe 
to thermal processes the generation of an ample number of 
electrons to account for the observed electron cloud. A 
thermochemical treatment of these data with application of 
the Saha equation suggested that values for the total number 
of electrons and thermal ionization probability could best be 
obtained from conditions wherein the temperature of reac- 
tion would be in the vicinity of 3000 K and with the canister 
bursting pressure range between 300 to 500 atm. 

Fig. 3 was derived from the Saha relationships as applied to 
the specific system employed. It relates the moles reacted to 
the temperature of the reactants for different bursting pres- 
sures. With the aid of this figure it can be shown (6) that a 
chemical yield of less than 10 per cent was to be expected. 
This fact was further substantiated by optical and radar data 
from subsequent experiments in the series in which yields of 
between 2 to 5 per cent had been determined in an independ- 
ent manner. 

An initial cloud optical radius of about 300 m and an initial 
electron density of the order of 3 X 107 em~* were determined 
from an analysis of optical and radar data, respectively. 
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From these determinations it was experimentally concluded 
that the total electron content of the cloud N, was of the 
order of 3 X 1021. Theoretical evaluation of VN, = 1 X 10?! 
was determined from subsequent Cs releases; a discussion 
will be found in another section of this paper. An N, value 
of between 1 to 3 X 10?! was adopted. 

For Cs a value of 10~* for the thermal efficiency c was de- 
termined from the results of this expériment in thermal release. 
Further, through utilization of the Cs data it was possible to 
extend these findings to other members of the metal alkali 
family. It was thus, that a corresponding value of 5 X 10~¢ 
was estimated for potassium. These c values are employed 
in later data analysis. 
Analysis of C-3 Radar Data: Diffusion, Chemical Yield, 
Ionization Efficiency 


This section will outline some of the methods which were 
employed to analyze radar data. The results of this analysis 
yielded values for several important parameters. 

Of all the radar data obtained, the C-3 ionosonde data were 
considered the most compatible to analysis. The analysis 
is restricted to the K and Cs releases at 121 and 128 km, 
respectively. It has previously been demonstrated (5,7) 


that for these cases, the significant factor causing a decrease 
in the electron concentration is diffusion. 


Decay processes, 


such as attachment, recombination and others, are not suf- 
ficiently fast for the time constants of the effects under discus- 


sion to materially affect electron concentration. 


Thus, con- 


siderable simplification is possible; in particular, the 
governing equations (7) are given by f 

= D,V*n — k {1 

= + kn [2} 


KN + Ne. — no 


nm = concentration of neutrals, em* 

Ne = concentration of electrons, em* 

D, = neutral diffusion coefficient, em?/sec 
D, = ambipolar diffusion coefficient, cm?/sec 
k = photo-ionization probability, 


The solution of Equations [1 and 2] have been given e se- 
where (7). The center point density solution can be writ!en 
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Fig. 3 Pressure within canister as a function of moles reacted 
and temperature of the reactants 
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Fig. 4 Plot of the C-3 experimental data for the release of ce- 
sium at 128 km, depicting the best fit values of the ambipolar 
and neutral diffusion coefficients D, and D, respectively. The 
two closely spaced curves serve to illustrate the sensitivity of the 
technique 
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1 exp (—kt 
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kp 
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where 
n= (ro? + 4D,l) 
Pp (70? + 4D,t) 
Ne = initial electron center density, which includes the initial 


thermal and the photo-ionization contributions 


It may appear that a large degree of flexibility is available 
for proper fitting in view of the several parameters involved. 
However, this is probably not the case; in fact, the require- 
ment of utilizing several reasonable values simultaneously for 
proper fit can be highly restrictive. For example, it was 
shown earlier that approximate N values for the cases of the 
Nike-Cajun Cs and Aerobee K releases should be about 1 X 
1074 and 3 X 1074, respectively. Further, it was ascertained 
that thermal ionization efficiencies (c values) of Cs and K are 
about 10-*and 5 X 10~4, respectively, for the chemical release 
system employed here. These values should be utilized for 
proper fit of the data to the theoretical curves. 


POTASSIUM RELEASE AT 121 KM 
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Fig. 5 Plot of the C-3 experimental data for the release of po- 


tassium at 121 km, depicting the best fit values of the ambipolar 
and neutral diffusion coefficients 


ARS Journal 


in 
rele: 
the: 
coell 
em- 
121- 
10° 

TI 
indic 
radal 
D,, v: 
nique 


Opti 


Fo 
it is a 
oecur 
1%. | 
densi 
ance | 


where 


D 


Iflum 
in the 
tered 

section 
the lin 
atom ; 
hence, 
¢] 
sential! 
sional « 
the clo 
wave | 
remain 
ambien 
consur 
sec at 1 


JUNE ] 


— 
wi 
I 
an 
we 
un 
pr 
not 
( 
pro 
dat 
reli 
sub 
. 
: tio! 
agr 
D, 
T(°K) 
1000 
2000 
1000 
100 Atmos. 
| 
= > 


For the estimation of the ro values to be employed, optical 
width data were employed from the Nike-Cajun Cs and Na 
twilight releases at altitudes of 128 and 116 km, respectively 
(7). These data were then theoretically adjusted to the corre- 
sponding gaussian half-widths rm. The 7, values of 300 
and 200 m were deduced for Cs and K, respectively. These 
were applied directly for data interpretation. Although some 
uncertainty exists in the use of this method for estimating the 
proper gaussian half-width, this parameter is fortunately 
not too sensitive to fitting. 

Concerning the choices of k values, solar photo-ionization 
probability values of 6.6 K 10~4see~! for Cs and 1.8 X 10% 
sec ' for K were derived from published experimental 
data (4,9). These values were considered to be the most 
reli:ble of the several parameters employed and, as such, not 
subject to variation, thereby imposing an additional restric- 
tio! upon the available flexibility for proper curve fitting. 
There now remained the problem of determining the degree of 
agr-ement achievable between the theoretical and experimen- 
tal ‘alues; application of various reasonable values of D, and 
D,, the ambipolar and neutral diffusion coefficients, respec- 
tiv: ly, were used toward this end. Indeed, it was found that 
in both eases, D, and D, values and the ratio of D,/D, were 
in close agreement with theory (10). For the 128-km Cs 
release, the best fit values were: yield, 1 Xx 1074 atoms; 
thermal efficiency, 10~%; neutral and ambipolar diffusion 
coefficients D, and D,, 2.5 X 108 em?/sec and 4.2 X 107 
em® ‘see, respectively. The corresponding values for the 
121-km K release were 3 X 1074 atoms, 5 X 1074 and 6.2 X 
10’ em2/sec and 9.2 X 106 cm?/sec. 

The resultant fit to experimental data using these values is 
indicated in Figs. 4 and 5, which include the appropriate C-3 
radar experimental data. In these figures two closely spaced 
D,, values are shown to illustrate the sensitivity of the tech- 
nique. 


Optical Data: Diffusion and Chemical Yield 


For obtaining diffusion and chemical yield determinations, 
it is assumed that the initial, very rapid, explosive expansion 
occurring at release leads to a gaussian cloud of half-width 

It is at this point that we set ¢ = 0. Subsequently, the 
density is expected to vary with time and position in accord- 
ance with the equation 


N r? 
4 
(r, 4Dt + ( 4Dt + 3) 

where 
n(r,t) = number of particles per cm’ at distance r from the 
center at time / 
N = total number of particles in cloud 

tT = half-width of gaussian (1/e value), em 

D = diffusion coefficient for the species at cloud height 


(em?/see ) 


If luminosity is due to resonance scattered sunlight and viewed 
in the backscattering direction, light will, therefore, be scat- 
tered to the detector as a function of the absorption cross 
section a, and in accordance with 7 the number of atoms in 
the line of sight. The absorption cross section ¢, of the Na 
atom and its resonance level at 5890 A is about 10~' cm?; 
henee, when the number of atoms in the line of sight is about 
10'', this portion of the cloud may be considered as being es- 
sentially optically dense. Since at these altitudes the colli- 
sional deactivation is negligible, the optically dense portion of 
the cloud will have a reflectivity of unity at the characteristic 
wave length. Further, the total number of alkali atoms will 
remain constant, since the rate of chemical consumption by 
ambient oxygen is so slow. For example, the time for half- 
consumption ti. is greater than 10’ sec and is greater than 10° 
sec at 130 and 115 km, respectively. [See (4), fig. 6.] 
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The total number of particles in the line of sight (7) at dis- 
tance d (em) from center is given by 


r? 

= exp ( 5 

+ P ( 4Dt + =) (5) 

If it is assumed that the optical size of the cloud is deter- 

mined by the distance from the center of that isophote for 

which = (since the intensity drops rapidly for > o), 
we have 


—(4Dt + ro?) In — + 1?) [6] 


d? = 


= |7 ] 
eT 
ta [8] 
Ver 


1 (No 9] 
= — ¢ 
4D Ne : 


A 
where 
a 

To effectively relate the optical model to the observed or 
photographed cloud size with the contaminant density and 
profile throughout the cloud, requires consideration of the de- 
tector response to cloud luminosity under the varying twi- 
light levels of background light. Many factors unite to make 
quantitative measurements from a photographic image dif- 
ficult. ‘Thus, for quantitative work it is necessary, in prac- 
tice, to combine careful instrument calibration with meticu- 
lous densitometry techniques to plot contour lines of constant 
density and relate these to a theoretical model before dimen- 
sions can be accurately ascertained. This was not accom- 
plished in the present work, and the absolute values derived 
suffered accordingly. 

As an illustrative example of the crude optical model em- 
ployed here, we analyze the Beattie-Coleman camera data 
for the May 21 flight (7). Here we are concerned with de- 
termining the parameters: 7, the initial gaussian half-width; 
N, the number of alkali atoms; and in this case D,, the ob- 


dmax = Maximum radius of optical radius 
tang, = time to achieve this maximum radius 
tay = time for cloud to disappear 
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Fig. 6 The optical determination of the number of atoms N and 
neutral diffusion coefficient D,. The time variation of the 
minor axis optical diameter is indicated by the circles. The 


dashed and solid curves are theoretical curves for two different 
yields but the same diffusion coefficients _ 
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served neutral diffusion coefficient. From Equation [6] for 


t = 0, we obtain 
dy? = In [10] 


where dy is the initial optical radius. 

Utilizing Equations [8 and 10] we obtain the fit illustrated 
in Fig. 6 where the size (2d) is plotted against the theoretical 
curve for two closely spaced, different values of N. This 
serves to show that the fit is surprisingly accurate, since small 
variations in both N and D, cause large changes in the curve. 
The data indicate a diffusion coefficient value of D, = 7 X 
10’ cm2/sec, and also that N is equal to 1.5 X 10? atoms (ora 
sodium yield of approximately 3 per cent). The value of D, 
is consistent with the theoretically extrapolated values of D, 
at the altidude in question (10). Furthermore, since tmax = 
1/D,, small changes in D, shift the maximum and the curve 
along the ¢ axis in an extremely sensitive fashion. ap: 


Wind Determination: Velocity and Shear Data a 


The formation of a small visible photographic cloud at alti- 
tude can yield wind velocity data by means of simple optical 
triangulation techniques. A method of obtaining shear in- 
formation from optical data is given here. In addition, a 
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Fig. 7 Upper atmospheric horizontal wind velocity determin- 
ation. The circles represent the C-3 radar slant range vs. time 
forthe May 21 experiment. The solid line is a best fit from Equa- 
tion [11] and a constant, horizontal wind velocity of 140 m per sec 
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Fig. 8 Schematic representation of geometrical determination 
of winds from C-3 slant range vs. time plot 


Equation [12]. 


geometric method of extracting accurate horizontal wind 


velocity values from C-3 radar data is discussed. 


The pertinent radar data for obtaining horizontal wind 
velocities is the change in electron cloud slant range vs. time. 
The data from the C-3 ionosonde are shown as circles in Fig. 
7. The smooth variation of the range with time for this 20- 
min period is strong support for the presence of an isolate 
electron cloud. From the initial slant range of about 150 
km, the subsequent change in slant range in time is presun- 
The character of the curve 
with a shallow minimum indicates the geometry depicted \n 
Fig. 8, in which the positions of cloud release and C-3 statiois 
are included. This geometry can be described mathemati- 


ably due to ionospheric winds. 


cally by the expression 


ws 
S = + wi(wt — 28,cos%)}* 


where 


So = initial slant range 


S = slant range between cloud trajectory and C-3 line >f 
sight 

6 = initial angle 

w = wind velocity 


A best fit to the data was obtained for a w value of 140 m p:r 
sec; the theoretical curve is shown as the solid line in Fig. 7. 
The excellent fit between the experimental points and tlie 
theoretical curve appeared to corroborate the notion of hori- 
zontal winds with a constant velocity of about 140 m per see 
This (140-m per sec) value for the velocity is considerably 
greater than those previously determined in a similar manner, 
namely, 60 and 78 m per sec at 95 km (1) and 121 km (5), 
respectively. Concerning wind direction, it is clear from Fig. 
8 that a direction ambiguity exists on the basis of these data 
alone. However, in this instance visual observations dictated 
the choice of direction as 87 deg rather than 293 deg. This 
case suffices to illustrate the geometric method of obtaining 
wind information via these and similar data. 

Wind shear in the upper atmosphere has recently been 
measured by a technique in which a long sodium filament is 
ejected into the upper atmosphere, between 70 and 200 
km (11). Careful analysis of the optical data leads to accu- 
rate shear data. Clearly, in a similar manner a point release 
can be utilized to determine the wind shear at the release alti- 
tude. For example, photographs of the contaminant cloud 
indicated an initially spherical cloud becoming ellipsoidal 
with time due to the prevailing wind shear. 

The photographed image of a cloud illuminated by reso- 
nance scattering is a complicated function of viewing angle, 
optical density of the cloud, film sensitivity and exposure time 
as well as background light intensity. Thus, for the experi- 
mental determination of shear by the employment of these 
clouds, a simplified model has been previously derived (7) 
and leads to the expression 


L = Y{at + (1/2)bt?] + constant [12] 
where 
a,b = constants of linearity of the shear 
L = elongation, em 
Y = shear coefficient, sec™ 


_ Camera data from the May 22, 1958 flight are plotted in Fig. 9. 
_ The variation with time of the minor axis is shown on thie 


left-hand graph, and on the right-hand graph is shown the 
variation with time of the major axis together with the plot 
From these, a shear value of Y = 1.1 X 
10~? sec was obtained for the 116-km release. 

Similar plots have been applied to Lumicon data (7), for 


which data from the May 22 flight yielded a shear constant 


of 1 X 10~* see in agreement with the foregoing. The 
value obtained for the May 21 Cs release at 128 km was 6 X 
10-* sec™!. These values were found to be consistent with 
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those obtained by sodium streamer experiments (11). It is 
recognized that the rate of diffusion will be a function of the 
cloud shape itself, but this effect is ignored in the given model. 
Another assumption is that shear does not affect the growth 
of the minor axis. Yet despite these limitations, a reasonable 
plot was obtained for the concomitant variation of the major 
and minor axes. 


Low Altitude Data Analysis; Chemical Consumption, 
Mutual Neutralization and Other Electron Decay Processes 


A series of three experiments will be discussed here, which 
were undertaken for the twofold purpose of studying the ef- 
fect of low altitudes (below 100 km) upon the formation of 
artificial electron clouds and the roles of a number of electron 
decay processes after the generation of these clouds. It was 
felt that the data acquired would afford not only an investiga- 
tion of the effects of individual decay processes, but would al- 
ow further investigation of an earlier conclusion (4) that there 
probably existed some lower altitude below which the forma- 
tion of detectable, high density, long lived and solar induced 
artificial electron clouds could not be achieved within the 
framework of the experimental techniques employed in this 


More pertinent to this conclusion, however, is the 


fact that several electron decay processes (e.g., chemical con- 
sumption, recombination, electron attachment, etc.) predom- 
inate at low altitudes in contradistinction to other processes 
which exert maximum influence above the 100-km region. 

For the present analysis, therefore, information was ob- 
tained concerning some of these decay processes. Of particu- 
lar inport was the derivation of coefficients for chemical con- 
sumption and mutual recombination. 

Of the several radio-radar probes utilized in the Cs releases 
at 91, 82 and 69 km (8), a 40.7 pulsed radar located 15 miles 
east of the launch site yielded the most consistent—and there- 
fore usable—data. The analysis applied to these data, indi- 
cated a total cloud duration (for 40.7 me echoes) of about 70 
sec at 91 km, 4.3 sec at 82 km and 1.5 sec at 69 km. 

The short radar reflecting lifetimes observed for the low 
altitude generated clouds are due to at least two important ef- 
The rapid destruction of Cs atoms by chemical con- 
sumption which prevents subsequent photo-ionization, and 
the rapid electron destruction through attachment followed 
by mutual neutralization. It was found (for the types of 
releases discussed in this paper) that the diffusion rate of the 
charged particles is much too slow to account for the observed 
cloud history. Also, the photo-ionization production term 
remained, relatively, not influenced by diffusion of the neu- 
trals, except at about 90 km. 

We present here a brief summary of the mathematical 


Details have been given in a previous paper (4). We 


dn,/dt = k no’ exp (—bt) — an? [13] 


= Ne 

= chemical consumption rate of consumption of the 
cesium by atmospheric oxygen, sec! 

= m/(1 + A) 


= ratio of negative ions to electrons at equilibrium — 


a=a,+ ha; [14] 


coefficient of mutual neutralization 
Cs-ion recombination coefficient 
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Fig. 9 Data from the May 22, 1958 flight. The solid dots and 

circles represent the time variation of the optical diameters of the 

minor and majoraxes, respectively. The twocurves show the best 
theoretical fit and indicate a wind shear of 1 X 10~? sec? 
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Fig. 10 Theoretical time variation of center point electron den- — 
sity in accordance with Equation [13] for 69, 82 and 91 km 


Fortunately, this nonlinear equation is susceptible to the 
complete formal solution (8) 


= 
bl, — + WE) + h(E) Ko) 
2a * + Lo(lE)Ki(l)] + — 


[15] 
where 
= (2/b)(2k 
BK = 2an(0)/b 
(1/2)bt 
= 
ye = 


Both expressions for y satisfy Bessels’ zero-order modified 
equation 
d?y 


1 dy 
—+4+— 16 
[16] 
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Unfortunately, the solution of Equation [15] is not susceptible 
to direct computation because of the large value of the argu- 
ments. However, numerical values have been obtained by 
the use of a variety of numerical approximations. These 
approximations and their range of validity are given else- 
where (8). 

Fig. 10 illustrates the fit obtained with the utilization of the 
best values for the several parameters involved. However, 
it is worth noting here that the deviation between the theoreti- 
cal and experimental data is greatest, as expected, at 91 km 
where no allowance was made for the effect of neutral dif- 
fusion on the photo-ionization term. Its direction is also as 
expected. Further, from the analysis it is deduced that a; = 
1077 — 10-* cm’ sec, in rough agreement with estimates of 
Bates (12) for similar processes. The value of 6 is in reason- 
able agreement with the values computed in an earlier paper 
(8) using the sodium rate of consumption. 

Using present techniques, the lowest altitude of effective 
cloud duration is about 70 km. A nighttime, low altitude 
experiment should give an independent measure of electron 
attachment rate of O. which can be meshed with the low alti- 
tude, daytime shots to obtain a more direct value of a;, the 
coefficient of mutual neutralization. 


Conclusions 


The foregoing is intended as a review of some of the main 
points in the systematic study of the physics of artificial elec- 
tron clouds. Its goal is the elucidation of the basic electron 
creation and destruction processes. 

This study has revealed that the theoretical photo-ioniza- 
tion cross sections are consistent with the data contained, and 
that cesium is the most effective of the alkali metals investi- 
gated for this purpose. It has been shown that the chemical 
yields have been disappointingly small, of the order of 2 to 3 
per cent. Thermal ionization has been shown to be an effec- 
tive process in the generation of artificial electron clouds, but, 
here again, thermal efficiencies of only about 0.1 per cent of 
the released number of atoms have been achieved. 

An optical method has been devised for the successful de- 
termination of the cloud initial half-width, chemical yield, 
wind shear and neutral diffusion coefficients. A simple 
method has been used for determining (with one degree of 
ambiguity) the wind velocity from radar slant range data. A 
mathematical model has been applied to those experimental 
releases above 120 km in order to describe the observed elec- 
tron decay history. The analysis gives measures of neutral 
and ambipolar diffusion, chemical yield, degree of thermal 
ionization and initial half-width. Reasonable correlation of 
radar data with the optical techniques employed has been 

demonstrated. 


ployed for the lower altitude (below 90 km) range which ac- 
counts for the dominant processes at these altitudes. These 
techniques also yield useful information on neutral and ambi- 
polar diffusion, photo-ionization cross sections, mutual neu- 
tralization, thermal ionization and the chemical consumption 
of the alkali contaminant by atmospheric oxygen. Infor- 
mation on electron attachment rates should be easily available 
by application of similar techniques. Additional work is in 
progress to solve the mathematical problems in the 90 to 
120-km region, where several processes are important. 
Finally, the foregoing summary has made obvious the po-- 
sible usage of artificial electron clouds for radar propagation 
purposes. In particular, a significant potential for this appli- 
cation is indicated by the small chemical yields and therm:| 
efficiencies resulting from the methods employed here. 
Clearly, these are capable of substantial increases. A system- 
atic chemical engineering study is being currently pursued 
to this end and should yield substantial improvement toward 
realizing the full potential of the utilization of artificial elec- 


tron clouds for radio-radar propagation (13). 
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FUNDAMENTAL equation governing the per- 

formance of an \-stage rocket in drag-free space is 


[1] 
where 
YY, [2] 
i=1 
Vi = + git [3] 


the latter being the dimensionless velocity increment of the ith 
stage, and 


N 
R= 
= i= 
where 
~ 


This is the mass ratio of the ith stage. These results are well 
known, having been derived, for example, in (1-4).? 

The purposes of this article are first to present the funda- 
mental performance equation in an especially useful form for 
the rocket systems designer and second to derive basic equa- 
tions showing the change in performance associated with 
changes in the characteristics of the individual stages. 

Referring to Fig. 1, adapted from (2), we see that the initial 
mass of any stage M)“ is simply the sum of the propellant 
mass M, and the empty mass MV,“ of that stage plus the 
initial mass Mo“*” of the next stage “| 


M. = M,® + + Mot 6] 
The burnout mass of any stage is 
Mo = M© + [7] 

and the propellant mass of any stage is the difference 
= — [8] 


Obviously the payload mass M,, is the initial mass of the 
(V + 1) stage 


M, = = [9] 
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The ideal performance of multistage rockets is presented in a graphical form especially suited for 
the systems designer interested in visualizing the intricate relationships between motor propellant 
fraction, size of the total propulsion system in relation to the payload, and burnout velocity. The 
effect of incremental changes in the propellant mass or inert mass of any stage and the effect of an 
| incremental change in the payload on the final burnout velocity are derived in the form of co- 
efficients that allow direct comparison of the relative importance of the characteristics of each 


and 
M,%t) = 0 [10] 
If the stage propellant fractien ¢; is defined by 
= 
so that 
1— = (12] 


then the mass ratio is given by 


Ri = 1/1 — [13] 
fi) 
e 
} | i+) 
= 
(i) 
M 
° 
wi 
p_4 
2227724 


M 
° % 


Fig. 1 Nomenclature for multistage rockets 


= 
(2) 
M 
p 
531 


_ Alternatively, we may define the motor propellant fraction 
A; by 
A; M,/(M, + [14] 


and write the stage ratio, i.e., the ratio of the total mass of 
any stage to the “payload” mass of that stage, as 


a; = [15] 


Then the ratio of the mass of the loaded motor of any stage to 
the ‘‘payload”’ mass it carries is 


a; — 1 = + [16] 


From Equations [14 and 16] it is apparent that the ratio of 
the propellant mass in any stage to the “payload’’ mass it 
accelerates is 

Ada; 1) = [17] 


and hence from Equations [11 and 15] the stage propellant 
fraction is given in terms of the motor propellant fraction by 

= Ala; — 1)/a; [18] 
Consequently, by Equation [13] the mass ratio for any stage 


can be written 7 
a,—1\]" 
R= A; (= aa [19] 


In the ideal optimum case where ¢;, A;, a:, R; and VU; are 
each constants unchanged from stage to stage, the total di- 
mensionless velocity obtainable with a given motor propel- 
lant fraction A; and a given stage ratio a; is 


Conversely, the stage ratio needed to achieve a given dimen- 
sionless velocity when the motor prorellant fraction is speci- 
fied can be found by solving Equation [20] for a; to give 


a; = Ai/ {Aj [1 — exp — (U/N)]}} [21] 


More important than the stage ratio is the total vehicle mass 
in relation to the payload mass 


a = M)/M, [22] 


=) 
° 


> 
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> 
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Fig. 2 Size of one-stage rocket vehicles as a function of the 
powerplant propellant fraction for various values of the dimen- 
_ sionless burnout velocity 
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d Vehicle in Relation to Payload, a 


which, from Equation [15], is ’ 


a=aN [23 


Hence from Equation [21] we get 


A; 
— [1 — exp — 


or, from Equation [20] 
VU = —N In [1 — Aa” — [25] 


Several interesting conclusions can be derived immediately 
from the last two equations. First, the minimum vehicle size 
in relation to the payload that is possible in the limit when 


1 is 


Conversely, the maximum velocity attainable with a given 
vehicle size for a given payload is 


VL [27] 


These limits are independent of the number of stages. It is 
also apparent from Equation [24] that a lower limit to the 
permissible value of the motor propellant fraction that can 
be used to attain a specified value of the stage velocity incre- 
ment U; = U/N is set by 


A; = 1 — exp — VU; [28] 


Finally, we may note that the limiting value of U for a given 
value of the motor propellant fraction A; as VN © is 


U=A,Ina [29] 


The most useful representation of the performance of N- 
stage rockets appears to be that afforded by Equation [24] 
which is plotted in Figs. 2-6 for N = 1, 2,3,4and5. These 
figures present, for a specified value of the total dimensionless 
velocity U, the total vehicle size a in relation to the payload, 
as a function of the parameter of most direct concern to the 
rocket designer, namely the motor propellant fraction <A,. 

Of considerable practical importance are the effects of small 


Loade 


Mass 


0.4 


Propellant Fraction of Power Plants, A; 


Fig. 3. Size of two-stage rocket vehicles as a function of the 
powerplant propellant fraction for various values of the dimen- 
sionless burnout velocity 
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Fig. 4 Size of three-stage rocket vehicles as a function of the 
powerplant propellant fraction for various values of the dimen- 


a sionless burnout velocity 


changes in the inert parts, propellant or payload masses on 
the overall vehicle performance. These effects can be deter- 
mined for the general N-stage rocket (where ¢;, Ai, ai, Ri 
and U; are not necessarily the same for each stage) as follows. 
Differentiation of Equations [1 and 4] gives 


= d&i/Bi [30] 


Thus the absolute change in the dimensionless velocity incre- 

ment of any stage is equal to the fractional change in its mass 

ratio. To put Equation [30] in terms of the actual masses 

involved, we start by differentiating Equation [5] to give 
dR; dM = dM 


MO” Me 


Utilizing dM) = dM,“ + dM, obtained from Equation 
[8] allows Equation [31] to be put in the form cee 


d&; Ee _ aM, 

We now recall that the burnout mass of the ith stage includes 
the propellant mass of all succeeding stages up to and includ- 
ing the Nth, the empty mass of the ith rocket stage and the 
empty mass of all succeeding rocket stages, including the 


(NV + 1) or payload a Hence, we may write — 


N+1 
b j=i [33] 
M, M,® 


In this equation terms of the form dM,/M,@ can be 
written 


— My dM,@ _ \ Mo dM, 
vhile 
MM 1 1 
— ..,— [35] 


My °° MyG-» 


r 


Therefore 
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‘Fig. 5 Size of four-stage rocket vehicles as a function of the 


powerplant propellant fraction for various values of the dimen- 
sionless burnout velocity 
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Fig. 6 Size of five-stage rocket vehicles as a function of the 
atti caeal propellant fraction for various values of the dimen- 
sionless burnout velocity 


J 6 
Similarly, terms of the form dM,“ /M, can be written 


dM) Mo 


( A; mo uo 
Therefore 


A; (2 (38] 
1, i l= Q; Ai+1 aj-1/ M® 
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With these transformations, it is now possible to write 


Equation [33] in the form q 


MY j=itl Ai+1 Qj-1 M,” j=i A; Ai+1 aj-1 


and Equation [32] becomes 


dR; dM 1 \ dM,‘ 4 1 
Ki M, fi j=i+l Qi aj-1/ L=— j=i Aj A; Ai+1 aj-i) M, a 


[40 } 


This equation shows the percentage change in the mass ratio 
of the 7th stage produced by a given percentage change in the 
propellant mass of that stage or any succeeding stage, or by 
a given percentage change in the inert mass of that stage or 
any succeeding stage including the payload. It should be 
noted that for 7 = N + 1 the factor ¢;(1 — A,)/A,; is unity. 

The absolute change in burnout velocity of the ith stage ob- 
tained by differentiating Equation [3] is 


dy = Vide; + — | [41] 


Consequently, in the most frequent circumstance where de; 
and d{git,“° | are both zero, we have 


dy = [42] 
and therefore, from Equations [30 and 40] 


1 1 \ aM, 
j=it+l aj-1/ M,% 


( 


The total change i in vehicle burnout veloc ay is found 7” sum- 
ming Equation [43] over all stages. When rearranged the 


result is 
N 
dM, dM dM, 
avy, Pi G) N M, [ 4] wy! 
where the coefficients are = 


= ati 


al — {fie 


( 

( 
( 


A — Qs O31 — S204 


A; Ay 1 C5 Ay 1 — C5 Gs; Os 1- C5 as, C5 


2 1— A; 


1 
— 


= ——_ - - - — 


A A 1 — Lies Ae al — te; Ay 1 — 65 mi fin 1 — ¢5 


| 


ARS JouRNAL 


ex] 


He 


39 | 
3 
eels Gn. n 

> 

i 

s 
h: 
Hi 
| 
in 
of 
at 
= 
mu 

534 JU 

= . 


1 & 1 & 
ty= (22... 2-8 (42.2 


It may be noted that the subscripts on p; and e; denote the 
contribution of the ith stage to dv, whereas the subscript on 
Ly denotes the total number of stages for the rocket being 
considered, 

In the special case where the stage ratio a;, the stage pro- 
pellant fraction ¢; and the effective exhaust velocity c; are the 
sume for each stage, it is apparent that the coefficient py, ex 
and Ly all converge to definite limits when VN ~ ©. By 
niaking use of the fact that the sum of the series 1 + x + 2? 
+ ...is1/(1 — x) when 2? < 1, we find these limits to be 


li is easy to show that p.. is equal toe... Hence for the upper 
stages ({> 3) of a multistage rocket, increasing the propellant 
mass or decreasing the inert mass by the same percentage will 
have an approximately equal effect on the burnout velocity. 

It is important to note that the coefficients p; can never be 
negative, since from Equation [18] 


Qi-1 1— ¢; a; — A,(aj-1) aA; — 1)/A; + 1 
[51] 


Hence p; is always positive. Physically this means that the 
increase in v, brought about by increasing the propellant mass 
of the ith stage can never be offset by the decrease in v, re- 
sulting from the associated increase in “payload” of the previ- 
ous stages. 

From the foregoing equations we can immediately deduce 
the effect of changing the propellant density in any stage. If 
the volume of propellant and its effective exhaust velocity 
are held constant in each stage so that the only change is in 


the propellant density p,;, then in [44] dM,@©/M,© = 
dp:/p:,dM© = 0, anddM, = 0; hence 
dy = pi — 
i=1 pi 


We may also deduce the tradeoff ratio between effective 


exhaust velocity and propellant density for any stage. From 
Equations [41 and 44] it is apparent that 
N 
day = Vide; + pi — [53 | 


i=] pi 


Hence if all stages remain unchanged but one, then dv, will be 
zero if 


‘hus a given percentage increase in effective exhaust velocity 
iust be greater than any accompanying decrease in propellant 
ensity multiplied by the factor (p;/c;0;) in order for there 
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[54 | 


to be a net increase in burnout velocity. Conversely, any 
attempt to increase vehicle performance by increasing the 
propellant density can be accompanied by a decrease in 
effective exhaust velocity no greater than that corresponding 
to Equation [54]. 

To illustrate the magnitudes of the coefficients p,;, e; and 
Ly, let us suppose that a multistage rocket is constructed with 
the motor propellant fraction A; = 0.875 for each stage, the 
stage propellant fraction ¢; = 0.70 and the stage ratio a; = 5. 
If the effective exhaust velocity is the same for each stage, 
then the coefficients are as follows 


= 0.700c, ps = 0.373c, ps = 0.3086, 

Py =0.295c,..., Po = 0.292c 
= = 0.280e, €; = 0.289c, 

= 0.29lc,...,e. = 0.292c 


L, = 0.560c, Ls = 0.579e, 
L, = 0.582e¢, . . 


L, = 0.467c, 
Le = 0.583¢ 


It is of interest to observe that in this example the effect of a 
given percentage increase in propellant mass is about twice as 
useful in the first stage as in any other stage. 

To extend the usefulness of Equation [44], recall that the 
maximum range of a ballistic missile over a nonrotating 
spherical Earth is 


R = — 0<R<ar [55] 


where ro = 2.09 X 107 ft = 3440 nautical miles is the radius 
of Earth. This equation may also be written 


v,2 sin (2/2ro) 


= [56 | 
1 + sin (R/2ro) 


and on differentiating gives 


gro [1 + sin (R/2ro)]? 2ro 


Division of Equation [57] by [56] gives | 


diy 
Up [1 + sin (R/2ro)] tan (R/2ro) 


or 


dR » Jtan [1 + sin (R/2ro))l dv» 
R/2ro 


[59] 


This equation may be used in conjunction with Equation [44] 
for estimating the effect of changes in rocket characteristics 
on missile range. 


Nomenclature 


c = effective exhaust velocity 

e€ = coefficient for inert parts, in Equation [44] 

g = effective component of gravitational acceleration — 
Ly = coefficient for payload, in Equation [44] 

M»y = mass of vehicle at propellant burnout 

M. = empty mass of rocket stage 

M, = mass of payload 


N 47 
[47] 
a 
x “Tey 
[57 | 
[58] } 
= 
} 
| 
|| 
| 


radius of Earth, 2.09 x 107 ft 
mass ratio, Equation [5] 


mass of propellant 
initial mass of vehicle stage 
number of stages in multistage rocket 
coefficient for propellant in Equation [44] - 
maximum range of missile over spherical Earth 
duration of burning - 
burnout velocity a 
dimensionless velocity increment, Equation [3] 
vehicle mass in relation to payload mass, Equation [22 


Qa 
ai stage ratio, Equation [15] 
oi stage propellant fraction, Equation [11] 


Orbits for Planetary Satellites From 
Doppler Data Alone 


pr 
» tele. &- 


HEN a satellite is put in orbit around some body other 

than Earth, it will then become possible to ascertain for 
that body the same kinds of data that terrestrial satellites 
have recently obtained. In this way we may expect soon to 
map the magnetic fields of the moon and planets, and the 
radiation belts associated with them. We may also obtain 
the oblateness terms in the gravitational potentials of these 
bodies, the properties of their exospheres and many other re- 
sults of high importance. 

However, in order to achieve results of this kind, it will 
usually be necessary to establish the satellite orbit with con- 
siderable precision. The vehicle will be thrown into a 
satellite orbit by some sort of end course guidance based, 
presumably, upon midcourse observations of its ballistic 
trajectory, but it is unreasonable to expect that the end- 
course guidance will be sufficiently exact to yield just the 
desired satellite orbit. A method is therefore needed for 
establishing the orbital elements after the application of end- 
course guidance. 

This requirement will not be easily satisfied by conven- 


Received Aug. 4, 1959. 

1 Consultant, Missiles and Space Systems Engineering; also 
Staff Member, Mt. Wilson and Palomar Observatories, California 
Institute of Technology, Pasadena, Calif. 


536 


motor propellant fraction, Equation [14], = 1 — eof ref. 
1) 
= propellant density ie 
Subscript 1 and superscript to 7th stage 


Ai 
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When an artificial satellite is put in orbit around another planet, it will be difficult or impossible 
from terrestrial observations to obtain its direction with sufficient precision to define its orbit. 
However, if the satellite carries a suitable transmitter, Doppler observations are capable of giving 
precise values for the line of sight component of the velocity vector. 

possible to obtain all the orbital elements and the mass of the planet. 
used by astronomers for the discussion of spectroscopic binary stars is here adapted to the interpre- 
tation of Doppler observations of a planetary or lunar satellite. 


determination for a satellite of Venus. 


From such data alone, it is 
The analytical apparatus 


An example is given of the orbit 


tional methods. At planetary distances, vehicles of the ex- 
pected order of size will be far too faint for optical detection 
by reflected sunlight. At the distance of the moon, the 
stellar magnitude of a perfectly reflecting sphere of 1-m 
diameter will be 16. But as a satellite of Venus, the same 
sphere will be of the 29th magnitude when Venus is near in- 
ferior conjunction, and a magnitude fainter at superior con- 
junction. As a satellite of Jupiter, the sphere will have a 
magnitude of about 37. The extreme limit of the 200-in. 
telescope is the 24th magnitude, and a practical working 
limit is the 22nd magnitude. 

Although optical observations are therefore out of the 
question in most cases, the vehicle may nevertheless be de- 
tected by suitable radio telescopes, provided it carries a 
sufficiently powerful transmitter. Actually, the power re- 
quirements aboard the satellite will probably not be exces- 
sive. Rechtin? has recently estimated the capability of a 
microlock-type system in 1960 as a 30-cps information band- 
width at the distance of Mars, provided the vehicle carries a 
10-w transmitter and an antenna of 16-db gain. 

The detection of continuous wave signals from a planetary 


2 E. Rechtin, colloquium at California Institute of Technology, 
Pasadena, Calif., March 20, 1959. 
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satellite will not necessarily solve the orbit determination 
problem, however. Let us suppose that a vehicle is put 
in orbit around Venus, with a semimajor axis of order 
10 per cent greater than the planetary radius. As seen from 
Earth, the maximum angular subtense of the semimajor axis 
will then be 70”; at superior conjunction it will be only 
11”. In the case of a Martian satellite, these angles be- 
come 27” at opposition and 4” at conjunction. Now, 
these numbers represent the whole angular excursion of the 
satellite relative to the planet; to establish an orbit of any 
utility, we should have to require that angular measurements 
discern the satellite position with an accuracy at least some 
10 times higher. But this requirement is far beyond the 
capabilities of any existing or projected high gain antenna; 
the Jodrell Bank paraboloid, for example, has a theoretical 
resolving power of 10’ at 1420 me per sec. 

Despite the evident inadequacy of angular measurements 
for orbit determinations of this kind, there remains the possi- 
bility that, if the satellite transmits a signal of high frequency- 
stability, an orbit can be accurately determined from Doppler 
dita alone. The orbital velocities V in the examples cited 
previously are 23,000 fps for the Venerian satellite, and 11,000 
fps for the Martian satellite. At speeds of this order, the 
velocity error 6V that corresponds to a frequency error 6A is 
given by the relation 6V/V ~ 10° (6\/A). For an error 

6\/X| —& 10-8, which appears to be a conservative figure 
for the frequency-stability that can be achieved, the corre- 
sponding error |6V | amounts only to about 10 fps. In the 
special case of a satellite orbit that is circular around a planet 
oi unknown mass, the measurement of V leads directly to a 
determination of a. Moreover, the fractional errors of a 
and V will be related as Aa/a = 2|Av/v|. For a Martian 
satellite with V ~ 11,000 fps, a frequency stability of 1078 
therefore would yield a with an error of 0.2 per cent. This is 
evidently quite sufficient precision for many important in- 
vestigations. 

An analogous problem of orbit determination is often en- 
countered in conventional astronomical research. From the 
Doppler shifts of the absorption lines in stellar spectra, as- 
tronomers have determined for many stars the component of 
stellar velocity along the line of sight to Earth. When this 
velocity component is referred to an origin at rest in the sun, 
it is said to be the radial velocity of the star. In many stars, 
the radial velocity varies periodically in a way which indi- 
cates that each such star is in Kepler motion around a close 
companion. The companion is usually so close that the pair 
cannot be resolved in the sky by any existing telescope, and 
it is usually so faint that its light cannot be detected in the 
observed spectrum. Techniques have been developed for the 
derivation of certain orbital characteristics of such a spectro- 
scopic binary star from its observed velocity curve. In 
particular, a so-called “mass function” can be obtained, 
which in some cases leads to a determination of the sum of 
masses of the two stars comprising the pair. 

The techniques appropriate for the discussion of spectro- 
scopic binaries will now be adapted for the case of a planetary 
satellite. We shall suppose that from pulse Doppler observa- 
tions of such a satellite, it is possible to obtain a velocity 
curve, which gives the radial velocity of the satellite as a 
function of phase in its orbit. In the first instance, we shall 
suppose that the satellite orbit is a Kepler ellipse fixed in 
space relative to an inertial frame, and that Earth and planet 
(say, Venus) are also stationary in this frame. The Doppler 
effect is supposed to be measured by a stationary observer at 
the center of Earth, and the orbit is supposed to subtend an 
infinitely small angle as viewed from this position. In later 
sections, account will be taken of the various effects due to 
departures from these ideal conditions. 


Doppler Orbits Observed in One Orientation 


Consider the unit sphere with center O at the center of 
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Fig. 1 Unit sphere with center at Venus. The pole of a satellite 
orbit lies in the direction of S, periplanet in the direction of P 


and Earth in the direction of E 


Venus. In Fig. 1, the termini are indicated for the following 
unit vectors: 


E, in the direction from Venus to Earth 

S, in the direction of the pole of the satellite orbit 

N, in the direction of the ascending node of the satellite 
orbit against the plane perpendicular to E (‘‘plane of the 
sky’’) 

X, in the direction of north along the hour circle through 
Venus 

Z, in the direction of —E 

P, in the direction of periplanet 


To define the orbit plane, it will suffice to give the inclination 
zi, which is the angle between E and S; and the position angle 
of the ascending node Q, which is the angle between X and 
N. The orientation of the orbit in its plane is then deter- 
mined by the longitude of periplanet w, which is the angle 
between N and P. The ranges of these three angular ele- 
mentsare: 0<1<7;0 <Q < 27;0 < w < 

Let r be the instantaneous radius vector of the satellite 
relative to Venus; then the true anomaly v is the angle be- 
tween P and r, measured in the sense of the orbital motion 
around S. The component of r along Z is therefore seen in 


Fig. 1 to be _ 


and, by spherical trigonometry, this is easily shown to be 
z=rsinisin (v + w) [2] 


As viewed from Earth, the radial component of the velocity 
of the satellite is then 


2 = frsin7zsin (v + w) + rd sin7 cos (v + w) [3] 
where the dot denotes differentiation with respect to time. 
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But from the theory of Kepler motion it may be shown (1)* 
that 


2m ae sin v 4 2m a(1 + e cos v) 
r= ro = [4] 


where 
a = semimajor axis of the orbital ellipse wie ¢ 
P = orbital period 


Introducing these expressions into Equation [3] and simpli- 
fying, we find that 


z= Kle cos w + cos (v + w)] [5] 
where the amplitude of velocity variation is i 
K = 2rasin i/PV1 [6] 


If the distance between the observer and Venus is increasing 
at the constant rate y, then this parameter will appear as an 
additive constant in the right member of Equation [5]. 

Let the time ¢ be measured from an arbitrary origin, and 
let the phase be defined in the relation 


® = (27/P)t 

The observed velocity curve then gives z as a function of ®. 
From this curve, it is possible to deduce the values of K, 
e, wand y. The techniques for obtaining these parameters 


are described, e.g., by Smart (1) and Aitken (2). The 


velocity curve also yields the epoch T of periplanet, which is 
an instant of minimum separation of planet and satellite. 
Knowing P, e and K, it is also possible to find from Equation 
[6] the quantity a sin 7. 

To obtain the mass function f(m) we now adduce the har- 
monic law 


* Numbers in parentheses indicate References at end of paper. 
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Fig. 2 Orthographic projection of the unit sphere at Venus. If 
the planetary mass is assumed known, then any one velocity 
curve for a satellite restricts S to two small circles s and s’, and 
_ P toa third small circle p 


In this relation, m is strictly the sum of masses of planet and 
satellite. For our application, it is essentially the planct 
mass. Again, a is strictly the semimajor axis of the satellite 
orbit around the planet as focus. For our application, it is 
essentially the same as the observed quantity, which is the 
semimajor axis of the satellite orbit around a focus at thie 
center of mass of planet and satellite. The mass function may 
now be defined as 


f(m) = m sin? i 
and, by Equation [8] it is 
4 
(a sin 
fim) = [10 


in which the right member contains only known quantities. 
When the mass may be considered known, the quantity 
sin 7 can be found directly from Equation [8] in the form 
(a sin 7)8 


[13] 


and a can then be found as well. Fig. 2 will illustrate tho 
geometry of this case. The circle represents an orthographic 
projection of the unit sphere centered on Venus. Equation 
{11] restricts the terminus of S to the two small circles 
designated s and s’ in Fig. 2. Similarly, the known value of 
w then restricts the terminus of P to the small circle desig- 
nated p. A particular case satisfying these restrictions is 
indicated by S* and P*, the corresponding orbit plane being 
designated f. 

This completes the discussion of the information that can 
be obtained from a velocity curve of a satellite under the as- 
sumption that its orbit is a fixed Kepler ellipse seen always in 
the same orientation. We shall next discuss the additiona! 
information that can be obtained when we take into account 
that the satellite orbit may be observed under different orien- 
tations as a consequence of the orbital motions of Earth and 
planet around the sun. 


Doppler Orbits Observed in Diverse Orientations 


_ Observations in Three Non Co-Planar Directions 


We shall now show that if velocity curves can be obtained 

for a fixed Kepler ellipse as viewed from three non co-planar 
directions, then all the orbital elements of the ellipse are 
unambiguously determinate, as is the mass of the planet. 
_In most of the foreseeable applications of Doppler techniques 
_ to satellite orbit determinations, the observations will, in fact, 
be very nearly confined to the plane of the ecliptic. More- 
over, the principal purpose of Doppler observations is likel) 
to be the determination of the shape, size and orientation of 


the satellite orbit, and for this purpose it will usually suffice 
_ to assume that the mass is known from independent astro- 
ate determinations. 


Nevertheless, it will now be appro- 
priate to consider the general non co-planar case, with un- 
known mass, before proceeding to the special cases which arc 
more likely to be of early practical importance. 

We shall use the same notation as in the previous section. 


We may take S and P to be constant vectors when referred to 


an inertial frame. However, due to the motion of Venus and 


Earth in their orbits around the sun, the vectors E, N and 


X will vary with time. Their values at the epoch ¢; will be 
designated E;, N; and X;. If E;, E., E3; are any three non 
co-planar unit vectors, we may write 


But from Fig. 1, we see that 
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Hence 


N, xX N, sin sin Is = x E,)-E.|S [14] 


By Equation [12] ae 


Sx E, = AoE» x E, + Q3E; xX E, 


Therefore 
N, x Ne sin sin 72 = a3[(E; x E,)-E2|S - {16} 
Moreover, Fig. 1 shows that 
N, x Nz = S sin — we) [17] 


Hence Equation [16] yields 


= sin (@ — Ws) sin 7, sin x E2) [18] 


and similar expressions can be obtained for a; and a, by 
cyclic permutation of the subscripts. 
Since S is a unit vector, we find by squaring Equation [12] 


that 


2D aa; cos = 1 [19] 


j=] ixj 


where I’,; is the (known) angle between E; and E; 
cos | E;-E; [20 | 


Hence Equation [19] may be written as - 


3 
at = (a?a;)? +2 (a? a;)(a? a;) cos Ti; 
j=l 
and in this equation all quantities appearing in the right mem- 
ber may be taken as known. Having found the semimajor 
axis in this equation, we may now proceed to obtain sin 2; 
from the known values of a sin 7;.. Equation [18] then yields 
the coefficients a; and therefore determines the vector S. The 
a; also determine the quadrant of each 7;, for Equation [12], 
and Fig. 1 gives us 


cos 1 = E,-S = + COS + Qs; cos [22] 


etc., for cos 72 and cos 13. 
With S known, we may now compute N; in Equation [13]. 
Moreover, we see in Fig. 1 that 


P = N; cos w; + Sx N; sin w; wb. [23 ] 
and this enables us to find the vector P. We also have 
E; sin Q; = X; x N; cos Q; = X;-N; [24] 


and since X; is known, these relations determine Q;. Finally, 
we may obtain the mass from the harmonic law, Equation 
[8]. In the practical case where m may be considered as 
independently known, Equations [8 and 21] provide a con- 
straint upon the observed values of a sin 7; and w;._ We shall 
see in the next section that when m is known, there is an 
additional constraint upon the parameters a sin 7; and o; 
that are associated with any pair of directions E;. 


Observations From Co-Planar Directions 


In the practically important case where the satellite orbit 
can only be observed from directions lying in a plane, the 
argument of the last section fails because S cannot in general 
be analyzed into components along any three E;, as in Equation 
[12]. To discuss this case, we define a unit vector Ey which 
is perpendicular to the plane of E, and E> 


E, sin = E, xX E, [25] 


The orbit is then never observed from the direction Eo, but 
only from £;, EZ, and (possibly) other directions in their plane. 
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sin (W; — We) sin 2; sin 72 
= 
E.) 


? sin (a — sin 2; Sin (26 
sin 
We also have 
2 « 
S= 27] 
j=0 


and therefore 


= a, + ae cos = cos 
E.-S = a, cos Ty. + a2 = cos t2 
It follows that 
a, sin? = cos 7; — cos cos 
a sin? = cos — cos cos 
But the square of Equation [27] is 
Qo? + + + cos = [30] 


and if we substitute the values of the a; into this equation, we 

obtain after some simplification 

[((2 — sin? — (Ay? + Ae?)a? + sin? — — 
4a‘ — A,2)(a? — As?) = 0 [31] 

where 


A; = asin? 
Ay = asin 


The observed values of Az and — therefore con- 
strain the semimajor axis a to satisfy an equation of the eighth 
degree, namely Equation [31]. In case we may regard m as 
independently known, then a follows directly from Equation 
[8], and Equation [31] becomes simply a constraint upon 
sin 2), sin 7% and w; — w:. In either case, however, there will 
be two values of 7; compatible with the observed .1,; and two 
values of 7 compatible with the observed 2. We may 
write these as 

or 

Io = I, or Io = I,’ => 

Suppose that when we compute (a, a2) in Equations [29], 
we use the pair (J), /2), and that this produces an (a1, a2) which 
satisfies Equation [30]. Then it is clear that will also 
produce an acceptable (a), az). In general, however, the pairs 
and (Jo, will then not correspond to possible solu- 
tions. It follows that the orbit plane may have one of two 
possible orientations, as shown in Fig. 3. 

To find the effect of this ambiguity on the vector P, we 
write it in the form 

2 


E,-P = B, + cos Ty. = —sin 7, sin 
E,-P = B, cos + Bs = —sin sin 


| 
|= As in Equation [18], we may then wri 
uquation , we may then write 
t 
| 
| 
| 
| 
Referring to Fig. 1, we see that 
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Corrections for Diurnal Motions, etc. _ 


It is now necessary to discuss the various effects which may 
be produced as the result of departures from the ideal condi- 
tions specified earlier. In the first place, we note that, 
during one cycle of the satellite, there may be an apprecialle 
change in that component of the observer’s diurnal motion D 
which lies along the line of sight to the satellite. But since 
D is known, and also the direction —E to the satellite (or, 
at least, to its primary planet), it is a simple matter to remove 
from all the observed radial velocities the contribution — D E 
of the diurnal motion. 

In the idealized theory, where the observer was supposed 
to be at the center of an infinitely distant Earth, it was not 
necessary to differentiate among the following vector d s- 
tances in Fig. 4: 


pE, from the center O of the primary body to the center € of 
Earth 

\F, from the satellite S to the center of Earth 

lG, from the satellite to the observer @ on the surface of 
Earth 


When the observations are sufficiently precise, however, it 
will be necessary to take into account that they yield tie 
quantity /. Certain small corrections must then be applicd 


Fig. 3. Degeneracy of the solution when all observations are to | to obtain the quantity z, for which the foregoing theory 
confined to one plane, e.g., the ecliptic. If one of the illustrated has been developed. In the case of a lunar satellite, for ex- 
orbits satisfies the velocity curves, the other orbit will do so ample, the orbit may subtend an angle of more than a degree 
equally well at Earth, and the radius vector R to the observer may sub- 
tend an angle about as large at the satellite. ; 
To find the corrections that must be applied to / in order 
to obtain 2, we note in Fig. 4 that 


r= p? 2pz [39] 


= \2+ R2+ 2\F-R [40] 


_ where the notation is the same as in the previous sections 


B sin? Py, = cos sin sin — sin sin and Fig.4. Therefore 
Be sin? Ty, = cos Ty: sin sin — sin sin we 
The ambiguities of Equations [33] therefore leave 8; and G2 -* p p p 
unambiguously determined. Moreover, it is clear that and 


and by Equations [27 and 34] this is ; a 1+ (1/A)F-R 


Fig. 4 Corrections of the velocity curve of a lunar satellite for diurnal motion of a terrestrial observer, parallax 
and finite angular subtense of the satellite orbit 
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Bo + + A282 + + a26:)cosT2 = 0 [38] But-to a close approximation, we shall always have 
which determines {8 and therefore P, when S is known. S 2 L~rA—p 2/\ <1 r/\ <1 
differs from S’ only in the signs of a and az (see Fig. 3), and Ri R/\<X1 E~F~>G 
therefore reverses its sign with a; and a». 
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TIME (MINUTES) 
Oo 60 120 180 240 


30 i i i i i i i i i 
_ DECEMBER 1, 1959 
10-4 
2 
—20- 
ve 
—30 T T 
0 90 180 270 360 90 180 
MEAN ANOMALY (DEGREES) 
e = 0.30 a SIN i = 4,772 NAUTICAL MILES 
T = DEC 1.683 U.T. @ = O DEGREES 


Fig. 6 Velocity curve for a satellite of Venus, Dec. 1, 1959. 
The dashed curve is for the image orbit 


TIME (MINUTES) 


Fig. 5 Schematic representation of a Venus satellite orbit. 30 


Velocity curves are supposed to be observed on Dec. 1, 1959 

ind June 1 and Dec. 1, 1960. The positions of Venus and Earth JUNE 1, 1960 
on these dates are indicated, together with the transfer orbit 
corresponding to a launch date at Earth of June 7, 1959 


nN 
1 


1 


RADIAL VELOCITY x 10-3 (FT/SEC) 
a 


Cherefore a first approximation to Equations [41 and 42] 


yields 
=h-p 1 + (d/dt\(E-R) [44] 
‘The second of these equations gives \ in terms of the measured 0 90 180 270 360 90 180 
radial velocity and the (known) correction for the component MEAN ANOMALY (DEGREES) 
of diurnal motion along the line of sight to the primary. e = 0.30 a SIN i = 4,552 NAUTICAL MILES 


T = JUNE 1.592 UT. @ = 134 DEGREES 


The first equation then gives a first approximation to 2, 
Fig. 7 Velocity curve for a satellite of Venus, June 1, 1960 


in terms of \, and the (known) radial velocity of the pri- 

mary. From 2, we may then find a first approximation to i rh i 

the orbital elements, by the methods detailed previously, and i) ; 

the corresponding approximation to r as a function of phase. a a 
Let this approximation be designated 7. From it, we : 


TIME (MINUTES) 


may then compute an improved approximation of each of the a. 
= pE—n +R (45) DECEMBER 1, 1960 
— 4“ 204 
, (lo/A2) — (d/dt)(F2-R) 46] 
= 2 
8 
If 2. differs significantly from 2 it is then possible to reiterate 8 -10+ 
this procedure until the observations are adequately repre- 2 7 
sented. 8-204 af 
There still remain to be discussed all those effects which = 7 
arise because of perturbations, i.e., because the orbit is not 
really a fixed Kepler ellipse. In order to view the orbit of a ) 90 180 270 360 90 180 
planetary satellite in appreciably different orientations, it MEAN ANOMALY (DEGREES) 
will usually be necessary to wait for intervals of several e = 0.30 a SIN i = 4918 NAUTICAL MILES 
months. However, during such intervals the satellite may T = DEC 1.501 UT. @ = 267 DEGREES 
have made a thousand cycles, and in the most interesting and Fig. 8 Velocity curve for a satellite of Venus, Dec. 1, 1960 
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ate 

24 
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informative cases the orbital elements will be appreciably per- 
turbed during these intervals. The techniques required for 
the analysis of perturbed satellite orbits are a natural exten- 
sion of the ones developed here. They will be given in detail 
in a later paper. 


Example: Doppler Orbit for a Satellite of Venus 


To illustrate the kinds of velocity curves that will be ob- 
tained for planetary satellites, let us consider the following 
hypothetical example. On June 7, 1959, a vehicle is launched 
in a Hohmann transfer orbit to Venus, as illustrated in the 
schematic diagram of Fig. 5. The vehicle then arrives at 
Venus’ orbit on Nov. 1, 1959; at this time, the last-stage 
motor provides an impulse sufficient to throw the probe into 
a satellite orbit around Venus. 

The probe carries a suitable transmitter to permit Doppler 
measurements at a terrestrial receiver. These observations 
are first made on Dec. 1, 1959 (Fig. 5). After correction 
for the known radial velocity of Venus and for the observer’s 
diurnal motion, the observations yield the (solid line) velocity 
curve of Fig. 6, with a period of 173.4 min. The other spec- 
troscopic elements that can be deduced from this velocity 
curve are also indicated in Fig. 6. If we adopt the astronom- 
ically determined value for the mass of Venus, we can also 
obtain a and sin 7. separately. In any ease, the quadrant of 7 
is indeterminate from this velocity curve, and so is the direc- 
tion of the line of nodes against the plane of the sky. 

Velocity curves are determined again on June 1 and Dec. 1, 
1960, and the corresponding elements are found in each case, 
as indicated in Figs. 7 and 8. When these three sets of spec- 
troscopic elements are analyzed with the apparatus pre- 
viously described, it is found that the pole of the satellite orbit 
must lie in the direction of the unit vector with components 


S, = 0.300 S, = 0.400 S. = 0.866 


The x-axis here points to the vernal equinox, as usual, and the 
z-axis to the north pole of the ecliptic. In addition, we find 
that the unit vector towards periplanet has the components 


P, = —0.377 P, = 0.884 P, = —0.278 


The values of inclination 7 and position angle of the node Q, 


respect to the plane of the sky in the following table 
1, deg at Q, deg 


Dec. 1, 1959 67 
June 1, 1960 119 253 
Dec. 1, 1960 71 ‘ 288 


If these orbital elements are now referred to the plane of th« 
ecliptic, as is customarily done for the orbital elements ot 
planets, they are found to be 


a = 5200 nautical miles = 1.55 Venus radii 


e = 0.300 

i = 30deg 
Q = 143 deg 
w = 326 deg 


P = 1734 min 
T = Dee. 1959, 1.683 UT 


The inclination of the orbit of Venus relative to the plane ot 
the ecliptic is only 3.4 deg, and the geocentric latitudes of the 
planet on the three dates of observations are, respectively 
2.2, 5.1 and 2.2 deg. The configuration is, therefore, not far 
from the degenerate case treated earlier, where the velocity 
curves are equally well satisfied by this orbit and by its mirror 
image in the plane of the ecliptic. To indicate, in this case, 
the discrimination between the true orbit and its image, the 
velocity curve corresponding to the image orbit has also 
been computed for Dec. 1, 1959, and is shown as a dashed 
curve in Fig. 6. 
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Detonation Wave Hypersonic Ramjet 


tonation wave are presented. 


ty 


YPERSONIC ramjets employing a standing detonation 
wave were first mentioned by Roy in 1946 (1).4 Al- 
though detonations had never been stabilized in any device 
at that time, their properties could be computed and potential 
performance assessed. Dunlap (2), Nicholls (3), Roy (4), 
Partel (5), Chinitz (6) and Weber (7) have recently presented 
papers which contain discussions of such a device. 

Several years ago a number of laboratories in the United 
States began to explore techniques for releasing chemical 
energy steadily in supersonic flow. Research groups at the 
Fairchild Engine Division and the University of Michigan 
have succeeded in producing steady detonations in the 
laboratory. A summary of this work on stationary detona- 
tions has recently been prepared by Gross and Nicholls (8). 
Other pertinent work is referred to in (8, 9 and 10). 

In Fig. 1 is shown a typical Schlieren photograph of a 
standing strong detonation wave of the type found experi- 
mentally by Gross and Chinitz (10). The detonation wave 
shown is a hydrogen-air reaction, stable at Mach 3.15 in a 
supersonic combustion research tunnel. The oblique shock 
Waves originate from side wall wedges and, under zero fuel flow 
conditions, intersect in a small normal Mach reflected shock. 
Upon the introduction of fuel upstream, and with the proper 
stagnation conditions for ignition [see (10) ], a strong detona- 
tion wave forms. This wave is stable, steady and repro- 
ducible. It has relatively simple properties, considering that 
there are large gradients in the inflowing fuel-air ratio and 
consequently large temperature and pressure gradients on 
the downstream side of the wave. The local conditions im- 
mediately across the wave are predictable using classical 
methods (11). The regions of stable operation are discussed 
in (10 and 12) and consist of the “strong’’ branch (M2 < 1) 
solutions. An ignition hysteresis effect was found (10) 
which, although not completely understood, may permit 
detonation wave stabilization at low inlet stagnation tempera- 
tures. 
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air specific impulse, specific fuel consumption and thermal efficiency. 
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Performance of a hypersonic ramjet having a detonation wave combustion process is described. 
Flight speeds from Mach 2.5 to 10 are examined. 
Engine cycle performance characteristics are computed including 


Thermodynamic property changes across a de- 


The effect of component 


efficiency on cycle performance is examined and the optimum place in the cycle for the detonation 
wave is found. The subsonic burning ramjet is shown always to have superior thermodynamic 

performance compared to a detonation wave engine, but the differences are often minor. Other 
; differences and some design features are discussed. 


With stationary stable detonations in the research labora- 
tory and an increasing awareness and understanding of their 
actual properties, assessment of the performance of a pro- 
pulsion device employing such a combustion system was de- 
sirable. This formed the basis and motivation for this*de- 
tailed study of a detonation wave ramjet cycle. 


Engine Cycle 


The engine cycle which is studied here is a typical one of 
compression, combustion and expansion. The basic differ- 
ence is the condition under which combustion takes place 
namely, a standing detonation in supersonic flow. 

The parameters which we determine for this cycle are: 
Air specific impulse, specific fuel consumption, thermal 
efficiency and exit to inlet area ratio. 


Schlieren photograph, strong detonation wave. Mach 


3.15, hydrogen-air, from Gross and Chinitz (10) 
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\\— STATION NO. 3 
DETONATION 


Altitude 130,000 feet, ICAO atmosphere ae 
Hydrocarbon fuel 

Chapman - Jouguet conditions at nant 
Mach number ahead of detonation - 4.0 

Non - dissociative case 


Fig. 2 Schematic of detonation wave ramjet engine, flight Mach 
number 6.0. Altitude 130,000 ft, hydrocarbon fuel, C-J detona- 
tion at Mach 4 
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Fig. 4 Equivalence for C-J points vs. pre-wave static tempera- 
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These parameters have been determined as a function of: 
Flight Mach number (2.5 to 10), detonation Mach number 
(2.5 to 4.0), inlet efficiency and exhaust nozzle efficiency. 

The calculations were all performed for a typical hydro- 
carbon fuel whose heating value AH is 8.78 X 105 Btu/mole 
and for flight conditions corresponding to 130,000 ft, as 
defined by the ICAO atmosphere. 

A typical engine configuration is shown in Fig. 2. The 
inlet consists of a passage in which supersonic compression of 
air occurs adiabatically from the flight Mach number to the 
combustion or detonation Mach number. The kinetic effi- 
ciency of compression 7x is defined as the ratio of the kinetic 
energy in the diffused stream (when that stream is isentropi- 
cally expanded to ambient free stream static pressure) to the 
free stream kinetic energy. Values of nx ranging from 1.) 
(ideal case) to 0.94 were employed and Yo was treated as « 
constant equal to 1.4. Fig. 3 shows the total pressure re- 
covery as a function of flight Mach number for several valucs 
of inlet kinetic efficiency. 

Inlet efficiency as used in this report applies only to diffu- 
sion from flight Mach number to the pre-wave detonation 
Mach number. This usage is at variance with usual practice, 
in which the inlet efficiency refers to diffusion from free strea!n 
to a low subsonic Mach number. In our context, an inlet 
efficiency of 0.94 is a pessimistic assessment of performance 
and represents a lower boundary for design expectations. 
For the flight speeds from Mach 5 to 7, inlet kinetic efficiencies 
of approximately 0.96 to 0.98 are realistic. The static ten- 
perature at the end of diffusion and just upstream of the 
detonation wave is 


4 


1+ — 1)/21Met 
> = 


Pi = Pull + [(Yo — [3] 


where 


Py | (7 
=—<]— 1 — {| — 
11 P, NK P,, 


_ Fuel injection is assumed to take place upstream of the in- 
let, so that a uniform fuel-air distribution exists at the detona- 
tion wave. Losses accruing from fuel injection have not been 
included in these computations. For an assessment of these 
losses the reader is referred to (2). 

Conditions across the normal detonation wave are treated 
by the solution of the classical conservation equations. The 
methods and computer program were adapted from (11). 
Most of the calculations presented here employ the non- 
dissociative case, with thermodynamic properties a function 
of temperature. The reactants and products of combustion 
result from the chemical change 


apCrHy + + + aA 
BiCOz + B2H20 + B;02 + BiN2 + BsA [4] 


and are treated as ideal gas mixtures. Combustion efficiency 
is assumed always to be 100 per cent. In general, it is ad- 
vantageous to have detonation occur under Chapman-Jouguct 
conditions (7), i.e, M. = 1. For those conditions which 
cannot be satisfied by this requirement, a strong detonation 
solution is employed. The conditions across the wave that 
are of interest in this analysis are shown in Figs. 4 through 8 
These data represent the nondissociative solution of a 
hydrocarbon fuel reacting with air in one-dimensional steady 
flow. The nondissociative calculation is presented here in 
detail because of its relative simplicity and because in a com- 
parison study it represents only a small error. In (11) dis- 
sociation as it affects the wave is discussed. The dissociative 
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equilibrium calculation may also be at variance with reality 
in that reaction rates may not permit the attainment of 
equilibrium. A dissociative equilibrium cycle calculation at 
Mach 6 gives, for example, an increase in specific fuel con- 
sumption of about 7 per cent. 

In a rapidly moving gas mixture where chemical equilibrium 
does not exist and where the constituents are either dissoci- 
ating or recombining, the extent of dissociation or recombi- 
nation in an arbitrary region of the mixture will depend 
upon the local mean temperature of the mixture, the partial 
pressure of the constituents and the mean residence time of 
the molecular species. Since the residence time in the com- 
bustion zone in the conventional engine is far longer than in 
the wave engine, it might be expected that in whichever 
way the equilibrium is shifting, the shift will be more com- 
plete for the long residence time flow. 

Figs. 6 and 7 show the variation of the wave total and static 
pressure ratio with respect to pre-wave static temperature. 
‘The rapid increase in the static pressure ratio should be noted 
in the regime of strong detonation solutions. The rise occurs 
because of the very rapid decay of the wave exit Mach num- 
ber, as shown in Fig. 8. The total pressure ratio in the 
strong detonation solution region is not heavily affected by 
pre-wave static temperature except at the extreme Mach 
number conditions. It is interesting that the total pressure 
ratio actually moderates at the lowest pre-wave Mach num- 
ber in the strong detonation solution region. In general, 
however, the loss increases in magnitude and gradient with 
increasing pre-wave Mach number. 

The exhaust nozzle expands the combustion products adi- 
abatically with a discharge coefficient C'p of either 1.0 or 0.98. 
The nozzle exit temperature is given by 

T; = T{Po/P2)~ 9/7 [5] 
where 
= [¥(T2) + ¥(Ts)]/2 


and is computed for a gas composed of the products of com- 
bustion as specified by Equation [4]. The exit velocity is 
computed from 


where 
hu = CnT2}1 + [(v2 


The cycle performance can now be evaluated. The 
parameters of interest with their associated equations are 
given in the following: 

Air specific impulse 


Fy V3 =. Vo 
I 
Specific fuel consumption 
SFC = (f/a)s3600/Yl4 
Thermal efficiency 
(f/a)s( AH/m;) as 
to inlet area ratio 
[10] 


‘quation [7] applies for an expansion process to ambient 
pressure. This may be undesirable because of the large 
exit area and hence drag. The desired degree of expansion 
is a detailed question concerning an engineering balance be- 
tween thrust, weight and external drag. We assume a fully 
expanded nozzle for comparison purposes. 
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Fig. 8 Wave exit Mach a vs. pre-wave static tempera- 


Results 


Figs. 9 and 10 show some basic data for evaluating engine 
performance. In these plots, air specific impulse, specific 
fuel consumption and thermal efficiency are examined as 
functions of flight Mach number, pre-wave Mach number, 
and inlet and exhaust efficiency. In each figure where the 
performance characteristics for a constant pre-wave Mach 
number are defined, a base reference curve has been included. 
This curve represents the performance of an ideal conven- 
tional engine operating on the same fuel-air schedule as the 
associated wave engine. The ideal engine is described in a 
later section. At a flight Mach number of 6, the detona- 
tion wave engine makes its nearest approach to the idea 
engine. 

In the range of pre-wave Mach numbers from 2.5 to 3.0 
the difference between ideal wave and conventional engine 
performance becomes great. For example, at M, = 2.5 and 
M, = 6.0 the wave engine delivers 96 per cent of the impulse 
of the conventional engine, whereas at M,; = 3.0 the impulse 
ratio falls to 87 per cent. From M, = 3.0 to M, = 3.5, 
approximately the same ratio between conventional and wave 
engine impulse is maintained. Beyond M, 3.5 the ratio 
again falls rapidly to a value of 71 per cent at M, = 4.0. 

Increasing the fuel-air ratio increases the specific impulse, 
whereas wave loss decreases the impulse. At the lowest 
pre-wave Mach numbers the slope of the fuel-air ratio with 
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Fig. 9 Performance parameters, detonation wave ramjet, 
M, = 3.0 


pre-wave Mach number is not as steep as the slope of the wave 
losses. This causes a rapid decline in specific impulse. 
In the middie region of pre-wave Mach numbers both rates 
of change steepen, but the fuel-air ratio is more drastically 
affected than the wave loss. This results in a mild impulse 
plateau which then drops off as the wave loss again assumes 
control. 

The inflection of the impulse ratio is reflected in the shape 
of the thermal efficiency curves shown on Fig. 11. The 
variation of air specific impulse and thermal efficiency can be 
examined with respect to pre-wave Mach number. Optimum 
impulse at flight Mach number 4.0 is achieved with a pre- 
wave Mach number of 3.0. At flight Mach number 6.0, 
optimum impulse occurs at M, = 3.75, and at flight Mach 
number 8.0 the optimum point is again M, = 3.0. The shift 
in optimum pre-wave Mach number results from the inter- 
play of pre-wave temperature and wave losses. Ata constant 
pre-wave Mach number the permissible (Chapman-Jouguet ) 
fuel-air ratio increases with increasing pre-wave tempera- 
ture. At low flight speeds, if My = 4 can be considered low, 
minimizing wave losses in the low fuel-air ratio range is essen- 
tial. Ata flight Mach number of 6.0 the optimum pre-wave 
Mach number is increased to take advantage of a more favor- 
able fuel-air ratio. At a flight Mach number of 8.0 the pre- 
wave temperature is sufficiently high to permit a large fuel- 
air ratio at a reduced pre-wave Mach number, thus increas- 
ing the optimum specific impulse further. 
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Fig. 10 Performance parameters, detonation wave ramjet, 
M, = 4.0 


Fig. 12 contains a plot of exit to inlet area ratio as a func- 
tion of flight Mach number for a fully expanded nozzle. 
From a mechanical design standpoint this figure summarizes 
a very significant parameter. The conventional perfect 
ramjet operating at constant fuel-air ratio has an area ratio 
curve with a mildly negative slope. As can be seen from Fig. 
12, the conventional engine operating on the wave engine 
fucl-air schedule has a gentle positive slope. The losses of 
the wave engine in conjunction with its inherent fuel schedule 
cause the magnitude and variation of the exit to inlet area 
ratio to be quite severe and from the designer’s viewpoint 
undesirable. However, burning at higher pre-wave Mach 
numbers during acceleration will have the effect of easing the 
exit to inlet area ratio requirements. Since the losses are 
greater with a high pre-wave Mach number, this will require 
a larger exit to inlet area ratio, and it may be possible to hold 
the actual design area ratio to more reasonable limits. Par- 
tial expansion nozzles would probably be used. 

Pre-ignition of the fuel-air mixture is a matter for concern. 
lig. 13 shows a plot of pre-wave temperature vs. flight Mach 
number. Included is a “possible pre-ignition area,” extend- 
ing from 1600 to 2500 R. The reason for the uncertainty 
lies in the fact that time is also involved in this phenomenon. 
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Fig. 11 Impulse and thermal efficiency vs. pre-wave Mach 
number 


Since the fuel-air mixture is moving very rapidly, it may be 
able to experience temperatures which are in excess of its 
static ignition point and yet enter the detonation wave before 
undesirable effects can take place. 

The inlet static pressure ratio is a function of flight Mach 
number. Maximum pressure ratio is reached at My = 10 
and is (at 130,000 ft-altitude) on the order of 3 to 4 atm. It 
is apparent that the wave engine static pressures are lower 
than those for the conventional engine in the inlet region. 
Rapid expansion following the wave will further minimize 
the regions of the engine experiencing high static pressure 
loads. The benefits of the resultant weight reduction are 
obvious. 


Ideal Conventional Ramjet Cycle 


To form a fair comparison, we have determined the per- 
formance of an ideal conventional ramjet. The same basic 
thermodynamic data are employed. This ideal engine con- 
sists of an isentropic inlet which diffuses, in the limit, to zero 
Mach number. The combustion process is at constant 
pressure and the exhaust nozzle is isentropic. The equations 


for this ideal machine are given in the following. 
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Fig. 13 Pre-wave static temperature vs. flight Mach number 
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Fig. 14 Ideal ramjet performance characteristics 


Burner inlet temperature 


T; = mihy/Cp 


where 
= CyTofl + — 
Ts) 


Burner inlet static pressure 


4 


[14] 


= 


_ Using these equations and [7 through 10], the performance 
of the ideal conventional engine was determined. In Figs. 
9 and 10, this ideal subsonic engine follows the same fuel air 
schedule as the detonation wave engine. In Fig. 14 is shown 
the ideal subsonic engine at several different fuel-air ratios. 
The best detonation wave ramjet (DWR) is taken as that 
which has the maximum air specific impulse. 


Improved Computations 


As a check, a more elaborate inlet program was devised 
which accounts for the variation of the specific heats in the 
inlet resulting from the changes in temperature and the 
addition of fuel. In general, the results of the more exact 
program tend to yield static pressure ratios higher and tem- 
perature ratios lower than are derived from the original as- 
sumptions. Use of this improved supersonic inlet digital 
program (S. IN.) will yield a small improvement in all param- 
eters. 

As with the inlet, a more elaborated computation was de- 
rived to evaluate the expansion process. The more exact 
method again slightly improves the engine parameters but 
not to the extent achieved by the inlet program. Impetus 
to embark on both computer programs, supersonic inlet and 
supersonic expansion (S. EX.), comes from the realization 
that the engine parameters as computed here are based upon 
relatively small differences between very large numbers (1.e., 
inlet velocity and exhaust velocity) and ae potentially 
sensitive to simplifying assumptions. 


Concluding Remarks 


Our results are similar to those of Weber (7) and Dunlap 
et al. (2), although different in detail. The results pre- 
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sented in (2) are strongly influenced by their arbitrary choice 
of a fixed pre-wave temperature. Neither of the aforemen- 
tioned works considered strong detonations, since their 
actual existence was uncertain until recently. 

Basically, the subsonic burning ramjet, when compared 
with a detonation wave engine, always has superior thermo- 
dynamic performance. However, the difference in some 
instances, such as the Mach 6 flight speed regime, are minor. 
The specific fuel consumption of the detonation engine is bet- 
ter than one might initially suspect from elementary con- 
siderations (7), because, even with large inlet and wave total 
pressure losses, at high Mach number the overall cycle pres- 
sure ratio is well-matched to the temperature ratios obtain- 
able by combustion. 

The results presented here are basically for design point 
consideration. There is, however, adequate margin for ac- 
celeration if one employs the strong detonation branch. 

Experimental experience [see (10)] has shown that these 
detonation waves have a very wide stability range which is 
predictable. The chemical reaction takes place in an ex- 
tremely short time, less than about 5 microsec. Combustion 
efficiencies, which for all engineering purposes are 100 per 
cent, have been achieved in flow distances less than 1 in. 

The wave engine should be able to exploit the wave proper- 
ties so as to achieve a shorter engine than the subsonic ma- 
chine. The lower static pressures in the cycle will permit a 
lighter structure design. Also, these low pressures reduce 
the nozzle throat heat transfer, and its short length reduces 
drastically the heat transfer area. 

The designer should take advantage of the wave properties 
as they feed into the exhaust nozzle. A very rapid expansion 
nozzle may alleviate the dissociation problem and will permit 
most of the nozzle interior wall surface to ‘“‘see’”’ the outside 
world and hence cool itself by radiation. 

On the other hand, such an engine has a fuel injection prob- 
lem. The seriousness of this problem has been estimated 
by Dunlap (2) but actual experimental experience has not 
yet been obtained. Both engines have a serious variable 
geometry problem, particularly with regard to “‘off design”’ 
operation, and here the detonation machine is the more 
serious offender. The control problem is very serious in 
either case, since inlet buzz and regurgitation may be disas- 
trous at high flight speed. 

Further data may be obtained from (13), of which this 
paper is a summary. The authors wish to acknowledge 
helpful discussions with J. Yampolski, C. L. Eisen and H. Van 
Valkenburg. 


Nomenclature 


J = velocity, fps 
= area, ft? 

I = static pressure, lb/ft? 
= static temperature, R 
= ratio of specific heats, Cp/Cr 


= universal gas constant, ft lb/lb-mol R 
M = Mach number 
AH = heating value of fuel, Btu/lb mol 


moles of ith reactant 


Bi = moles of ith product 


f/a = fuel air ratio 
y = equivalence ratio, stoichiometric fuel-air ratio/fuel air — 
ratio 
mi = molecular weight of ith constituent 
W, = weight flow, lb/sec 
Q = Q/C,T 
"y = net thrust, lb 
nk = kinetic efficiency of inlet = 
= inlet pressure recovery 
Nth = thermal efficiency 
Cp = coefficient of specific heat at constant pressure — 
Q = heat added to flow, Btu/lb sec 
h = enthalpy, Btu/lb mol F 
Cp = coefficient of discharge, actual nanan ae velocity 


g = gravitational constant, 32.2 ft/sec? 4 
C,H, = hydrocarbon fuel 


Oz = oxygen 

Ne = nitrogen 

CO, = carbon dioxide 
H.O = water 
Subscripts 

0 = free stream 

1 = wave entrance 

2 = wave exit 

3 = exhaust nozzle exit 
t = total state 

7 = arbitrary constituent 
f = fuel 

I = inlet 

K = kinetic 

S = stoichiometric 
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HE GREATLY increased communication ranges and the 

antenna pointing accuracy required for present-day 
space vehicles have created a need for greatly improved rocket 
telemetry. Systems using correlation techniques based on 
recent advances in information theory appear to be the an- 
swer. However installation of such systems into rocket 
vehicles creates immediate problems. Thus, when correla- 
tion communication equipment is installed into rocket- 
powered space vehicles, the phase stability of the airborne 
units becomes a major design criterion. Although many 
available electronic components perform satisfactorily in the 
environment, it has been found that ultra high frequency 
(UHF) cavities of conventional design cannot be used. The 
essential requirements for the UHF cavities are phase sta- 
bility, light weight and high efficiency. 

A typical phase-coherent communication system uses phase 
modulation of the carrier in order to transmit intelligence 
over the radio link. A phase detector at the receiver accom- 
plishes the desired demodulation. If the transmitter is to 
operate in the microwave region, a combination of very high, 
ultra high and microwave frequency multipliers may be used 
to obtain the desired output frequency for transmitter or local 
oscillator purposes. The multiplier chain is driven from a 
highly stable VHF oscillator source so as to insure phase co- 
herence for the system. 


Phase Stability Criterion 


Mechanical vibration of components used for establishing 
RF resonance in a phase-coherent system produces incidental 
phase modulation (IPM) of the carrier or local oscillator sig- 
nals, which, if sufficiently high, may seriously impair system 
performance. With correlation detection, these IPM signals 
appear as unwanted low frequency compone nts at the output 
of the phase detector, with frequencies and amplitudes 
directly related to the rocket vibration as shown in block 
diagram form in Fig. 1. If vibration produces phase modula- 
tion which exceeds the 90-deg limits of the phase detector, 
system lock may be lost, thus, temporarily destroying the 
intelligence transmission capabilities of the system. For- 
tunately, on the premise that the rocket vibration may be 
assumed to have a spectral distribution similar to white noise 
over the band of interest, it is possible to build more sophisti. 
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Desion of UHF Multipliers for Space 
Probe Communications 


tion, an equivalent tuning stability of one part in 10°. 
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_ The basic problem of phase stability in a correlation system is discussed, and data are included to 
demonstrate the design of highly stable UHF modules. It is shown that lightweight, efficient 
UHF cavities may be built which have a desired phase stability of better than 0.1 deg per G of vibra- 


‘ated correlation systems which can extract the desired phasv- 
modulated components in the presence of high level IP\I 
Hence, a relatively large degree of IPM, but generally much 
less than 90 deg, may be tolerated without impairment of tlie 
highly accurate information such a system may normally 
carry. A figure of +30-deg peak at the phase deteetor has 
been found to be acceptable, representing a level of some 3 
deg per G for typical rocket environments. An RF component 
design, to achieve the desired phase stability, becomes most 
critical for a transmitter employing a large number of fre- 
quency multiplier stages because then, with the 30-deg cri- 
terion, the nth frequency multiplier may only have an IPM 
of 30/n deg (n is the ratio of the output frequency to the 
multiplier-module frequency). Thus, for a microwave phase- 
coherent system, with mechanical vibration of the order of 10 
g, only phase deviations of less than 0.1 deg per G can be 
allowed in the early stages of the UHF multiplier band if a 
tolerance level of 3 deg per G at the output frequency is as- 
sumed. 

The tuning stability requirements for mechanical excitation 
of this nature may be determined from the relationship be- 
tween phase angle g and the detuning of a resonant circuit 


¥ (see Fig. 2). Whe re 
RF input frequency 
cavity resonant frequency + Aw 
then 
= tan“! Q(1/y? — 1) 
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Fig. 1 Simplified coherent transmission system illustrating 
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an approximation for a high Q parallel resonant circuit, and 
Aw is the detuning caused by vibration. The cavity is assumed 
to be in precise tune with the RF input frequency without the 
vibration. 

For phase shifts of the order of 0.1 deg, which is well within 
the linear portion of the phase characteristics of a tuned cir- 
cuit 


Ag = tan~! (—2Q6) 


W here 
6 = detuning caused by mechanical vibration 
Ag = the corresponding phase excursions 


Thus, the permissible amount of cavity detuning for a Q of 
100 must be less than one part in 100,000 for a phase stability 
of 0.1 deg per G. This approaches the stability of a quartz 
crystal. The long-term frequency stability need not, of 
course, be of this order, since phase tracking circuits which 
form part of the system may readily accommodate slow 
phase changes produced by small order cavity tuning drifts. 
Expression [3] shows that the IPM will be minimized with 
lower Q. As thermal tuning drifts become serious with high 
Q cavities, a low “actively loaded” Q gives greater tuning 
stability, and, since the modulation techniques used in correla- 
tion systems may also require wide bandwidths, it is ob- 
viously advantageous to operate with the lowest possible 
louded Q in the UHF resonant cavities consistent with effi- 

‘ney. It will be shown that it is possible to design highly 
eflicient UHF multipliers which will meet the phase stability 
requirements and have the desired property of low operating 


Q 
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Payload Design Requirements 


rhe design of a UHF frequency multiplier which will sue- 
‘essfully meet both the system and the environmental speci- 
fications is complex, and such diverse problems as mechanical 
resonance with vibration, electron transit-times in vacuum 
tubes and the reduction of thermal gradients must be con- 
sidered. Desirable properties for UHF frequency multi- 
pliers to be operated in an airborne phase-coherent system 
may be outlined as follows: 

Freedom from phase modulation with vibration. 

2 Lightweight and compact overall assembly. 

3 High thermal stability and ability to dissipate heat 
directly to the payload frame. 

{ High plate efficiency and maximum RF power gain. 

5 Stability of tuning and ease of adjustment. 

6 Adequate bandwidth for the svstem under considera- 
tion. 

7 Minimum RF leakage and low external RF fields. 

S Ability to change tubes readily and to use tubes within 
the manufacturer’s range of input and output capacitance 
variation. 

9 Plate voltage requirements to be within the range of the 
normally available unpressurized low voltage airborne supply. 

10 d-e insulation to insure freedom from voltage break- 
down at high altitudes. 

It will be shown that these requirements may be met by a 
much closer integration of mechanical and electrical design 
parameters than has been customary in the past. Well- 
established operating techniques for Class C radio-frequency 
multiplier operation or the basic design of coaxial cavities 
will not be discussed except insofar as they are related to the 
design of highly efficient, phase stable UHF multipliers. 


Tube Considerations 


Overall tube considerations for a complete multiplier chain 
operating from 30 me to X-band indicate that the region 
between 200 and 2000 me must be covered by coaxial-type 
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Fig. 2 Phase characteristics of tuned circuit showing in- 
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Fig. 3 UHF triode GL-6442 showing cutaway section 


It is fortunate that the basic requirements for high 


tubes. 
UHF performance, fine wire grids and closely spaced rigidly 
supported electrodes, also produce a tube which has desirable 


phase stability characteristics. One UHF tube of this type 
which has been successfully used in missile cavities is the 
GL-6442. Fig. 3 is a cutaway section showing the internal 
structure. The manufacturer’s stated characteristics are: 
Heater power, 5.7 w; plate dissipation, 6 w maximum; trans- 
conductance, 13,500 umhos; and frequency multiplier opera- 
tion, up to 2500 me. Although the major electrodes of this 
tube are ruggedly mounted, vibration tests have disclosed 
that the heater is insufficiently secured within the cathode, 
and hence mechanical resonances cause changes in heater-to- 
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Fig. 5 Schematic of x8 UHF frequency multiplier, 250-2000 mc 
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Fig. 4 x8 UHF frequency multiplier, 250-2000 mc 


cathode capacitance. Phase modulation so produced has 
been minimized by the use of a bifilar-heater choke that 
raises the heater assembly to the same RF potential as the 


cathode. 
UHF Multiplier Module Design 


A typical UHF module for a phase-coherent system operat- 
ing between 200 and 2000 mc may be expected to have three 
cavity stages operating x3, x2, x2 or x2, x2, x2. Multiplier 
efficiency is generally too low for factors greater than 3 

Using lumped LC circuit constants instead of a tuned line 
at the lowest or input frequency, /4 cavities for the inter- 
mediate stages, and a ? \ cavity, by necessity, at the 2000-nic 
output frequency, an overall length of the order of 6 in. may 
be obtained for the three multiplier stages. Fig. 4 shows 
the overall assembly of an x8 UHF multiplier, and Fig. 5 
is the corresponding simplified schematic. 

A high degree of mechanical stability combined with ad:- 
quate conduction cooling is obtained by constructing the 
cavity bodies from rectangular blocks. Holes bored into the 
blocks form the outer conductor of the coaxial resonators. 
Having thus secured a rigid outer coaxial construction, the 

_ problems of phase stability now become directed at the inner 
coaxial assembly, i.e., the mounting of the tube and inner 
coaxial conductors, the contacts to the tube elements and 
the methods employed for cavity tuning. 

The internal construction of a cavity is determined largely 

_ by the circuit arrangements, which inevitably form part of the 
basic mechanical structure at ultra high frequency. The only 
logical RF circuit for a UHF triode is the grounded-grid con- 
figuration. By adopting an end-to-end assembly for the 
associated plate and cathode cavities, and by making these 
two sub-units completely independent, mechanically and elec- 
trically, the tube grid flange, which is at RF ground, may be 
clamped rigidly between the two independent structures with 
a strong flat spring interposed to grip the tube flange. The 
grid electrode and, thus, the tube are now rigidly clamped to 
the basic structure. Mica spacers may be interposed to 
allow for grid biasing. For the other tube electrodes, con- 
siderable success has been achieved by fabricating the contact 
springs from a solid beryllium-copper rod, slotting and re- 
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Fig. 6 1000- to 2000-mc cavity multiplier internal assembly showing construction and plate cavity tuning studs 
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moving metal as required in order to provide long, para- 
bolically curved springs which form, simultaneously, the 
inner coaxial conductors and the tube electrode contacts. 
Where the length requires additional support, a Teflon spacer, 
inserted so as to preserve the constant characteristics of the 
coaxial resonator, provides increased mechanical rigidity. 
Fig. 6 illustrates the construction of the inner coaxial assem- 
bly. 


Cavity Tuning 


Cavity tuning is required to provide desirable system fre- 
quency adjustments and to cover tube-capacitance variations. 
The relationship between electrical length 6, coaxial line 
impedance Zp, and tube or terminating capacitance C, which 
determine cavity resonance, is given by the transcendental 
equation 


tan = Xc/Zp [4] 


funing may obviously be accomplished by a variation of 
any one of the three parameters, and accepted methods of 
tuning vary either the capacitive reactance or the electrical 
length. Neither of these methods were considered suffi- 
ciently stable for the considered environment. Instead, a 
tuning method which varied the Z) was evolved. The pos- 
sible tuning range is more limited because of mechanical 
considerations, by this means, but it has the consequent ad- 
vantage of being more stable. In order to insure sufficient 
tuning range the electromechanical design is held to tight 
tolerances; then, with median capacity tubes, the cavities 
will tune to the desired median frequency. Also, the coaxial 
cavity is assumed to terminate at the tube grid flange instead 
of at the anode, and the outer cavity diameter is regulated to 
avoid the pronounced discontinuities normally associated 
with the tube. The required variation of Z) may now be 
obtained to maximum advantage by the insertion of a number 
of adjustable screws into the shorted or high current end of 
the cavity, which effectively control Z. As the mean im- 
pedance of the resonant line is changed, the electrical length 
is varied, as shown by Equation [4]; hence, the resonant fre- 
quency is varied. Tuning by this means gives accurate, 

stable and noncritical tuning control. 


The requirement for UHF cavities of high plate efficiency is 
inextricably bound up not only with the parameters enumer- 
ated in the introduction, but also with freedom from voltage 
breakdown, adequate bandwidth, maximum RF power gain 
and low plate operating potential. 

The overall efficiency of a tube-driven UHF multiplier 
plate tank circuit is 


Plate Cavity Design 


overall = [5] 
where 
n, = efficiency of the cavity 
n2 = efficiency of the tube as a Class C amplifier | 
nz; = efficiency of the tube as a Class C multiplier 
ns = efficiency of the tube as a function of frequency 
The cavity efficiency 
(Qu = Q1)/Qu [6] 


where Q, and Q, are the unloaded and loaded Q of the plate 
avity, respectively. An unloaded Q of well over 1000 may 
he expected; whereas, the desired loaded Q is of the order of 
100, which gives an efficiency of 0.90. The Class C amplifier 
efficiency mn, may be expected to be of the order of 0.75, 
i median figure for Class C operation. The multiplier effi- 
iency n3 represents the efficiency with which the tube ex- 
racts the harmonic component from the highly nonlinear 
fundamental frequency tube current. A figure of 0.65 for the 


CATHODE 


GRID 
had 
An | ANODE 
PLATE | 


| (0) ELECTRONS OF VELOCITIES 


CAVITY 
| BUNCHED AT GRID, DIFFERENT ARRIVAL 


TIMES AT PLATE 


t (b) PLATE CURRENT PULSES AT GRID PLANE ~ 


t 


f Bn BP» bm, (c) PLATE CURRENT PULSES ON ARRIVAL 
ip . AT ANODE PLANE 


Fig. 7 Effect of grid-plate transit time on shape of current 
pulse applied | to cavity showing electron debunching 


second harmonic and 0.40 for the third represent typical 
operation based on unity for Class C fundamental frequency 
operation (3). The parameter ng is a function which is de- 
pendent on the electron transit time between the grid and plate 
electrodes of the tube and is related to frequency. Bunched 
electrons at the grid plane arrive at the plate at different 
times as a result of velocity differences caused by the finite 
time of travel from cathode to grid. As a result of the dif- 
ferent arrival times, the plate current pulse is broadened (1), 
and, thus, the efficiency of multiplication is reduced. This 
electron debunching or dispersion is illustrated in Fig. 7. 
Neglecting the expected efficiency for a doubler, 
0.90, 0.75, 0.65 = 0.44. Transit time debunching may reduce 
this overall efficiency to 0.15 for a UHF doubler operating 
with an output frequency of 2000 me. This would make ng 
equal to 0.34, which at 2000 me would thus become the 
dominant term. 

With the efficiency considerations established, it can be 
seen that the external cavity bodies may be constructed, for 
example, from aluminum with considerable weight saving, 
paying only a slight penalty for its somewhat low conduc- 
tivity. 


Input and Output Coupling 


It must be assumed that for any type of Class C operation, 
the drive to the tube, the coupling and the d-c plate potential 
are all optimized. Failure to optimize coupling with result- 
ant mismatch will exact penalties not only in efficiency, but 
also in reduced bandwidth and tuning stability. Optimum 
output match is obtained when the load impedance is made 
to present an impedance equal to the effective tube impedance 
seen as a generator of harmonic components, which is 
considerably higher than when the tube is used as an ampli- 
fier. With optimum coupling, optimum bandwidth and thus 
desirable low Q for phase stability are simultaneously realized. 
The bandwidth, however, is also dependent on the terminating 
capacitance C, of the coaxial line and the line impedance Zo 
(2) as shown by the following 


fo 1 
f Q 2rC.-R: 
1 6 R, Tv 


which is an approximation for the equivalent lumped circuit. 
Therefore, in order to insure maximum bandwidth, C, must 
C, may be minimized 


be minimized and Zp) maximized. 
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is effective tube impedance as a harmonic generator. 


by making the internal plate-grid capacitance of the tube 
the only major discontinuity in the coaxial structure, and 
Zo may be maximized by making the energy storage in the 
distributed line capacitance a minimum, that is, by increasing 
the ratio of outer-to-inner conductors. 

The commonly used loop-coupling to cavities presents many 
difficulties with tandem operation of tube multipliers. A 
more effective and also simpler method of coupling is a capaci- 
tive coupling to the inner coaxial conductor. This may be 
compared to the impedance-transformation method commonly 
used at broadcast frequencies, as shown in Fig. 8. The capaci- 
tive coupling, which is more suited to UHF operation, enables a 
precise exact match to be made to the generator impedance 
and, at the same time, gives the required open circuit d-c return 
for multistage operation. A direct tap into the cathode line 
of the following multiplier is almost essential to insure an 
accurate match into the low impedance of the grounded-grid 
input circuit. Fig. 8 also shows the equivalent circuit for 
impedance transformation of R,; to match the effective im- 
pedance of the tube as a harmonic generator 


Ri’ = (Ri? + Xe,?)/Ri [9] 
with R; small as compared with 
X¢,2/Ri [10] 


Design of Input Cavity 


If Q of a cavity is sufficiently low, and the mechanical 
dimensions are accurately maintained, a relatively wide 
band of operation may be obtained without recourse to ad- 
justable tuning means. Thus, the input cavities, which are 
terminated with the low input impedance of a grounded-grid 
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Fig. 8 Matching of plate cavity to 50-ohm load using direct capacity coupling. Lc represents resonant cavity; RK, 
(a) Conventional low-frequency circuit; (b) UHF equivalent 
circuit of (a); (c) impedance transformation of load at UHF 
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tube have a low Q and can be fixed tuned. At UHF, thie 
input cavity terminating impedance also has an additional 
shunt component due to the transit time of the electrons ‘n 
the cathode-to-grid space. That is, the input admittance is 


Y, = Yr+Gr [11] 


where 
Yr = tube input admittance 
Gr = conductance term resulting from transit time loading 


For small signal inputs, the input admittance of a grounde:- 
grid amplifier, neglecting transit time loading is 


Zt) 
where 
rp = tube parameters - 
Zi = impedance in the plate circuit eh ge 


As Z_, is essentially zero for a multiplier, ¥7 would be ex- 
pected to approximate the g,, of the tube for a high p triode. 
Class C operation makes this quite inaccurate. The param- 
eter Y7 can be approximated for high level, nonlinear opera- 
tion by determining the average loading of the tube over the 
operating cycle. This may be simplified by assuming a 
linear transfer function for the e,/7, characteristics and a 
linear relationship between g,, and 7,. This is illustrated in 
Fig. 9. The plate current pulse for Class C operation may 
be regarded as a clipped sine wave, for which the average 
value from Fourier analysis is 


E (0/2) — (0/2) cos 


1 — cos (6/2) 


[13] 


where @ is the angle of current flow. The maximum g,, may 
be determined at tmax from the tube characteristics. For the 


GL - 6442, assuming a orange gm of 16, 000 wmhos and « 
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Fig. 10 Output impedance of Class C multiplier compared with 
‘small small signal grounded-grid amplifier 
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The input conductance resulting from transit time loading is 


Gr = gm(wt)?/20 [14] 


where ¢ is the transit time of the electrons from cathode to 
grid. 

For the GL-6442 at 1000 me, Gz is 2.5 X 10~* mhos, so that 
from Equation [11] 


Y, = 3.5 X 1073 + 2.5 X 10-3 mhos = 6 X 1073 mhos 


Thus, using Equation [7], with a tube input capacity of 5.5 
uuf, the expected Q will be 5.7 giving a 3-db bandwidth of 
175 me. This Q is sufficiently low to permit accurate match- 
ing over a wide band using a fixed tuned cavity. 

The output impedance of the tube for the highly nonlinear 
operation of frequency multiplication may be arrived at by 
an analysis similar to the one used for the determination of 
the input impedance. For a small signal grounded-grid 


amplifier 
Zou = Rp + + 1)R 


where R, is the equivalent input impedance, approximating 
zero for multiplier operation (see Fig. 10). The average 
tube conductance, where Gayg = 1/R, avg, for a half-sine 
wave is given by 


[16] 


based on a tube resistance which varies in a half sinusoidal 
manner over the operating cycle. The tube impedance for 
second harmonic pulse-type operation is thus approximated 


‘pave ~ p [17] 


giving an equivalent generator impedance R, of 22,000 ohms 
for the GL-6442. The tuning and coupling arrangements 
are illustrated in Fig. 11.0 


d-c Plate Voltage 


For RF power outputs of a few hundred milliwatts and 
highly efficient cavity operation, d-c plate potentials of the 
order of 150 v become possible, so that pressurization for 
high altitude operation in Earth’s atmosphere is not re- 
quired. The inefficiencies which may arise from too high a 
voltage are readily seen from Fig. 12, which illustrates that 
with a low plate voltage swing, the continuous dissipated 
power is held to a high level. However, too low a plate 
voltage may give an insufficiently high ratio of plate current 
to grid current and will also increase transit time losses 


June 1960 


MISSILE SUBFRAME 


Fig. 11 Cross section of UHF cavity multiplier 
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Fig. 12 Illustration of importance of matching plate voltage to 


RF voltage swing in order to minimize plate dissipation and in- 
crease efficiency with Class C multipliers 
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Fig. 13 Environmental performance data, vibration of x12 UHF 
cavities 
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resulting from lower electron velocities. In addition, for 
emission-limited tubes of the type used, the plate voltage 
must always be maintained in order to secure the required 
RF power output at the higher operating frequencies without 
exceeding the average emission capabilities of the planar 


cathode. 


With the tube operation adjusted for maximum efficiency, 
the power to be dissipated by heat conduction to the missile 
subframe becomes established. A steady heat flow to the 
missile subframe will insure minimum operating tempera- 
tures for the cavities. However, the required circuit d-c 
isolation of the tube elements necessitates the insertion of 
insulation materials between the high temperature electrodes 
of the tubes and the cavity body. Heat conduction between 


Thermal Considerations 


two parallel faces is determined by aii. 
Q dé 
where 
: 
- Q/T = amount of flow in a given time fps 
d0/dx = temperature gradient over the spacing x 
A = surface area in contact oe 


Obviously, minimum operating temperatures are obtained 
when contact areas are maximized and insulation kept at a 
minimum thickness. Also, for airborne application, it is 
highly desirable that the insulation should be so inserted that 
all external parts of the cavities are at low d-c potential. 
These d-c, RF and thermal insulation requirements are best 
met at UHF by discoidal-type capacitors having a plate 
diameter approaching the maximum cross-sectional dimen- 
sions of the cavity. A capacitor of this type has sufficiently 
low reactance to replace the conventional resonant line short 
circuit and thus allow the external cavity bodies to be at d-c 
ground potential. Accurate machining of all contact surfaces 
and the use of thin sheet mica as a dielectric reduce estimated 
and measured thermal gradients to a few degrees. The 
advantages of mica insulation are high heat conduction, 
mechanical rigidity in compression, high dielectric constant 
and low RF losses. Fig. 11 shows the manner in which the 
discoidal capacitors are built into the UHF cavity. 


Performance 


UHF multipliers constructed in accordance with the design 
criteria outlined here have demonstrated their ability to 
operate successfully in an airborne phase-coherent system. 


pe 


Fig. 13 shows the results of an environmental test with an 
applied vibration level of 4 g/rms. An IPM of less than 
0.5 deg per G was observed at the output frequency using « 
specially developed IPM instrumentation system. The 
wide band, low frequency peak at 150 eps is the result, almost 
entirely, of residual IPM from the heater-to-cathode vibra- 
tion of the tube. Table 1 gives the pertinent performance 
data, and the curve in Fig. 14 illustrates the input-output 
characteristics. In estimating overall efficiency, heater power 
is taken into account. Several successful flights have been 
obtained in which these UHF cavities have operated satis- 
factorily as an integral part of a phase-coherent system. 
It is not considered that the design illustrated is optimum 
either in weight or performance, but it is believed to point 
out a new direction which the design of such components must 
take if success is to be achieved in the penetration of some of 
the new space frontiers. 
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Table 1 Performance data of x8 UHF multiplier 

RF input 

RF output 

Total d-c plate input 

* Includes heater power. 
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Technical Notes 


Determination of Diffusion Coefficients 
From Asymmetric Concentration 


Profiles 


A. W. BLACKMAN?» 
& 
Research Laboratories, United Aircraft Corp., 
; East Hartford, Conn. 


N THE investigation of mixing problems which occur in 

rockets, ramjets and turbojets, e.g., those associated 
with fuel injection and combustion chamber cooling by the 
injection of cool gases adjacent to the walls, it is sometimes 
necessary to calculate the rate of change of asymmetric con- 
centration profiles when diffusion coefficients are known, or 
alternately to determine diffusion coefficients from measured 
asymmetric concentration profiles a known distance apart. 
A procedure is here developed: for the calculation of these 
quantities. 

Considering a control surface of height dy and length dz 
with unit width in the z direction, combination of the diffusion 
equations with a material balance for one of the components — 
in the stream results in the following equation 


a (af 
p (fU,)dady p ay (fU,)dxdy pd, oy 


drdy + pD, — (2 dzdy 
Ox \Ox 
Assuming that the volume of material entering the control 
surface is to equal the volume leaving during the same period, 
then 
aU, aw 
+ —= 0 [2] 
ox oy 


Neglecting diffusion in the axial direction and combining 
Equations [1 and 2] yields 
fa) 
[3] 


Neglecting the y component of velocity and replacing f by a 
new variable defined by 


F = (f — fo)/fo [4] 

Equation [3] is transformed to 

U, 

oy? 
Boundary conditions will be specified as 
at z=0 F = Fly) 
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A solution of Equation [5] satisfying these boundary condi- 
tions can be written as 


F = — > 


n=1 


This sine series will vanish at the end points y = 0 and y = h, 
but will represent /'(y) in the interval 0 < y < hif a sufficient 
number of terms is taken. F(y) at x = 0 can be determined 
from a measured concentration profile and a mass balance 
of the injected flows, and the constants A, can be found by 
solving Equation [6]. A convenient way of solving the result- 
ing set of equations is to arrange them in matrix form. Ex- 
pansion of Equation [6] at x = 0 into six terms, for example, 
would give the following set of equations for six values of F 
measured at equal intervals of y 


+ +... Asyis = Fi 


[7] 


+ 
which can be written in matrix form as 


yie F, A, 


Yu 
[8] 

Yor Yer. Fe Ag 

The inverse matrix can be determined by the method of 
triangulation? and the constant A;, Az... Ags can be deter- 
mined for different initial profile shapes by multiplying the 
matrices, provided the measured values are always taken at 
the same points. 

Diffusion coefficients can be determined by solving Equa- 
tion [6] utilizing values of F measured at two stations a 
known distance apart: or alternately, if the diffusion co- 
efficient is known, Equation [6] can be solved for the rate of 
change of a known initial concentration distribution with 
distance. 


Nomenclature 

A = constant 

D = diffusivity 

f = weight fraction of diffusing species — 


h = height of duct in y direction 
U = velocity 


zx = variable axial distance 

y = variable vertical distance 
p = weight density of stream 
Subscripts 

x = 2x direction 

y = y direction 


2 Dwyer, P. S., ‘‘Linear Computations,’’ John Wiley and Sons, 
N. Y., 1951, p. 183. 
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Angular Displacement of a Magnetized 
Sphere by the Terrestrial Magnetic Field 


During Vertical Descent 


Ss The angular displacement bounds of a sphere with a 
magnetic dipole field are given for free-fall vertical descent 
in Earth’s gravitational and magnetic fields. Effects of 
the atmosphere are neglected, and Earth is assumed 
spherical. The terrestrial magnetic field is treated as the 
field from a dipole situated at Earth’s center. 


MARVIN H. HEWITT! 


The Martin Co., Denver, Colo. 


SPHERE in free-fall vertical descent has a magnetic 
dipole field. The torque acting on this sphere in 
Earth’s magnetic field is given by 
L = —MH sin B [1] 
where the minus sign indicates that the torque is in such a 
direction as to decrease 8, the angle between the axis of the 
dipole and Earth’s field and 


LL = the torque on the sphere 
the sphere’s dipole moment _ 


M = 
= Earth’s magnetic field 


The sphere’s center of gravity is taken coincident with its 
geometric center and the dipole is centered there. Let the 
magnetic axis of the sphere be initially along the gravitational 
plumb line. The Earth’s dipole field has moment a and is con- 
sidered to be situated at Earth’s origin and parallel to the 


axis through the poles. 
Let 

~ Ay r = central distance to sphere 

@ = (1+ 3sin? 

6 = latitude q 
Consequently 
L = —(Mad/r?) sin B (2 


Also, let a be the angle between the line of magnetic force and 
the plumb line. We obtain for vertical descent 


sin a = cos 0/d = constant [3] 
Recalling the definition of 8, we obtain 

a-—-B=6 [4] 
where 6 is the angle between the sphere’s magnetic axis and 


the plumb line. Therefore a 
[5] 


— (Ma@¢/r) sin(a —6) [6] 
pone 


sin (a — 6) = sin B 


The sphere’s momentum A y ields” 


dA/dt = [7] 
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Va 


dt Ir 


where the point 0 about which the moment is taken is at the 
center of gravity of the sphere. Consequently = 


(d/dt)Ilw = L [8] 
where 
I = constant moment of inertia of the sphere | 
w = the body’s angular speed 


a sin (a — 6) 


Moreover 


Equation [10] yields 


$= sin (a — 6)dt? 


where the integration is over the time of descent. 
It is apparent that the maximum value of | @sin = — 6)/r| 
isgivenby 


abi iw 


| — RR 
where R is Earth’s radius. 


[12] 


Also, t here is a continuous variable 


of limited total variation; therefore a4. 
5 < Mat,2/IR* 

Initial values of w and 6 have been taken equal to zero. The 


moment of inertia of a sphere about a diameter is 


= (2/5)mb? [14] 


We obtain the bound 
[15] 


If an iron sphere with a magnetic moment of 3 X 10 
gauss-cm* and a radius of 50 cm, fell for 10? see in Earth’s 
dipole field (8.06 X 10% gauss-em*), 6 would be less than 0.2 
deg. 

Deflection 6 can be given a sharper bound than that given 
by [15] if a — 6 is small (large latitudes). Then, since sin 


where } is the radius of the sphere. 


6 < 5Mat,2/2mb?R3 


(a — 6) is a weak function of time or 7 
Ma(a — 6) 
6 = ol ff $ dt? [16] 
We obtain 
5 < Mal (a — [17] 
For a sphere [17] becomes 
6 < 5Ma(a — [18] 
Inequality [18] may be put in the form 
a(5Mat [19] 


1 + (5Mat4?/2mb*R?) 
where the bound on 6 is exhibited as a function of latitude 
through a. It is apparent that bound [19] has sharpened 
the range of 6 under, of course, the restriction on latitude 
that has been ee 
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Calculation of Electrical Conductivity 


off Tonized Gases 


4. SHERMAN! 


General Electric Co., Flight Propulsion Laboratory, 
Cincinnati, Ohio 


N ANY practical application of magnetohydrodynamic 

concepts a knowledge of the electrical conductivity of the 
working fluid is essential. Unfortunately, however, for an 
ionized gas of varying pressure and temperature the evalua- 
tion of this property by experiment or analysis is, at best, a 
difficult task. Some experimental measurements have been 
mace (1,2,3),? but, in general, they are not easily carried out, 
since there is no readily available method of containing uni- 
form samples of extremely high temperature gases for reason- 
able lengths of time. In order to avoid this experimental 
problem, at least temporarily, it would be very desirable to 
have some accurate analytical method for electrical conduc- 
tivity calculations. The purpose of this note is to review 
existing methods and call attention to a general calculation 
procedure. 

A large number of theoretical investigations presented in 
recent years deal with the question of electrical conductivity 
of ionized gases, but none treat the problem completely. 
That is, a number of authors have derived equations for 
conductivity assuming a slightly ionized gas, and others 
have done the same for a completely ionized gas. There 
have been no basic investigations of the conduc tivity of a 
gas with an arbitrary degree of ionization. 

A careful analysis of conductivity for a slightly ionized 
gas is given by Chapman and Cowling (4). They consider a 
binary mixture of neutral particles and electrons and show 
the electrical conductivity to be oe 


= (ne?/kT)Dy 


where 
ne = electron number density .) 
e = electronic charge 
7’ = absolute temperature > 
k = Boltzmann constant Py 
Dy» = diffusion coefficient 


The diffusion coefficient is derived for a binary mixture of 
rigid elastic spheres in which one species is much lighter and 
much less plentiful than the other. It is found to be 


This equation is valid for a gas which is so little ionized that the 
influence of the ions is small, and the ions may be counted as 
neutrals. That is, the ion number per unit volume is so small 
that an electron rarely collides with one, and such collisions 
are unimportant. 

Two different theoretical approaches have been used in the 
derivation of formulas for the electrical conductivity of a 
fully ionized gas. The first assumes the Boltzmann equation 
to be valid and solves it by the Enskog procedure. Using this 
method a number of calculations have been made (5,6,7). A 
result typical of the rest is found by Burgers (7) and is given 


(kT)? 


= 0.582 In A 


= . 


where 


3 1/2 
2e* \ rn, 


On the other hand, when the forces between the components _ 
of the gas are Coulomb forces, it becomes harder to justify the — 
basic assumption underlying the Boltzmann equation. That 
is, the ratio of the time spent during a collision to the time 
spent between collisions is very small. In a fully ionized gas, a 
charged particle will experience many small deflections rather 
than several sudden collisions. As a result a fresh approach 
to the problem was developed by Spitzer (8,9), and has been 
used by him as well as Allis (10) to derive an equation for 
electrical conductivity of fully ionized gases. A summary of 
their results is given ; 


(kT)? 1 


In A 


o = 0.591 


where A takes on different values. From (8) we have 


31 
A, = 
2 e? \ rn, 
and from (9) we find 


1 31 A 
533 (=) 
22 


TNe 


finally, from Allis (10) 


31 1/2 A, 


The equation derived by Spitzer (8) is exactly that given by 
Burger except for the constant. 


3.396" [ kT 

Dy = SroQ EA ¥ [2] It is of interest to note that Equations [4] or [5] are inde- 
pendent of the nature of the gas. In addition, the conduc- 
where tivity depends weakly on n,, through the logarithmic term, 
1 and strongly on the temperature. Actually, if 7’ is high enough 
n = original particle number density = (n. + n,) 7 to have a fully ionized gas, then o will increase with gas pres- 
n, = neutral number density wl 7 “a 7 sure. At lower temperatures the reverse is true, since a, the 
= collision cross section = 10,2" = =f degree of ionization, depends very strongly on pressure and 

712 = radius of collision cylinder -— inereases with decreasing pressure. 
me = electron mass One approximate procedure for calculating the electrical 
Combining this with Equation [1] gives Chapman and conductivity of a gas of arbitrary degree of ionization has 
Cowling’s formula been developed by Lin, Resler and Kantrowitz (1). They 
state that the resistance (electrical) due to neutral particles, 
ea ae? 1 3 and the resistance due to ionized particles, should be additive. 

(m Q Thus they write 
o On 
a = degree of ionization = n,/no ite 
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and using Equation [3] for o, while taking Equation [5] for 
go; with A from (9) this becomes 


= 1.69 ——— A 


le 
| 
* 
a 
| 


where 
3 (== 
A \ an, ) 


and Q, refers to the collision cross section for collision of elec- _ 


trons with neutral atoms or molecules. 
Genera] Calculation Procedure 


Another approach to obtaining a formula valid for arbitrary 
degrees of ionization, the so-called ‘free path kinetic theory,” 
has been employed by a number of authors. The basis of this 
theory is that the particles of a gas are assumed to make in- 
stantaneous collisions with each other but move freely be- 
tween collisions. This differs from the more refined ap- 
proaches in that it does not inquire into the actual form of the 
forces acting between the particles (ion, electron, neutral 
atom). However, it does have the virtue of greater generality. 

If the electron mean free time 7 is assumed to be constant 
for all electrons under given conditions of pressure and tem- 
perature, the following formula can be derived (11) 


o = (n,e?/m,)T [8] 
where 7 is given by 
C (QT /rm)" 
A = electron mean free path a 
e = average electron velocity suitable for a Maxwellian 
distribution [for 7 = constant, it can be shown 


(11) that the distribution function for electrons is 
indeed Maxwellian ] 


In the development of Equation [8] the ionic current, pressure 
and temperature diffusion phenomena, and the influence of a 
magnetic field have been neglected. These could be included 
but really constitute a separate investigation. Finally, in 
order to bring Equation [8] into agreement with Equation [3] 
for a slightly ionized gas we write 


o = 0.85(n.e2/m,)r [10] 


Before use can be made of Equation [10] we need some way 
of evaluating for gas mixtures. An expression for this 
quantity is derived by Jeans (12) from free path theory and is 


A= 1/20; {11] 


where 
n; = number density of neutrals or ions 
Q; = cross section for collision of electrons with neutrals or 


ions 


As far as neutral species are concerned, considerable experi- 
mental data are available and are summarized by Massey and 
Burhop (13). A collision cross section for electrons and ions 
can be obtained if we consider a fully ionized gas and compare 
our present equation to Equation [5] with A obtained from 
(9). The justification for this choice for A is the experimental 
evidence with argon (1). Making the comparison gives 


\2 15 
12 
(=) (=) ] 


If the collision cross sections for neutrals are obtained from 
Ramsauer-type experiments, then a suitable average cross 
section is given by (1). 
kT 


= 
0 n 


where C is the electron speed. When the cross sections are 
obtained from electron diffusion experiments, this averaging 
is unnecessary. 

Next, in order to calculate o we have to calculate n,. For 
the special case of a monatomic gas at temperatures low 
enough so that electronic excitation above the ground level 
may be neglected, the Saha equation may be used to obtain 
n,. For other situations more complex calculations are necvs- 
sary. 

One final point would be worth discussing. It is often stated 
that the terms in A, Equation [11], due to electrons colliding 


with neutral particles (Dinn:Qni) can be neglected com- 


pared with the electron-ion term (n,.Q;) if the degree of ioniza- 
tion is only a per cent or so. As an example of one case in 
which such a statement is incorrect we consider mercury 
vapor (14) at l-atm pressure. Jn this case, we find that thie 
two terms in X are equal at 7400 K and 2.7 per cent degree of 
ionization. However, the contribution due to ionic interac- 
tions is larger by a factor of 10 only at 10,200 K and 38.3 per 
cent degree of ionization. Thus, if the relation for a fully 
ionized gas, Equation [5], has been used to compute the con- 
ductivity of mercury vapor at 1 atm and 7400 K, the calcu- 
lated value would be too large by a factor of 2. The assump- 
tion which we have just shown to be in error in general is, 
however, more accurate when the gas being considered ionizes 
at a lower temperature and when its electron neutral cross 
section is small. In other words, lower pressures and/or lower 
ionization potentials tend to make such an approximation 
more valid. 
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Approximation of Optimum Thrust of the line segments MQ and QE can be made to obtain closer 
tangent points. The parabola thus constructed yields a 


Nozzle Contour close approximation to the contour computed by the method 
of characteristics. For example, in the case of a nozzle of 
area ratio 25:1 and length ratio 12:1, the parabolic approxi- 
Rocketdyne Div., North American Aviation, Inc., mation does not differ from the computed contour (3) by more 

Canoga Park, Calif. than 3 per cent in radial dimension. The maximum devia- 
tion in this example occurs at one third of the distance from 
throat to the exit. 

The values of 04 and 8z shown in Fig. 2 are the results of 
calculations using y = 1.23 and throat conditions such as 
shown in Fig. 1. A different value of y does not appear 
to appreciably change the nozzle contour when the area ratio 
and length are prescribed (3). Similarly, one can expect that 
minor changes in the throat region would not cause large 
changes in the contour. Hence, the values shown in Fig. 2 
and the parabolic contours can be used in general to approxi- 
mately describe the optimum thrust nozzle contours. had 


G. V. R. RAO! 


)p- A prescribed length and area ratio the nozzle contour 
' producing maximum thrust can be uniquely deter- 
mined (1,2).2. The thrust performance of such optimum con- 
toured nozzles, for various area ratios and lengths, was dis- 
cussed in (8). Even though such nozzle contour is unique, 
determination of it is by the method of characteristics and 
requires the use of high speed computing machines (1,2). 
Besides the nozzle thrust, a design engineer would like to 
know its weight and cooling requirements to enable him to 
make the proper choice. For this purpose an approximation 
to the optimum thrust contour is enough, and a simple geo- a 


m: trie method is described here. 1 Guderley, G. and Hantsch, E., ‘“‘Beste Formen fiir Achsensymme- 


The first step in defining a nozzle contour is the description trische Uberschallschubdusen,”’ Zeitschrift fiir Flugwissenschaften, vol. 3, 
of the throat region. After a certain amount of initial ex- no. 9, Sept. 1955, pp. 305-313. 
ribed th ll. the cont 2 Rao, G. V. R., ‘‘Exhaust Nozzle Contour for Optimum Thrust,” Jer 
pansion along the prescribe t roat wall, the con our 18 SO PRoPULSION, vol. 28, no. 6, June 1958, pp. 377-382. 
formed as to turn the expanding gases to a near axial direc- 3 Rao, G. V. R., “Optimum Thrust Performance of Contoured Nozzle,” 
tion. For the family of nozzle contours considered here and 


in (3), the throat region is assumed given by two circular arcs 
as shown in Fig. 1. With these initial throat conditions, 
let us consider a nozzle having a certain area ratio A./A; 
and length. ratio L/R, The nozzle contour optimized for 
thrust will have unique values for 04, the wall slope at the in- 
flection point M, and @z, the wall slope at the end point E. 
These values of 6% and 6z can be evaluated directly (3) prior 
to computing the contour itself, and are shown in a para-— 
metric form in Fig. 2. For the sake of clarity, even values — 
for 64 and odd values for @z are indicated in the figure. 
The nozzle end point given by (Rz/R,, L/R,) and the respec- 
tive values of 64 and 6z are sufficient to describe a parabola. — 
In Fig. 1 is shown a simple geometric construction for the : 
parabolic contour. Through the nozzle end point £, line 

EQ is drawn inclined at an angle 6g to the nozzle axis. Line 
MQ is drawn inclined at an angle 04 to the nozzle axis and 
tangent to the prescribed wall profile in the throat region. 
The line segments MQ and QE are then divided into an 
equal number of parts by points a, 6, ¢ and e, f, g, respec- 
tively. The required parabola is the envelope of lines join- 
ingatoe,btofandctog. If necessary, further subdivision 
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Fig. 1 Sketch of a nozzle showing nomenclature and construc- Fig. 2 Parametric values of 9y and 9; corresponding to various 
tion of the parabola exit radius ratios and nozzle length ratios 
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Estimating Transient Temperature 
Distributions During Ablation 


M. RICHARD DENISON! 
Aeronutronic Diy., Ford Motor Co., Newport Beach, Calif. 


T IS often important to investigate the penetration of heat 

in a solid whose surface experiences ablation. Because 
of the removal of material, a simple transient heat conduc- 
tion analysis may lead to over-conservative estimates of the 
penetration distance. Integral methods similar to those of 
boundary layer theory have been used recently to solve the 
steady two-dimensional conservation equations for an ablat- 
ing material with a liquid layer (1, 2).2 Also, these methods 
have been applied to transient one-dimensional heat con- 
duction (3-5). In this note it is shown that the integral 
technique is suitable for studying the combination of these 
effects. In addition, the heat penetration distance is com- 
puted for a simple special case which may be useful for design 
purposes. 

The important physical regions which take part in the abla- 
tion process are illustrated in Fig. 1. In this note, attention 
is focused on the condensed phases. When the liquid layer 
is present, the thermodynamic change of phase from solid to 
liquid may play a minor role in condensed-phase heat and 
mass transfer phenomena. Instead the liquid layer may be 
defined as a region in which the velocity parallel to the surface 
decreases, due to the strong increase of liquid viscosity with 
temperature drop, from its value at the interface u; to zero. 
The normal mass flux relative to the surface at the cooler 
boundary of the liquid layer is ms. The difference between 
ms and the mass flux m; at the gas-liquid interface represents 
a liquid runoff rate. The thermal layer extends from the 
liquid-gas interface where the temperature is 7; to the un- 
disturbed solid where the temperature is 7). Attention is 
restricted to cases in which the heat does not penetrate to 
the back surface. Chemical reactions and heat of fusion 
in the condensed phase are neglected. 

Integral equations for conservation of mass and energy on 
two-dimensional or axisymmetric bodies are given by 


-3 
udy + (ms — = 

T-T 2 rf,” Td 


— + = 0 [2] 


Conservation of momentum is represented by the assumption 
that the shear is constant across the liquid layer. For several 
materials the variation of viscosity with temperature can be 
represented by 


If, in addition, the liquid layer 
6 is considerably thinner than the thermal layer 67, a linear 
temperature distribution in the liquid layer is adequate. 
Under these circumstances a velocity distribution with respect 
to temperature in the liquid layer can be ise wal 


where 8 and A are constants. 
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Successive integration by parts results in the asymptotic e\- 
pansion 
TK 


Mids n=i - 


For purposes of describing the temperature distribution in 
the thermal layer, similar profiles which satisfy the boundary 
conditions may be applicable. In this case , 


sy 
T; — To f 
Use of the distributions of Equations [5 and 6] in the con- 


servation Equations [1 and 2], leads to the following ex- 
pressions 


— ty =~ — 
ms — m; (T;) [7] 


Gi = — msTo) + cT (ms — + 


(KT. 
(EE! F(T) Ts) + Te)br [8] 


N 
FAT;) = (F(T) > (—1)"+ ni (n +m + 2)! 


n=l m=1 


where 


n! n!(n +m + 1)! 1)! 
> —])ntm 
n=1 m=1 (n +4 2)! (n 


In general, Equations [7 and 8] must be matched to solu- 
tions of the gas-phase boundary layer equations. If A/T; 
is sufficiently large, F.(7;) is approximately unity. This 
condition applies fairly well, for example, to the ablation of 
silica. The third term on the right-hand side of Equation 
[8], which involves the x gradient of temperature, is zero at 2 


ee stagnation point and sometimes may be neglected in other 
‘regions. 


Then Equation [8] becomes identical to that which 
4 would be obtained if the heat convected parallel to the sur- 
face in the liquid layer were neglected. Under these circum- 
stances, if 7’; and ms are known functions of time, Equation 
{8] alone is sufficient for estimating the depth of penetra- 
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tion é6r as a function of time. Because of the similarity 
assumption of Equation [6], the heat flux qg; is related to the 
surface temperature and 67 through 


= a2K(T; — [9] 


where a = f’(0). 
If a simple parabolic temperature distribution is assumed, 
the boundary conditions are satisfied by 


T—T% y y 
(+) @ (10] 


Therefore 


Sometimes the surface temperature 7; quickly rises toa = 


rither constant value determined by vaporization conditions 
If, in addition, the mass flux ms is constant, Equation [8] can 
be solved analytically. The result is 


aK ark araipK 


For small values of time or low ablation rate, retention of 
only the leading term in a series expansion of the left-hand 
side of Equation [12] and use of Equation [9] for a parabolic 
temperature distribution results in 


= 


The exact solution yields (6) under these conditions . 
A 
Gi Kpc 
= 14 
— To V Tt 
In the steady state, Equation [12] becomes 
msc(T; — To)/qi = 1 [15] 


Equation [15] is identical with the exact solution, whereas 
Equation [13] differs only by V/3/m. A comparison be- 
tween the approximate and exact solutions for the entire range 
of time is shown in Fig. 2. It can be seen that the maximum 
deviation is on the order of 10 per cent. In order to get an 
indication of the adequacy of the thermal thickness 67 
as determined by the approximate method, the temperature 
rise at the location 67 was computed from the exact solution. 
The result is also shown in Fig. 2. The temperature rise 
varies from about 2 to 13 per cent. Therefore, in the case 
studied, the simple parabolic distribution for determining 67 
is fairly good. However, when similarity is not valid, e.g., 
in cases in which a material is first heated then cooled, a modi- 
fication of the method should be developed. 

Some general conclusions can be drawn from Fig. 2 which 
may be useful for design purposes. The dimensionless 
parameter mg%ct/pK plays a role analogous to that of the 
Fourier number for transient heat conduction in a finite slab. 
If ms2ct/pK > 3, then steady-state ablation conditions hold, 
and 67 is approximately equal to2K/msc. On the other hand, 
if ms2ct/pK < 0.01, ablation does not affect the temperature 


distribution, and 67 is approximately equal to 20/3Kt/pe. 


ie 
Lo STEADY STATE 7 
APPROXIMATE 
(Eqn.i2; 
0.8} 


EXACT C(Ti- To) 
(Ref. 6) 
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Nomenclature 


c = specific heat 

f = temperature distribution function (see Eq. [6] ) 

k = unity for axisymmetric and zero for two-dimensional bodies 
K = thermal conductivity 

m = mass flux per unit surface area 

q = K(0T'/dy), heat flux 

ro = distance from axis of symmetry to body surface 

t = time 

7 = temperature 

To = undisturbed temperature 

u = velocity parallel to surface 

x = distance along surface 

y = distance normal to surface in outward sibieiaen’* _ 


= thickness of liquid layer 


67 = thickness of thermal layer 
= viscosity | 

= A/T (see Eq. [3]) 

p = density | 


= pw(du/dy), shear stress 


Subscripts 


i = gas-liquid interface condition 
6 = condition at inner edge of liquid layer 
5p = condition of inner edge of thermal layer _ a 
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Isentropic Compression of Shock Tube 
Driver Gas 

<a 

‘be a National Research Council, Ottawa, Canada 3% 


R. J. STALKER? 


SING the constant area shock tube, with air as a test 
gas, it has usually been found convenient to heat the 
driver gas in order to substantially exceed a shock Mach 
number of 12 (1).2. This heating has generally been effected 
by combustion or by electrical means, but may, of course, re- 
sult from isentropic compression immediately prior to initia- 
tion of the shock tube flow. The present note is intended to 
draw attention to a method of doing this, and considers the 
resulting shock tube performance when helium is used as 
driver gas. 

The apparatus involved is represented schematically in 
Fig. 1. The method of compression, following a technique 
used in fundamental physical studies of gases at high pres- 
sures and temperatures (2), employs a free piston driven by 
air pressure into a cylinder containing the gas to be com- 
pressed. The piston is released at station A, accelerates, and 
comes to rest again at station B when all the energy imparted 
to the piston by the air from the reservoir has passed into the 
compressed gas. If the volume of the reservoir is sufficient 
for the pressure therein to change negligibly during the com- 
pression process, and the maximum velocity of the piston is so 
low that not only is the gas compressed isentropically, but 
the pressure driving the piston remains approximately equal 
to the reservoir pressure, then 


For large compression ratios p4/ps < 1 and therefore wi 


= -1=(—) =]1 —-1)—h-i 


which enables the state of the driver gas at peak pressure to 
be determined from its initial state and the initial pressure 
ratio across the piston. Arranging so that the diaphragm rup- 
tures near the peak pressure, and that the mass of the piston 
is such that the time rate of pressure variation is small near 
this value, the subsequent unsteady flow may be taken to be 
identical with that in a conventional shock tube, and the per- 
formance computed on this basis. 

Results of such computations, based on perfect gas assump- 
tions, with helium as driver and air as test gas, are shown in 
Fig. 2. Two intersecting sets of curves are presented; one 
set displays the variation of the velocity behind the shock 
wave with the overall pressure ratio Pr/P; for differing 
volumetric compression ratios p/p, and the other shows the 
same variation for selected values of the maximum diaphragm 
pressure ratio Pp/P,. 

The curves indicate that, using this technique, high shock 
Mach numbers may be obtained. However, it must be 
nee out that in an apparatus of practical dimensions, the 

rge ratio between the volume of the compression cylinder 
and that of the driver section of the shock tube implies that 
the driver section must be short, and therefore that the as- 
sociated test time of steady flow between the passage of the 
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shock wave and the contact surface is limited. For example, 
using a 3-in. diam compression tube 10-ft long, with a 
volumetric compression ratio of 75 and a 1-in. diam shock 
tube, the theoretically available test time at M, = 25 is 
approximately 50 microsec, when the air in the low pressure 
section is at 5mm Hg and allowance is made for real gas effects 


on the density ratio a across the shock wave. 


Nomenclature 
P = pressure 
p = density 
a = speed of sound 
y = ratio of specific heats in driver gas 
Per = reservoir pressure 
P, = initial pressure of test gas in mane tube 
M; = shock Mach number 
Subscripts 
A = conditions in driver gas before 
B- = conditions in driver gas at peak pressure 
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An Approximate Analytic Solution of 


Re-Entry Trajectory With Aero- 


dynamic Forces’ 
Polytechnic Institute of Brooklyn, Freeport, N. Y. 


KENNETH WANG? and LU TING: 


HE PROBLEM of re-entry trajectory for space vehicles 

with aerodynamic forces has been investigated by Lees, 
Hartwig and Cohen.‘ In Lees’ paper an analytic solution is 
obtuined when the entry angle is small and the entry velocity 
is close to the circular satellite velocity (26,000 fps). 

In this note, the method of Lees et al. is modified so that an 
approximate analytic solution is obtained without the restric- 
tion on the entry velocity. Numerical results are computed 
for entry velocity equal to 35,000 fps and compare favorably 
with the machine solution presented in Lees’ paper. 

Most of the nomenclature used in this paper is the same as 
in Lees et al. 

The equations of motion in the directions tangential and 
normal to the trajectory, asshownin Fig.1,are 


[1] 
m 


In deriving these equations, g is assumed constant and the 
velocity of air due to the rotation of the Earth is neglected. 
For most entry trajectories, 3 is small, and these equations 


may be simplified 
dV 2 


dd L v2 V°CLA yt 


2m r 


Terms in the parentheses in Equation [4] were neglected in the 
analysis of Lees et al., for the case where the entry velocity 
is close to the circular satellite velocity. 

Equations [3 and 4] yield 
2m dV 


[5] 


From the exponential approximate atmosphere p = po exp 
(—Gh), the following relation is obtained by differentiation 
with respect to time 


dp/dt = BpVd 
With Equation [3] it follows 


The change in velocity from the entry to the peak acceleration 
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has been shown to be moderate.‘ Furthermore, the quanti- 
ties in parentheses in Equation [8] represent the terms due to 
gravitation and centrifugal force. They are predominant 
only at the earlier stage of entry and become relatively in- 
significant as compared to the lift term as the velocity de- 
creases. Therefore, a good approximation is achieved by re- 
placing the velocity V in the gravitation term with the entry 
velocity V.. Since r remains practically constant, Equation 
[8] may now be integrated with Cz held constant. 


v= [». mB (p — Pe) (2 ze) In 4, [9] 


For C, varying as a function of p, Equation [9] may be 
modified. The acceleration vector felt by the pilot, in gravita- 
tional units g, is given by Lees et al. as 


pV°C,A _| 
= — 4/1 — 10 
2mg y + Co [10] 


Thus, the condition for the peak acceleration or G,, for con- 
stant Cp is 


Om = (CrA/mB)pm (11) 


Equation [9] together with Equation [11] gives the density pm 
and the angle 3, at the peak acceleration. Equation [7] may 
be integrated using Equation [9]. At the peak acceleration 
it becomes 


om dp 


= 12 
Vn 2mB Pe 


Note here that the integrand has #(p) as the denominator. 
The integral depends predominantly on the minimum value 


of 3(p) which is #,, at p = pm. Hence, the term In (p/pe) in 
in Equation [12] may be evaluated aa 
Vn  2mBV/C 


3(p) or Equation [9] is expanded in terms of p/pm. By re- 
taining only the first two terms in the expansion, the integral 
Ve 
In — = 
[<4 + Bom + Con? + + 13) 
VA + Bp. + Cp? + + (B/2VC) 


4 Lees, L., Hartwig, F. W. and Cohen, C. B., ‘‘The’ Use of Aero- 
dynamic Lift During Entry Into the Earth’s Atmosphere,” 
ARS JourNAL, vol. 29, no. 2, Sept. 1959, pp. 633-641. 


aa 
x 
A 
and Equation [5] be J 
es ENTRY TRAJECTORY 
@ 
r 
O x 
I 
Fi 
ig. 1 Inertia’ 
1 coordinate system ae 
565 


C./Cp 
30 
O AND RESULTS OF THE APPROX. ANALYTIC } 
SOLUTION 
° —— RESULTS OF THE MACHINE CALCULATION fo) 10 
IN REFERENCE | 
=x 
w 
a 
4 
i?) 3 6 9 12 


ENTRY ANGLE, @¢ (DEG) AT 400,000 FT. 
Fig. 2 Peak acceleration with constant C,/Cp. Ve = 35,000 


i g Pm 


where 


A 


ll 


26 
O AND A RESULTS OF APPROX. ANALYTIC 
SOLUTION 
24+ —RESULT OF THE MACHINE 


CALCULATION IN REFERENCE | 


3 6 9 12 
ENTRY ANGLE, @¢ (DEG) AT 400,000 FT. 


22 


VELOCITY AT MAXIMUM DECELERATION (Vm) x FT/SEC 


Fig. 3 Velocity at peak acceleration. V,. = 35,000 - 


In Figs. 2 and 3, the values of the peak acceleration and 
the velocity, respectively, for C,/Cp = (1/4) and 1 are com- 


pared with those in Lees’ paper® for the case of V. = 35,000 
fps. 
Applications and generalizations of present analytic solu- 
tion will be presented in a later report. 7 


5 These data are obtained from an ASME preprint, no. 59-V-36 
of the forementioned paper by Lees, Hartwig and Cohen. 


Pressure-Density-Temperature 


EUGENE A. TKACHENKO! 
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Two methods, one based on thermodynamics the other 
on a generalized equation of state, are employed in con- 
junction to calculate the pressure-density-temperature 
relationship of liquid oxygen in the temperature range of 
110 to 210 R, and in the pressure range of saturation con- 
ditions to 2500 psia. Since no compressibility data for 
liquid oxygen are available at temperatures higher than 
160 R, the results of this investigation should be of special 
use in rocket engine applications, where the temperature 
interval 160 to 210 R is of greatest interest. The pressure- 
density-temperature relationship of liquid oxygen is pre- 
sented in tabular form, whereas the coefficient of isother- 
mal compressibility is shown graphically, and is compared 
with previous calculated values. 


HE EXTENSIVE use of liquid oxygen as a rocket pro- 
pellant places special emphasis on knowledge of its 
volumetric properties. To date, no measurements of the 
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Relationship of Liquid O, 


pressure-density-temperature relationship in the liquid 
region of oxygen have been made except at saturation condi- 
tions, although a definite need exists for such data for 
large temperature and pressure ranges. In order to furnish 
liquid oxygen volumetric data for rocket applications which 
require better accuracy than the assumption of an incompres- 
sible liquid, certain basic thermodynamic relationships to- 
gether with a generalized equation of state were employed, 
resulting in a refinement and extension of previous analytical 
results (3).? 

The two methods which were employed conjunctly may 
be described in short as follows. One method requires 
saturation liquid density, vapor pressure, liquid enthalpy at 
saturation conditions and sonic velocity as functions of tem- 
perature. By means of basic thermodynamic relationships 
the coefficient of isothermal compressibility is obtained from 
these properties. The pressure-density-temperature relation- 
ship, in turn, may be computed on the assumption that the 
coefficient of isothermal compressibility is independent of 
pressure, or on the basis that it varies with pressure charac- 
teristic of normal liquids. The second method requires the 
following minimum number of physical properties: The 
three critical constants, the normal boiling point and the 
density of the saturated liquid at two temperatures. This 
method is formulated in an equation of state of the Van der 
Waals type based on a modification of the principle of cor- 
responding states, and was developed by Hirschfelder and co- 
workers (6). This equation of state is semi-empirical in 


2 Numbers in parentheses indicate References at end of paper. 
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nature, with constants which can be determined from gen- 
eralization of empirical data for many substances, or from 
experimental properties of the substance studied. 


Method 1 

Utilizing certain thermodynamic relationships it is possible 
to arrive at isobaric and isothermal conditions from satura- 
tion states 


D = f(T, P) 


dD = (2) dT +( 
\or]p 
(2 
D \oT /p 
dP 
a= Br aT = 
| H = f(T, P) 
oH oH 
aH = dT ——} dP 
(2 (2), 
oH oH 
Cp = = — 


From thermodynamics 
Ta? 
Br = += [3] 


De? CpD 
From Equations [1, 2 and 3] the value of 87 is 


Br 


_ — — (1 + 
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Fig. 1 Liquid O. compressibility 


(dH/dT) — (1/D)(dP/dT) + (T/D*)(dP/dT)(dD/dT) — (T/D*c?)(dP/dT )? 


It is possible to arrive at a pressure-density-temperature 
relationship by assuming that 


Br # f(P) 


L (22) 


and integrating at constant temperature 


AD = D,(e®T4P — 1) [5] 


If the substance studied can be considered a normal liquid 
the following relationship (5) can be employed pane 


(01ln Br/OP)r = —10Br 


Integrating at constant 7 
Br = (10P + 
Substituting Equation [7] in Equation [4] and integrating 


AD = D,{(10P2 + C/10P; + — 1) [8] 


Method 2 

The development and some applications of Method 2 
are described in (6 and 7). Rather than recapitulate the 
details of Method 2, only essential equations and comments 
will be presented 


Pi = —[ko + (B + — [(1/2)(1 — ko — + 28)(1 — + 1 + Bp? + (8.44 + 4.508 — 0.3636?)(p — 1)%(1/p) + 
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pi(p, t) = pr(p, t) — prlpr(t),t] + [10] 
where subscript I signifies the dense gas region and II de- 
scribes the liquid region, and 

B(3B — 1)\(1 — = (P./RT.D.) = Z. [11] 
ko, B, a, ho, hi, he and hz are constants characteristic of a pure 
substance and are best evaluated from empirical data. Since 
Method 2 is based on the principle of corresponding states, 
the variables, pressure, density and temperature are expressed 
in the reduced form. The reduced saturated liquid density 
p, and the reduced vapor pressure p, may be expressed in 
any accurate functional or tabular form convenient for 
computing purposes. 

Equation [10], a generalized equation, quite useful for the 
representation and prediction of volumetric properties in the 
liquid region, is of a form suitable for computing by high 
speed digital computers. _ 


Application of Method 1 and Method 2 


At present the need for liquid oxygen volumetric data for 
rocket applications extends over an approximate temperature 
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density values, if 8, is considered independent of pressure. 
If one considers By = f(P), as given by Equation [6], an 
improved relationship is obtained. The procedure followed 
in this investigation was to employ Method 1 with By = 
f(P) in the temperature range 110 to 160 R and to use the 
resulting pressure-density-temperature values to evaluate 
the constants of Method 2. 3 
| 

Experimental Data on Which Table 1 Is Based Sy, P 

The data presented in Table 1 are based on graphically 
smoothed density values stemming from (1 and 10), on 
the critical constants recommended in (8), vapor pressure 
data presented in (8), enthalpy of the saturated liquid as 
given in (11) and heat capacity data presented in (4). The 
temperature dependence of sonic velocity originates from 
(9 and 12). 
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Table 1 Pressure-density-temperature relationship of liquid O- 
Liquid O, density, lb/ft? 
Pressure, 
Temp,R—~ 110 120 130 140 150 160 170 180 190 200 210 psia 
80.14 78.53 76.98 75.33 73.63 71.87 70.05 68.09 66.05 63.93 61.81 saturation 
80.16 78.56 77.01 75.36 73.66 71.91 70.09 68.11 are gig ih 50 
80.19 78.59 77.04 75.40 73.70 71.96 70.14 68.18 66.12... pie 100 
80.24 78.64 77.10 75.47 73.79 72.05 70.26 68.31 66.29 64.16 62.00 200 
80.28 78.70 77.16 75.54 73.87 72.15 70.37 68.44 66.45 64.36 62.23 300 
80.33 78.75 77.23 75.61 73.95 72.24 70.48 68.57 66.60 64.54 62.46 400 
80.38 78.81 77.29 75.68 74.03 72.33 70.58 68.70 66.75 64.73 62.67 500 
80.43 78.86 77.35 75.75 74.11 72.43 70.69 68.83 66.90 64.90 62.88 600 
80.48 78.92 77.41 75.81 74.19 72.52 70.79 68.95 67.04 65.07 63.07 7 
80.52 78.97 77.47 75.88 74.27 72.61 70.90 69.07 67.18 65.23 63.26 800 
80.57 79.02 77.53 75.95 74.34 72.69 71.00 69.19 67.32 65.39 63.45 900 
80.61 79.07 77.59 76.02 74.42 72.78 71.10 69.30 67.46 65.55 63.62 1000 
80.85 79.34 77.88 76.35 74.79 73.20 71.58 69.85 68.09 66.27 64.43 1500 
81.08 79.59 78.16 76.67 75.15 73.61 72.03 70.36 68.66 66.91 65.15 2000 
81.30 79.84 78.44 76.98 75.50 73.99 72.46 70.84 69.19 67.50 65.79 2500 
range of 160 to 210 R, and over a pressure range of saturation Nomenclature 7 
conditions to 1000 psia. Method 1 requires sonic velocity, : 
experimental values of which are available in the range of D/aT 
only 108.8 to 162.8 R (9). In order to arrive at the desired 
pressure-density-temperature relationship one of two ap- - = coefficient of thermal expansion 
proaches may be followed. It is possible to extrapolate sonic Br = coefficient of isothermal compressibility 
velocity to the desired temperature and employ the equations dP /dT = slope of vapor pressure vs. temperature 
of Method 1, or use only the experimental points of sonic H = enthalpy 
velocity, arrive at the pressure-density-temperature relation- dH /aT ha slope of enthalpy vs. temperature at saturation 
ship by Method 1, and utilize these results to evaluate the : “a conditions 
constants of Method 2, by means of which the relationship ef ‘4 oi specific a at constant pressure 
for the entire range of interest can be evaluated. 
Applying experimental sonic velocity data in Method 1, —_ = D/Dz, reduced density — 
the relationship can be computed considering Br # f(P) t “(ae /Tc, reduced temperature a 
or By = f(P). The assumption By ¥ f(P) is fairly good at . = dp,/dt, slope of reduced vapor pressure vs. reduced 
low temperatures and in a wide pressure range above 7000 : temperature at the critical point 
psia (2), but it does not hold for low pressures and high a . reduced density of saturated liquid 
temperatures. Saturation 6,7 increases rapidly with tempera- Po = reduced vapor pressure . 
ture and is infinite at the critical point. Thus at tempera- Pc, De, Te = critical pressure, density and, respectively, tem- 
tures near 260 R, in the case of liquid oxygen, at a pressure perature 
of 1000 psia, (0D/OT)p becomes positive, giving unrealistic References 
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Analog Network to Convert Surface 
Temperature to Heat Flux’ 


G. T. SKINNER? 


Cornell Aeronautical Laboratory, Inc., Buffalo, N. Y. AL 


A simple analog network has been devised, for use with 
thin film resistance thermometers, which carries out the 
mathematical operation of transforming surface tempera- 
ture into heat flux. The method is applicable to shock tube 
and shock tunnel work. A calibration unit completes the 
system which eliminates the necessity for complicated 
computational procedures. 


N A recent review article by Hall and Hertzberg* surface 

temperature thermometry was discussed and an outline 
was given of the data reduction procedure. With sometimes 
as many as 24 thin film resistance thermometers on a model 
in a hypersonic shock tunnel, data reduction becomes a major 
problem, and has raised the question of whether a suitable 
analog network could be devised to convert the gage output 
voltage, which is proportional to surface temperature, into a 
voltage proportional to heat flux per unit area. The con- 
version is usually carried out by some sort of numerical pro- 
gram which solves the one-dimensional heat flux equation 


|Kep F(+0) F(t — 7) 


Q(t) = heat flux per unit area i 

F(t) surface temperture 

surface temperature derivative 

K, c, p = thermal conductivity, specific heat and density 
of the material on which the gage is plated 


where 


Using lower case letters to denote Laplace transforms of Q 
and F, Equation [1] becomes 


qs) = VKep Vs f(s) [2] 


The analog network required to solve the problem is simply a 
filter whose steady-state response (by substituting jw for s) is 
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Fig. 1 Five-section network (for 10-millisec maximum run) 


proportional to Viw. That is, the amplitude response is 
proportional to Vw and the filter has a constant phase shift 
of 45 deg, or in other words a 3-db/octave filter. Over any 
finite range of frequencies, the Vs response can be ap- 
proximated by a transfer function - 


_ (8 + ai(s + a2)... (8 + Gn) 
(s + bi)(S + be)... (8 + 


g(8) [3] 


where the constants ai, b;, de, be, ..., b, are separated by octave 
intervals. 

In practice, a better approximation results if a small over- 
shoot is allowed at each end, and, if there are more than three 
pairs of terms, this is done by extending a; and 6, one half an 
octave. By experiment it has been determined that a five- 
section filter is most satisfactory. With fewer sections ac- 
curacy suffers, and with more sections the signal-to-noise ratio 
becomes poor. The transfer function used is 


(s + 0.692)(s + 4)(s + 16)(s + 64)(s + 256) 


(s + 2)(s + 8)(s + 32)(s + 128)(s + 739) 4] 


g(s) = 


When the resistors and capacitors have been determined, it is 
a simple matter to adjust these to obtain suitable input and 
output impedance levels, and to shift the response to the de- 
sired position on the frequency axis. The network shown 
in Fig. 1 has been designed for use on a hypersonic shock tun- 
nel where the maximum anticipated duration of a run is 10 
millisec. The starting transients affect the results; therefore, 
the optimum positioning in frequency has been determined 
experimentally. Amplitude response is within 0.1 db from 
200 cps to 2 ke per sec, and within 1 db from 2.8 eps to 12 ke 
per sec, while the phase shift is within 1 deg from 22 eps to 1.6 
ke per sec. Those familiar with hot wire anemometry may 
recognize the similarity to a very common compensator for 
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Fig. 2 Typical pressure, temperature and heat flux record from the hypersonic shock tunnel 
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thermal lag, where only a single section is used in the ap- 
proximation. 

One of these filters.is used as a feedback network to carry 
out the inverse operation on a step (of amplitude propor- 
tional to the gage current) and provide a calibration signal 
which, when passed through the amplifier containing the 
analog network, gives rise to a step output. A tuning fork 
oscillator provides an accurate spot check on the calibrator. 

The complete system, which has now been thoroughly 
tested, will be described in detail in a forthcoming Cornell 
Aeronautical Laboratory report.‘ A typical record of a hyper- 

4 This report is now published: Skinner, G. T., “Analog 
Temperature to Heat Flux,’’ Cornell Aeronautical Laboratory 
Rep. no. CAL-100, Feb. 1960. 


sonic shock tunnel run is shown in Fig. 2, where tunnel stagn:- 
tion pressure, model surface temperature and analog output 
are reproduced. The stagnation pressure shows clearly the 
steady portion of the run. The surface temperature record 
shows a rise approximately proportional to (time)!/? during 
this steady part. The analog network, operating on this tem- 
perature function, shows the constant rate of heat transfer 
during this steady part. Fig. 3 shows the comparison be- 
tween the heat flux as computed from the temperature record 
of Fig. 2 by the standard IBM 704 program currently in use, 
and the heat flux obtained by reduction of the heat flux 
record of Fig. 2 using the reference amplitude provided by the 

‘alibrator unit Discrepancies are well within the reading 
errors. 


Elliptical Orbit Characteristics 
G. S. GEDEON!' and B. J. GREER? 
Chance Vought Aircraft, Inc., Dallas, Texas 


To expedite satellite feasibility studies, the well-known 
orbital equations are presented in the form of convenient 
charts. These charts permit the establishment of the per- 
manent characteristics of the orbit and also allow follow- 
ing the change of the variable elements along the orbit. 


WO WORKING charts? are presented which permit the 
establishment of the orbital characteristics of a satellite 
after initial launch or after transfer from another orbit. The 
curves presented in these charts are based on two-body equa- 
tions and cover a range limited by the apogee radius equal to 
not more than twice the perigee radius. 
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where epi 


The charts —_ six basic parameters of an orbit which 
are defined on Fig. la. Knowing any two of these parameters, 
the remaining four may be read from the charts, and the 
effects of changing one of them may be readily determined. 
Thus, the charts serve as a simple and rather accurate way to 
minimize the time-consuming calculations connected with 
orbits. 

To maintain the accuracy of four significant digits, which 
can be read from the charts, it is assumed that a desk calcula- 
tor will be used to convert to and from the unitless quantities. 

The equations defining the relationships among the variables 
which are plotted in Fig. 1a are as follows: : 

Radius ratio 


r 
l-—e 


(7) [2 cos? @ = eccentricity 
VOW. v,) 


See (1).4 


* Numbers in parentheses indicate References at end of paper. 
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Fig. la was obtained by assuming a set of (V/V,) values for 
a constant ¢ and calculating the corresponding (r/R) and zx 
values. Now plotting (r/R) against xz, a constant ¢ curve 
is obtained, on which the defining points are labeled ac- 
cording to the assumed (V/V,) values. Having drawn in 
this manner the constant ¢ curves, the points of equal (V/V,) 
can also be connected to yield the constant (V/V,) lines. 

Fig. la serves to establish the orbit from the burnout con- 
ditions (launch or transfer), or to obtain the burnout condi- 
tions for achieving a required orbit. Note that the line seg- 
ment corresponding to an extremicity value represents the 
orbit. Any point on this line segment yields the radius, the 
velocity and the angle of elevation of the orbit if the perigee 
radius is already known. 

The plot of Fig. 1b was made by assuming constant (r/R,) 
values and by varying e, calculating @ from the ais 
equation for the true anomaly 

1) 


= | - 
Then 6 was plotted against x (converted from e) and the curve 
was labeled according to the assumed (r/R,) value. The 
second set of curves plotted on Fig. 1b was obtained by as- 


suming constant ¢/P values, and by again varying e, calculat- 
ing E from the following formula 


Calculate 


Rp, = 5000/1.07 = 4673 
Ra = 4673 X1.5 = 
Tmean = (4673 + 7010)/2 = 


For this value of rmean, read P = 2.52 hr. = 


Using Fig. 1b and (r/R,)) = 1.07 and x = 1.5, read @ = 
52.5 deg and (t/P)) = 0.10. 

Since the angle of elevation was positive, the (zero) perigee 
passage of the orbit is 4 = 0.10 X 2.52 = 0.252 behind the 
burnout point. 

It is assumed that a retrorocket which produces a velocity 
decrement of 900 fps will be fired 1.26 hr after burnout. Find 


the characteristics of the new orbit. eal 


First calculate 
0.252 1.26 
= = 0.60 after the zero perigee 


2.52 


—0.40 before the first perigee passage 


For this value and x = 1.5, one can read from Fig. 1) 
6, = —155 deg and (r/R,); = 1.465 and from Fig. la ¢ = 
—6 deg and (V/V,); = 0.91. 


t t t 
ee i. P P 5! 24 P 3* sin "> sin 7 > ... = eccentric anomaly 


See (2). 
From E and e, 9 was calculated by the following equation 
for the true anomaly 


= cos—[(cos — e)/(1 — e cos E)] 


Finally 6 was plotted against x (converted from e) and labeled 
according to the assumed ¢/P value. 

Fig. 1b is convenient for finding the time necessary to ar- 
rive at a certain point on the ellipse which is specified by a 
true anomaly, or vice versa. The initial anomaly can be 
found by transferring from Fig. la to Fig. 1b with the same 
radius ratio and extremicity values. 

In order to eliminate the necessity of calculating the circular 
orbit velocity and the period of an orbit, Fig. la is equipped 
with two additional curves which are based on the following 
equations 


= RoV qo/r 


wheres 
go = 32.146 ft per sec? bile : 
Ro = 20,926,083 ft, 
V, = taken for (r = rmean) 


Both figures represent only one half of the orbit. Due to 
symmetry, the second half can be handled by reversing the 
direction of motion. This will produce negative values for the 
anomaly, angle of elevation, and orbiting time. 


Sample Problem 


Determine the orbital characteristics of a launched satellite 
which had the following burnout conditions: 


= 5000 miles Vo = 24,706 fps ni i 
Solution: 
Using Fig. la and ro = 5000 miles, read (V,)o = 23,090. 


Calculate: (V/V.)o = 24,706/23,090 = 1.07. 


go = 8 deg 


1.07 


Thus the radius at the control point is 
r; = 1.465 X 4673 = 6846 
The corresponding circular velocity is 
(V.)1 = 19,730 
and the actual velocity is 
V; = 0.91 X 19.730 = 17,954 
The retrorocket decreases the velocity to 


V,’ = 17,054 
which yields 
V\’ 17,054 
(7) 19,730 = 0.8644 { 
For this (V/V,),’ and ¢; = —6 deg the following values can 
be read os 
= 1.745 (r/Rph’ = 1.72 
0,’ = —164 deg (t/P);’ = —0.425 
The new perigee radius is then 
R,' = 6.846/1.72 = 3980 
The mean radius of the new orbit is 2 


For this mean radius the new period is P’ = 2.28. 


Remarks 

Since the submission of this paper, a new set of detailed 
charts have been worked out to extend the eccentricity range 
up to 5. These are also available on request from Chance 
Vought Research Center, Dallas 22, Texas. 
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Comparison of Transpiration and _ 
Ablation Cooling 


CURTIS C. BEUSMAN! and JOEL WEISMAN? 


Nuclear Development Corp., White Plains, N. Y. 

ECENT summaries of ablation and transpiration meth- 

ods for providing thermal protection of re-entry satellites 
have pointed out the importance of the heat blocking effect, 
commonly called the transpiration factor, in increasing the 
effective heat capacity of a given coolant material. Adams 
(1)* presented data comparing the effectiveness of Teflon, 
fiber glass and foam quartz in providing protection and in- 
sulation for a recoverable satellite exposed to a rather mod- 
erate heat load during re-entry. Scala (2) on the other hand, 
presented ablation and transpiration data applicable to the 
higher heat transfer rates occurring during a ballistic missile 
re-entry. In relation to this study, Scala considered liquid 
air as a transpiration coolant and a quartz-like refractory 
oxide-as the ablative shield. For all values of the ballistics 
factor W/C,)A in the range 50 to 200, Scala concluded that 
ablation is more effective than transpiration cooling on a 
coolant mass consumption basis. 

We have investigated these conclusions, and extended the 
calculations to include other more promising transpiration 
coolants, i.e., water and lithium, to illustrate the added sig- 
nificance of the molecular weight of the transpiration coolant 
as well as the coolant heat of vaporization. Although the 
use of liquid air as a transpiration coolant does provide for 
important analytical simplification in the boundary layer 
analysis, the added engineering complexity and system 
weight of a eryogenic coolant storage system seems un- 
warranted. In view of the experimentally verified effective- 
ness of water (3) as a transpiration coolant, its application in 
this instance could prove fruitful. 


Ballistic Re-Entry Comparisons 


Scala expresses the mass transfer rate of the coolant m as 


Qaero + Qrad 
— hres) + AQ/Am = 


m= 


where 
a2 
Qaero = aerodynamic heat transfer rate 
Qrad = radiative heat transfer rate 
hw = enthalpy of the coolant at wall conditions 
Tiree = enthalpy of the coolant at reservoir conditions 


(AQ/Am) = “effectiveness quotient” representing the heat 
blocking effect, caused by boundary layer 
thickening 


Using the values of m, Qsero, Qraa and the wall and reservoir 
temperatures given by Scala for liquid air transpiration cool- 
ing, the values of (AQ/Am) used may be computed. The 
computed values of the heat blocking effect (up to 17,000 | 
Btu per lb) are appreciably greater than the enthalpy rise 
(hw — hres) for either liquid air cooling (hw — hres = 1100 
Btu per lb) or water with a reservoir temperature of 100 F 


(hw — hres = 2500 Btu per Ib). This fact would indicate that 
the choice of coolant is relatively unimportant, a fact ob-— Liquid air — transpiration cooling 
served early in the study of transpiration cooling (3). How- 2.0 ee) 
ever, this ignores the enhancement of the “effectiveness & -_ | ee Water — transpiration cooling 
quotient” through a reduction in the molecular weight of the 4 wat os 
coolant gas. Adams (1) points out that the blocking effect d 1.0 4 apt) | eee Vaporization cooling I= 0.25 
is inversely proportional to the coolant molecular weight 7 | Lithium — transpiration cooling 
to about the one third power. Sublimation, = 1 

Received Jan. 8, 1960. 0 100 180 200 
Fig. 1 Coolant mass required stagnation point of re-entry 

satellite 


’ Numbers in parentheses indicate References at end of paper. — 


The effect of the increased enthalpy rise and lower molecu- 
lar weight for water cooling is shown in Fig. 1. It is seen 
that the coolant requirements are reduced nearly 30 per cent 
below the requirements of liquid air cooling. This system 
is thus more nearly competitive with the low gasification ratio 
([) materials, e.g., quartz discussed by Scala. (The gasi- 
fication ratio represents the fraction sublimed of the total 
material ablating from the re-entry body.) 

As a further extension of this approach, lithium has also 
been studied as a transpiration coolant. Lithium is attrac- 
tive because its lower molecular weight further reduces heat 
transfer to the re-entering body and its very high enthalpy 
of vaporization provides an additional energy sink. The 
required mass transfer rate for lithium cooling (based on a 
reservoir temperature of 400 F) is also shown in Fig. 1. It is 
seen that the coolant requirements are less than one third 
those needed when liquid air is used. Lithium weight re- 
quirements are even less than those of low gasification ratio 
ablative materials. 

Both water and lithium can undergo chemical changes in the 
dissociated boundary layer; i.e., water can dissociate and 
lithium can react with dissociated oxygen at the wall. These 
additional effects have not been analyzed. In the case of 
water, the additional heat sink represented by dissociation 
(~20 per cent at 3500 R) can improve the performance of 
water in this instance (4). Lithium combustion could be 
included; if the blowing rate is high enough, no additional heat 
load is anticipated (5). In addition, ionization of the lithium 
can be significant, both from a cooling and from a telemetry 
standpoint. 


Recoverable Satellite Re-Entry 


We can also estimate the heat shield weight for satellite 
recovery using water as a transpiration coolant. For the 
case presented by Adams, severe weight penalties were in- 
curred in ablative shields of quartz or fiber glass because of 
the insulation required to keep inner temperatures sufficiently 
low to avoid damage to the payload. In transpiration cool- 
ing, the countercurrent coolant flow prevents heat from pene- 
trating through the porous skin (6) except for the outer- 
most portion of the re-entry surface. 

For this example (a satellite entering the atmosphere at 2 
deg from the horizontal and subjected to a total heating of 
7650 Btu per ft?) with an assumed outer wall temperature of 
800 F, the required weight of water coolant is 1.1 lb per ft?. 
This weight includes insulation in the sense that transpiration 
cooled surfaces do not require added interior thermal pro- 
tection. In Table 1 the weight requirement noted is com- 
pared with the ablative insulation schemes discussed by 
Adams. On this basis, the transpiration cooling scheme is 
the most effective of those considered. 

This estimate neglects the weight of the porous surface, as 
well as the supply reservoir and distribution system. Simple 
lightweight distribution schemes have been designed (7) to 
provide structural strength for a thin porous sheath, as well 


573 


=: 
see 
ne 
hy 
ty 
ty 
nil ‘3 
ks 
= 
ly 
1 
1 
3 

BES 


Table 1 Weight comparison of thermal protection schemes 
Total weight, lb/ft? 


Material Tsurtace, F (ablation and insulation) 
fiber glass ~3600 15.5 
foam quartz 3600 ra 
Tefion 800 
water (transpiration 
cooling) 800 | 


as adequate flow control for various re-entry trajectories. 
Other schemes, using capillary action to provide coolant dis- 
tribution or self-regulation of coolant injection via gels, can 
be presented. The present direct comparison on the basis 
of total mass jettisoned during re-entry is believed to be fair 


since both the direct ablative as well as the transpiration 
methods require added material to insure structural integrity 
at the conclusion of re-entry. 


References 


1 Adams, M. C., ‘‘Recent Advances in Ablation,”” ARS Journat, 
29, no. 9, Sept. 1959, pp. 625-632. 

2 Scala, S. M., “Thermal Protection of a Re-Entry Satellite,’"’ ARS 
Journa., vol. 29, no. 9, Sept. 1959, pp. 670-672. 

3 Jakob, M. and Fieldhouse, I. B., ‘‘Cooling by Forcing a Fluid throug! 
a Porous Plate in Contact with a Hot Gas Stream,’ Proc. 2nd Symposiur 
on Heat Transfer and Fluid Mechanics, ASME, 1949, p. 191. 

4 Sanger-Bredt, I., ‘‘The Characteristics of Hydrogen and Water as 
Working Gases for Reactor-Heated Rocket Motors,’’ Astronautica Acta, vol. 
3, no. 4, 1957. 

5 Sutton, G. W., “Adiabatic Wall Temperature Due to Mass Transfer 
Cooling With a Combustible Gas,”” ARS JourNaAt, vol. 29, no. 2, Feb. 195, 
pp. 136-137. 

6 Grootenhuis, P., ‘‘The Mechanism and Application of Effusion Coo!- 
ing,” J. Royal Aero. Soc., vol. 63, pp. 73-89. 

7 Beusman, C. C., “Feasibility Study of Transpiration Cooling for the 
Mark II Polaris Re-Entry Body,’’ Lockheed Missile and Space Division, 


Combustion Instability in Solid Pro- 


pellant Rocket Motors 
fw 
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JN A Technical Comment by Green and Nachbar (1),? the 
& subject of burning rate of solid propellants during high 
frequency oscillatory combustion was discussed further, and 
objections were raised to a comparison [by Price and Sofferis 
in (2) ] of experimental work with earlier statements [by Green 
in (3)].. The purpose of the present note is to cite new experi- 
mental data on this problem and to clarify, if possible, the 


sources of controversy in (1 and 2). A more detailed version 


of the present note is presented in (4). peatie 4 


New Experimental Results 


The new experimental results were obtained in a continua- 
tion of the program described in (2), and the details of the 
experimental procedure will not be repeated here. Tests 
have now been made on a wide range of double-base and 
composite propellants, totaling about 800 firings on 15 differ- 
ent propellants. Of the results obtained with conventional 
propellants,* the following generalizations relevant to the 
Green-Nachbar theory can be made: 

1 Burning rate of propellants always increases during os- 
cillatory burning, as indicated by the accompanying increases 
in pressure (Fig. 1). 

2 Partially burned charges show that the greatest in- 
crease in burning rate is in the region where the amplitude 
of the velocity oscillations of the gas is a maximum (Fig. 2). 

3 The observed increase in burning rate during oscillatory 
burning exhibits no dramatic dependence on composition of 
propellant or on frequency of the oscillatory behavior (in- 
vestigated for one propellant over a frequency range of 
1500 to 6500 cps). It has further been observed [the re- 
sults reported in (2) and cited by Green and Nachbar in (1) ] 
that consistently different behavior is exhibited by the family 
of mesa-burning double-base propellants.‘ 
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of burning rate on pressure. 

‘ Propellants characterized by negative pressure exponents 
of the burning rate over a limited range of pressure. Three 


mesa propellants were tested. 
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4 With these propellants, the burning rate (averaged over 
the charge surface) decreases during oscillatory burning [(5); 
also Fig. 3]. 

5 The maximum decrease in burning rate occurs in the 
region where the amplitude of the tangential velocity oscilla- 
tion of the gas is a minimum. 


Comparison With Green-Nachbar Theory and Earlier 
Notes 


With respect to the items listed, the following comments 
are made: 

1 A wide variety of propellants was tested, with varying 
degrees of sensitivity of burning rate to oscillatory behavior. 
However, all conventional propellants exhibited appreciable 
increases in burning rate. The observation of Green and 
Nachbar (1,6), that instances must be expected of no change 
or of decreases in burning rate, was not borne out in the 
present work. In the face of difficulties in performing quan- 
titative measurements of pressure oscillations in rocket 
motors, and in the face of the possibility of other mechanisms 
producing pressure peaks, it is doubtful that severe rate 
changes and pressure peaks are ever caused directly by low 
amplitude gas oscillations. 

2 The occurrence of enhanced burning rates predomin- 
antly in regions of severe velocity oscillations is suggestive of 


Pressure psi 


Pressure psi 


Fig. 1 Pressure-time record for firing with severe oscillatory 
combustion, showing variation in pressure (lower trace) with 
amplitude of pressure oscillations (upper trace). Conventional 
(r = cp") propellant. The pressure would normally increase 
gradually during burning even without instability, because of the 
increase in burning surface of the internal burning charge. The 
irregular increases in pressure are due to combustion instability 
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Fig.2 Variation in burning rate with position along the charge 4 as 
determined from interrupted burning tests. The ordinate is dis- 
tance burned into the charge web minus distance burned in at 
the end of the charge (velocity node, stagnation point), divided by 
the distance burned in at the end. The abscissa is distance from 


the end of the charge. JPN (conventional) propellant 


an erosion-like mechanism, and does not appear to be con- 
sistent with the Green-Nachbar theory of (1 and7). Further 
scrutiny of the theory on this score (4) leads the writer to the 
conclusion that, contrary to (1), the theory provides no basis 
for an unequivocal statement as to where the largest increases 
in burning rate should occur on the propellant charge. In- 
stead, the location where the burning rate is changed most by 
oscillatory behavior depends upon design, fluid dynamic and 
propellant variables, and even changes during a single test 
(according to the theory). This is hard to reconcile with the 
experimental observation in partial burning tests that the 
greatest amount burned from the charge web was always at 
the velocity antinode, no matter what conventional pro- 
pellant was used or when burning was interrupted during the 
test. 

3 In the series of tests with different charge lengths, a 
rather large range of oscillation frequencies was represented, 
with all other variables being the same at some times during 
each test. No singular dependence of burning rate on fre- 
quency of oscillation was evident. This result also appears 
to be inconsistent with the Green-Nachbar theory, and sug- 
gests the relevance of a mechanism similar to erosive burning. 

4 The results obtained with mesa propellants are probably 
not relevant to the Green-Nachbar theory, because of the 
complicated solid-phase decomposition process and because 
of the characteristic severe rippling of the burning surface 
during oscillatory combustion (Fig. 4). Thus, the decreased 
burning rate exhibited during oscillatory combustion is 
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Pressure psi 
Pressure psi 


_ Fig. 3 Pressure-time record for firing with severe oscillatory 


combustion, showing evidence of decreased burning rate of mesa 

propellant during oscillatory combustion (pressure rises when 

severity of oscillations decreases). In the absence of in- 

stability, this test would show a smoothly progressive pressure- 

time curve corresponding to an internal burning circular cylindri- 
cal charge 


Fig. 4 Surface ripples typical of a mesa propellant recovered by 

interrupted burning during oscillatory combustion. Gas flow 

and gas oscillation were oriented in a horizontal direction in the 

figure. Distance between crests of ripples is about 0.06 in. 

(wave length of the gas oscillation was about 24 in., displacement 
of the gas during a cycle of oscillation was about 4 in.) 


probably not a valid demonstration of the decreased burn- 
ing rate predicted by Green and Nachbar (1). The 
example of reduced burning rate reported in (2) as an 
exception to the theory, and subsequently interpreted 
by Green and Nachbar in support of their theory (1), 
was a result observed with a mesa propellant. Thus, 
the results of (2) are probably an exception to the theory 
as originally reported, but do not necessarily represent 
an irreconcilable contradiction to the theory, provided that 
the theory is suitably restricted. 

5 In view of the remarks 2 and 4, the arguments by Green 
and Nachbar (1) comparing theory with the results of (2) 
seem to be invalid with respect to the site of maximum burn- 
ing rate change as well as to the | oceurrence of decreases in 
burning rate. ws 


Conclusions 


In summary, the experimental results of the present pro- 
gram suggest that the burning rate increases during oscilla- 
tory combustion reflect a mechanism much like the more 
familiar erosive burning phenomenon of steady-state burning. 
The present results do not preclude the possibility of other 
mechanisms. However, they do seem to be in contrast to 
any theory that predicts critical sensitivity of burning rate 
to gas oscillation frequency or other propellant and environ- 
mental variables. It is suggested that the critical conditions 
predicted by the Green-Nachbar theory of resonance burning 
may be a consequence of the unsophisticated model of burn- 
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ing rate response to temperature and of the linear model of 
heat transfer from flowing gas to solid propellant. As such, 
the critical conditions predicted by the theory might not be 
relevant to severe instability and could not be identified 
unless quantitative tests with low amplitude instability were 
made. 
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XAMINATION of the new experimental data to which 
Price refers in the preceding Technical Note shows that 
the source of the difficulties, which have been encountered in 
applying the model of his Reference (7)* to the inter- 
pretation of these experiments, is in the linearization which 
was performed on the equations. The analysis of (7), p. 519 
specifically assumed small amplitude sinusoidal perturba- 
tions from steady-state values. Since only terms linear in 
AT, and AF are retained in the expansions, it then follows 
that there is no predicted change in mean burning rate with 
these small perturbations [(7), p. 520: This is easily seen by 
averaging Eq. (5f) over one cycle]. The only way in which 
finite amplitude changes in surface temperature and finite 
changes in mean burning rate could arise from this linearized 
analysis is if operating conditions are at or very close to the 
predicted critical values. Under these conditions, even small 
(infinitesimal) amplitude oscillations of chamber pressure 
could give rise to finite changes in mean burning rate. Be- 
cause of the limitations of the linearization, however, the 
amplitude of the mean burning rate at critical frequency could 
not be determined. 

The conditions of oscillatory burning under which Price 
has reported sizable mean burning rate changes are not those 
which were assumed in the analysis of (7). Price’s 
experimental results, reported in detail elsewhere, have 
shown that sizable mean burning rate changes occur in his 
self-excited burner only with measured peak-to-peak ampli- 
tudes of pressure oscillation of the order of the mean chamber 
pressure. 
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The predicted behavior of the linearized model is similar 
to that which had been observed in conventional solid 
rocket motors with internal burning tubular charges by ex- 
perimenters before Price (6). Price’s experimental results, 
at least those which have so far been reported, are in apparent 
disagreement with this earlier evidence, and, indeed, by 
implication they contradict the fundamental premise on which 
the linearized model is based. The occurrence and the 
amplitude of abnormal burning rates are correlated in every 
instance with the occurrence and the amplitude of chamber 
pressure oscillations (sonance). No critical or ‘‘resonance”’ 
conditions of burning are found. This experimental con- 
clusion is corroborated by the recently reported work of 
Brownlee.* 

The experiments of Price and Brownlee have all apparently 
been performed with great care and with superior instrumenta- 
tion. However, there still appears to be some conflict between 
these results (see Price’s generalization and comment 3) and 
what has been reported elsewhere as solid rocket motor experi- 
ence. It is not appropriate here to enter further into discus- 
sion of this, except to point out the possibility that these experi- 
ments may involve an artificially created situation not pres- 
ent in actual motor designs, a situation in which the acoustic 
damping has been so much reduced that high frequency os- 
cillations of finite amplitude readily occur. 

Price’s comments 1 and 4 are in no way surprising or in 
contradiction to the theory, for Price has not given any 
data for his propellants from which the algebraic sign of 
[(E/2RT..) — 1] can be precisely determined. Therefore, 
no theoretical prediction can be made relevant to his tests. 
However, rough estimates of typical values’ (EZ = 20 kcal per 
mole, 7. = 1200 K) will clearly indicate that, for conven- 
tional propellants, this quantity will be positive, and thus 
conventional propellants will be expected to undergo an in- 


4 Brownlee, W. C., “An Experimental Investigation of Un- 
stable Combustion in Solid Propellant Rocket Motors,” Doctoral 
Thesis, California Institute of Technology, Pasadena, Calif., 
1959. 

5 Summerfield, M., Sutherland, G. 8., Webb, M. J., Taback, 
H. J. and Hall, K. P., “Burning Mechanism of Ammonium Per- 
chlorate Propellants,’’? ARS preprint 737-58, Nov. 1958. 
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crease in mean burning rate in oscillatory burning [see 
(7), p. 524]. This behavior may be different for mesa- 
type propellants. At least the theory allows for the pos- 
sibility of this difference, although we have no data on 
which to base a numerical prediction. 

As to comments 2 and 5, the controversy which has arisen 
here is the responsibility of the present authors, who 
in (1 and 7) drew one conclusion which overextended 


the applicability of the linearized theory. No prediction 
can be made from this linearized theory regarding the location 
of the site of maximum burning rate change with respect 
to the standing wave of velocity perturbations. Conse- 
quently, the conclusion we reached, that the maximum burn- 
ing rate change necessarily occurs at the nodes of the velocity 
oscillation, is admittedly wrong, and we are indebted to Price 
for questioning this point. 


Visual Sensitivity of Residents of Other 
Planets 


ALLEN E. FUHS'! 
Space Technology Laboratories, Inc., Los Angeles, Calif. 


i ECENTLY it was announced (1)? that a search would be 

made for radio signals from other planets. Naturally 
this implies the existence of intelligent beings on other 
planets. One might speculate as to the nature of the eyes 
of these beings. Basic for the speculation are several facts 
about the eyes of human beings as they have evolved on 
Earth. The human eye as evolved on Earth is the result of 
the interplay of two factors which are illustrated in Fig. 1. 
There are several windows in Earth’s atmosphere, one of 
which coincides with the peak of the spectral distribution 
curve of radiation from the sun. The eye is sensitive to 
radiation that is at the peak of solar radiation and at this 
window of the atmosphere. 

Consider the planet of an M star in the Harvard spectral 
sequence (2). M stars have temperatures in the range 2500- 
3200 K, so that the peak radiation will occur in the range 
0.94-1.20 uw. If the planet in the M star system has an 
atmosphere similar to Earth, i.e., similar concentrations 
of CO,., ete., then the inhabitants will probably 
have eyes which have maximum sensitivity at 1 wu. There 
is a “window” near 1 y, as can be seen in Fig. 2 which is based 
on the results of Gebbie et al. (3). Both Earthmen and 
residents of the M star planet could see; however each 
would observe the other at different wave lengths. This 
is illustrated in Fig. 3. 

Next consider the planet of a late O star, when peak radi- 
ation occurs at 0.3 wu. Assuming that the atmosphere of this 
planet is similar to Earth, then its people should have similar 
eyes except for greater sensitivity for violet. 

It is conceivable that the atmosphere of a planet in some 
other “stellar’’ system differs from Earth’s atmosphere; in 
fact, it probably does differ. If there were strong absorption 
in wave length range from 0.4 to 0.75 y, the beings of that 
planet would probably have sensitivity to radiation in one of 
that planet’s windows. Hence an Earthman would be blind 
on that planet whereas the local residents could see. This 
is the condition prevailing on some of the other planets of 
the stellar system; however there probably are not any 
“local residents.” 
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Fig. 1 Maximum sensitivity of human eye occurs at the peak of 
solar radiation 
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Fig. 3 Maximum radiation occurs at a longer wave length. 
There is a window at the same wave length 
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Book Reviews__ 


Feedback Control Systems, by J. C. Gille, 


M. J. Pelegrin and P. Decaulne, 
McGraw-Hill Book Co., Inc., New York, 
1959, xx + 793 pp. $16.50. 
Reviewed by G. J. Murpuy 
Northwestern University 


According to the authors, this book is 
intended both as a textbook for seniors 
and first-year graduate students and as a 
reference book for the engineer working in 
the field of automatic control. The aim 
of the authors is: “By progressing from the 
elements of the subject to provide in a sin- 
gle book a treatment both of the overall 
theory of feedback control systems and of 
their components.” 

The text is divided into five parts, titled 
in sequence: Dynamics of Linear Systems, 
Linear Servo Systems, Nonlinear Servo 
Systems, Components of Servo Systems 
and Basic Design of a Servo System. The 
major concepts are presented very clearly, 
and the entire volume is written in a 
literary style that is unusually good. 

There are numerous references to sig- 
nificant French, German, Russian, English 
and American papers. In several in- 
stances in which different methods have 
been proposed for solving a given problem, 
the different methods are compared, and 
their relative merits are discussed. 

For some reason, the authors chose to 
employ standard symbols in nonstandard 
manner (e.g., r(¢) for output, instead of 
input; e(t) for input, instead of error), and 
this may be disturbing to the engineer 
accustomed to thinking in terms of the 
standard definitions. 

Another item to which many readers 
will take exception is the attitude ex- 
pressed (on p. 126) in the last sentence of 
the following passage: ‘‘One can study the 
H(jw) function, i.e., the behavior of 
H(s) on the imaginary axis in the Laplace 
plane. This is the harmonic approach. .. . 
Conversely, one can study the singularities 
of H(s), that is, study the locations of its 
poles and zeros in the Laplace plane. This 
is the pole-zero-configuration approach. . . . 
Technically speaking, the harmonic ap- 
proach is more general in its application.” 

Similarly, there will be many who dis- 
agree with the statement ‘flow graphs in- 
volve more advanced thinking than do 
block diagrams,” which implies a non- 
existent essential difference between the 
signal flow graph and the conventional 
block diagram. 

The accent throughout the book is on 
the application of control theory, rather 
than on the theory itself. There are 
numerous practical examples of consider- 
able value. However, the details are often 
lacking. One receives the impression 
that, in spite of the large number of pages 
in the book, the treatment of many topics 
is somewhat superficial. Inexplicably, 
an extended treatment (45 pages) is given 
to first-order and systems. 
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Further evidence of the greater interest 
in applications than in theory is the scar- 
city of proofs, the weakness of many of the 
proofs that are presented, and the general 
lack of rigor in mathematical statements. 
For example, on p. 116 one reads, ‘‘When 
the transfer function of a system S has a 
single pole equal to zero, it can be written 
in the form H(s) = P(s)/sQ,(s), where 
0,(0) ¥ 0,” no mention being made of a 
similar restriction on P(0). 

Perhaps more serious is the following 
misleading statement, made on p. 166, in 
reference to Routh’s criterion: “Jf all the 
coefficients in a row are zero, this indicates 
that there are two equal and opposite 
roots, in particular, two conjugate roots 
which are purely imaginary; i.e., the 
system is just oscillatory.” 

Also serious is the failure to distinguish 
between a stationary random process and 
an ergodic random process. This failure is 
followed naturally by erroneous state- 
ments, such as: “If y(/) is a stationary 
random function, its autocorrelation func- 
tion ¢(r) is defined as the mean value of the 
product y(t) y(t + 7). This definition 
leads to the equation 


T 
oat... y(t) y (t + 7) dt 


which is independent of the sample used.” 

The value of this book is enhanced by 
the inclusion of some topics that are 
unusual in a book on automatic control. 
Among these are a treatment of transfer 
matrices, as a generalization of transfer 
function, and Tsypkin’s method for study- 
ing oscillations in relay control systems, 
which is shown to be a generalization of the 
describing-function method of Kochen- 
burger. In addition, the book is ac- 
companied by a set of very useful graphs, 
including Bode gain-and-phase curves for 
first-order systems and (with damping 
ratio as a parameter) for second-order 
systems, and a Nichols chart. Numerous 
other graphs useful in design are included 
within the text. 

The material on components is largely 
qualitative and hence of limited value. 
However, an extended treatment of the 
design of servo motors is included. 

Without doubt, the book is impressive. 
Most control engineers will find in it some- 
thing of value, and it will probably be 
used widely as a handbook. 
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Advanced Propulsion Systems, edited by 
Morton Alperin and George P. Sutton, 
Pergamon Press, New York, London, 


1959, 237 pp. $6 


Reviewed by A. K. OPPENHEIM 
University of California 


The book represents the proceedings of 
the unclassified portion of a symposium 


Studies the Air Force Office of Scientifi: 


Research in December 1957. It contains 
a good deal of information of current in- 
terest on advanced propulsion systems. 
All of the papers are concerned with the 
fundamental aspects, and the exposition is 
quite elementary, so that an uninitiated 
reader can follow the text without much: 
difficulty. The subject matter is, of 
course, quite speculative in nature and 
consequently the book is full of imagina 
tive but, as yet, untried concepts. 

Since each paper represents essential] 
an entirety of its own, the best way to 
_describe the book is by a résumé of it: 
contents. The basic features of ion pro- 
pulsion systems are reviewed by S 
Naiditch, and design aspects of the ion 
rocket engine are described by Robert H 
Boden. Some results of experimenta! 
studies on the thrust from a continuous 
plasma jet are reported by Gordon Cann, 
Adriano Ducati and Vernon Blackman: 
the general subject of thrust from plasma 
is covered by Winston H. Bostick. Krafft 
A. Ehricke presents a comprehensive com- 
parison of various advanced propulsion 
systems; Rolf D. Buhler and Peter J. 
Gingo discuss the potential aircraft ap- 
plications of closed gas cycle nuclear 
power plants; and R. A. Koehler reports 
on speculations about future trends in 
light weight electric equipment. The 
problem of direct power conversion is 
presented by John H. Huth and George 
Safonov. The electrostatic generators 
applicable to advanced propulsion systems 
are described by A. John Gale; their mag- 
netohydrodynamic counterpart is eluci- 
dated by Richard J. Rosa. The status 
quo of liquid propellants is summarized by 
John F. Tormey. The particular features 
of the ozone-fluorine system are discussed 
by Edward N. Hall, metals and metal hy- 
drides by James M. Carter, the problem of 
stabilized free radicals by Julius L. Jack- 
son, and the research on free radicals as 
rocket propellants by George Moe and 
Don L. Armstrong. The book contains, 
also, a paper by Homer J. Stewart on 
human hazards of space flight, and one on 
system considerations affecting new pro- 
pulsion applications, by Allen F. Donovan. 


Quantitative Molecular Spectroscopy and 
Gas Emissivities, by S. S. Penner. 
Addison-Wesley Publishing Co., Inc., 
Reading, Mass., 1959, xv + 587 pp. 
$15. 

Reviewed by A. G. Gaypon 
Imperial College, London 
In the period between the two World 
Wars, the emphasis in spectroscopic re- 
search was on the determination of atomic 
and molecular structures. This involved 
analyses of spectra in order to build up 
energy level schemes and their quantum 
mechanical interpretation. This phase 


in research, at least in regard to the simpler 
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atoms and molecules, is nearing comple- 
tion, and the emphasis is now increasingly 
on the quantitative use of spectroscopic 
data. The lead for this type of work has 
come from astrophysics, detailed studies of 
intensity and line contour being used to 


itor calculate the physical and chemical con- 
stitution of stellar atmospheres. Similar 
; calculations are, however, necessary in 
other fields, including the spectroscopic 
tific study of combustion processes, tempera- 
ains ture measurements of plasma, study of 
in- radiation from shock heated gases and 
ms. radiative heat transfer to space vehicles on 
the re-entry. 
his Underlying the whole subject is a treat- 
tec| ment of the profiles of spectrum lines. 
ich These profiles are the result of a combina- 
a tion of the natural (uncertainty) line 
ni width with Doppler broadening, collision 
19 broadening, self-absorption, and with in- 
strumental resolving power and slit width. 
lly Professor Penner deals thoroughly with 
6 this subject; he goes on to discuss emis- 
it sivities of hot gases; there is an especially 
0 long (92-page) chapter on infrared meas- 
S urements; the emissivity of heated air is 
7 dealt with fully, as this is relevant to the 
1. re-entry problem; the chapter on emissiv- 
al ity of hydrogen atoms is mainly of astro- 
Is physical interest. The three final chapters 
n. are of especial value for combustion and 
1: rocketry; they are on measurement of 
2 flame temperature, on radiation, heat 
rt transfer and burning rates of solid pro- 
% pellants, and on liquid-fuel rocket motors. 
n This is an advanced monograph, as- 
j suming a good knowledge of molecular 
4 spectroscopy and quantum mechanics. 
r There is very little overlap with standard 
5 books on these subjects. As Penner says, 


the selection of subject matter has been 
strongly biased by the researches of his 
: own group, which have, of course, been 
; most valuable. Thus, the infrared work 
is well treated, both theoretically and 
experimentally. 

The discussion of electronic transitions 
is much less complete; there is only a 
casual reference to Hund’s coupling cases, 
and rotational intensity factors are not 
given explicitly. 

The book is not limited entirely to 
molecular spectroscopy; the other sub- 
jects include photo-detachment of elec- 
trons, the continuous spectra produced by 
ion recombination and by free-free transi- 
tions, and optical dispersion. 

Penner tends to gloss over the under- 
lying physical principles and to concen- 
trate on the mathematical formulation of 
the problems. This may add to the pre- 


cision but certainly not to the clarity of the 
treatment, and the reviewer, whose back- 
ground contains mostly experimental work, 
has found the book very heavy reading and 
has failed to understand many sections of 
it. A major weakness is the inadequate 
subject index; with a more detailed index, 
the book would have been much more use- 
ful for reference work. 

There is a recent tendency to publish 
research papers, especially those with 
aeronautical applications, in obscure re- 
ports of research organizations and uni- 
versities, rather than in established scien- 
tific journals. These reports are usually 
not abstracted and are not available in 
small libraries. Penner has apparently 
had access to a large number of these 
reports and has made good use of them, 
the book making readily available for the 
first time many new results and detailed 
calculations. There are few references to 
Russian work, though; and in some sec- 
tions, e.g., the calculation of vibrational 
overlap integrals, the bibliography is 
incomplete. 

From what has been said, it will be clear 
that the book is at an advanced level, and 
the statement on the dust cover that it is 
written on the advanced undergraduate 
level is not true. However, for the special- 
ist in the field, with adequate spectro- 
scopic and mathematical background, the 
book will be valuable. 


Absorption Spectrophotometry, by G. F. 
Lothian, Macmillan Co., New York, 
1958, viii + 246 pp. $10.75. 

Reviewed by A. E. Fuus 
Space Technology Laboratories, Inc. 


In the last two decades absorption 
spectrophotometry has increased in its 
usefulness as a research tool. New instru- 
ments have been developed, numerous 
significant papers have appeared in the 
literature, and new applications have 
evolved. In “Absorption Spectropho- 
tometry,”’ which is a second edition of a 
previous book ‘‘The Practice of Absorp- 
tion Spectrophotometry,’’ these topics are 
concisely discussed. The author has 
divided the book into three parts. 

Part I deals with the laws and principles 
of the subject. Integrated band inten- 
sities, pressure broadening, effects of 
finite slit-width, scattered radiation and 
absorption measurements on diffusing 
materials are among the principles pre- 
sented in detail. 

Part II considers applications. In this 
part of the book, problems which have 


been attacked successfully by absorption 
techniques are briefly outlined. The 
results of several papers specifically con- 
cerned with the application are then dis- 
cussed. Typical absorption spectra are 
frequently given as illustrations. If the 
reader has a deeper interest in a particular 
problem, sufficient references are provided 
to guide additional study. Many of the 
applications are chosen from the fields of 
biochemistry and metallurgy and, hence, 
consider complex molecules. In the usual 
jet propulsion and astronautical problems, 
the engineering scientist is primarily con- 
cerned with monatomic, diatomic and a 
few triatomic molecules. Although the 
book does not discuss specific applications 
selected from aeronautical and space 
sciences, it does clearly demonstrate the 
principles and techniques. 

Part III is concerned with the tech- 
niques. To cover wavelengths from 0.1 
to 1004 requires a wide variety of instru- 
ments; these instruments are discussed. 
The principles of operation of modern in- 
struments are described. Accessories such 
as cells, light sources, and standard refer- 
ence substances are included in Part ITI. 

The cross references between different 
parts do enhance the value of the book in 
spite of the many shifts from page to page 
required to fully utilize the cross refer- 
encing. Although the book is primarily 
descriptive, formulae are brought into the 
discussion as necessary. Units and sym- 
bols are clearly stated and defined. In 
reading the book, the reviewer noticed 
very few typographical errors. Drawings 
and plates are used generously throughout 
the book. 
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Experimental Investigationin an Altitude 
Test Facility of Burning of Excess Com- 
bustibles in a Rocket Engine Exhaust, by 
Harry E. Bloomer, Paul E. Renas and 
Robert J. Antl, NASA TN D-200, Jan. 
1960, 25 pp. 

Full-Scale Propulsion Testing in Wind 
Tunnels, by R. W. Hensel and H. K. Matt, 
Arnold Engng. Dev. Center, AEDC-TR-59- 
9, May 1959, 80 pp. 

Combustion and Propulsion, Third 
AGARD Colloquium (Palermo, 1958); 
Noise-Shock Tubes—Magnetic Effects— 
Instability and Mixing, Pergamon Press, 
N. Y., 1958, 614 pp. 

Development Problems in Large 
Liquid Rocket Engines, by R.S. Levine, 
pp. 3-13. 

Propulsion by Air Breathing Engines, 
by A. A. Lombard, pp. 24-41. 

Effects on Turbojet Combustors and 
Afterburners of Other Engine Com- 
ponents, by J. H. Child, pp. 55-77. 


The Problem of the Ducted Rocket, 
by E. Le Grivzés, pp. 86-113. (In 
French.) 


Jet Engine Noise Reduction, by N. D. 
Sanders and W. J. North, pp. 185-196. 

Some Thought on Noise Suppression 
Nozzle Design, by E. J. Richards, pp. 
197-222. 

Muffling of Aviation Engine Stands; 
Technical, Financial, and Physiological 
Aspects, by P. Ducrot, A. Riehn, G. 
Roumilhac and H. Souvras, pp. 224-262. 
(In French.) 

Role of Fundamental Research in the 
Study of Inlets, by C. Fouré, pp. 333- 

(In French.) 

Some Problems of Providing for Scien- 
tific Research on Rockets, by A. A. 
Blagonravov, ARS Journat, vol. 30, Jan. 
1960, pp. 22-26. 

Analysis of Ramjet Engine Performance 
Including Effects of Component Changes, 
by Richard J. Weber and Roger W. Lui- 
dens, NACA Res. Mem.E56D2-, Oct. 1956, 
46 pp. (Declassified from Confidential by 
authority of NASA TPA 8, 7-22-59.) 

Effect of Immersed Surfaces in Com- 
bustion Zone on Efficiency and Stability 
of 5-inch-diameter Ram-jet Combustor, 
by Thaine W. Reynolds and Donald W. 
Male, NACA Res. Mem. E54C25, June 
1954, 35 pp. (Declassified from Confiden- 
tial by authority NASA TPA 8, 7-22-59.) 

Development Problems of Rocket En- 
gines for Ballistic Missiles, by T. F. Dixon, 
Interavia, vol. 14, July 1959, pp. 818-819. 
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Some Remarks Concerning a Magneto- 
hydrodynamic Method of Propulsion, by 
Rudolf X. Meyer, Space Tech. Labs., Inc., 
Phys. Res. Lab., GM-TR-0165-00544, Dec. 
1958, 16 pp. 
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John J. Newgard, ARS Journat, vol. 29, 
Aug. 1959, pp. 603-604. 
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Magnetic Field in a Fully Ionized Gas, 
y M.S. Sodha and Y. P. Varshni, Phys. 
Rev., vol. 114, May 15, 1959, pp. 946-947. 

Ionic Theory of Plasmas and Magneto- 
hydrodynamics, by H. S. Green, Phys. 
Fluids, vol. 2, no. 4, July-Aug. 1959, pp. 
341-349. 
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Propagation of Weak Hydrodynamic 
Discontinuities, by J. Bazer and 
Fleischman, Phys. Fluids, vol. 2, no. 4, 
July-Aug. 1959, pp. 366-378. 

Rayleigh’s Problem in Magnetohydro- 
dynamics, by C. C. Chang and J. T. Yen, 
Phys. Fluids, vol. 2, no. 4, July-Aug. 1959, 
pp. 393-403. 

Laminar Steady-state Magnetohydro- 
dynamic Flow in an Annular Channel, by 
Samuel Globe, Phys. Fluids, vol. 2, no. 4, 
July-Aug. 1959, pp. 404-407. 

Radiative Transport within an Ablating 
Body, by Leo P. Kadanoff, Avco-Everett 
Res. Lab., July 1959, 37 pp. 

Aerothermodynamic and_ Electrical 
Properties of Some Gas Mixtures to Mach 
20, by W. Chinitz, C. L. Eisen and R. A. 
Gross, ARS JouRNAL, vol. 29, Aug. 1959, 
pp. 573-578. 

Meteoric Ablation through Interfacial 
Instability, by Sin-I Cheng, ARS Jour- 
NAL, vol. 29, Aug. 1959, pp. 579-587. 

Graphical Computation of Shock and 
Detonation Waves in Real Gases, by 
George H. Markstein, ARS JourNat, vol. 
29, Aug. 1959, pp. 588-590. 

Magnetic Boundary Layers, by H. F- 
Ludloff, ARS Journat, vol. 29, Aug: 
1959, pp. 590-591. 

Optimum Magnetic Field fcr Stagna- 
tion Heat Transfer Reduction at Hyper- 
sonic Velocities, by Robert Goulard, 
ARS Journat, vol. 29, Aug. 1959, pp. 
604-605. 

Analysis of the Aerodynamic Heating 
for a Re-entrant Space Vehicle, by M. J. 
Brunner, J. Heat Transfer (ASME 
Trans.), vol. 81, Series C, Aug. 1959, pp. 
223-229. 

A Photographic Study of Boiling in the 
Absence of Gravity, by R. Siegel and C. 
Usiskin, J. Heat Transfer (ASME Trans.), 
vol. 81, Series C, Aug. 1959, pp. 230-236. 

A Note on Magnetohydrodynamic- 
hypersonic Flow Past a Blunt Body, by 
William A. Bush, J. Aero/Space Sct., vol. 
26, Aug. 1959, pp. 536-537. 


The Effect of Vehicle Deceleration on a 
Melting Surface, by D. M. Tellep, //. 
Aero/Space Sci., vol. 26, Aug. 1959, 
pp. 537-538. 

Solutions of the Laminar Compressible 
Boundary-Layer Equations with Trans- 
piration which are Applicable to the Stag- 
nation Regions of Axisymmetric Blunt 
Bodies, by John T. Howe and William A. 
Mersman, NASA TN D-12, Aug. 1959, 
73 pp. 

Turbulent Boundary Layer Equations 
with Binary Diffusion, by N. Mess and 
H. G. Lew, Gen. Electric Co., Missile & 
Space Vehicle Dept., T.I.S. R59SD310, 
a 1958, 35 pp. (Aerophys. Res. Memo. 


Compressible Flow Tables k = 1.40; 
Stagnation Temperature Change With- 
out Friction, by C. Hoebich, Army 
Rocket and Guided Missile Agency, Ord- 
nance Missile Labs. Division, Design and 
Development Lab., ARGMA TN-1H1N- 
5, June 1959, 89 pp. 

On the Stability of a Shear Layer, by 
J. Mankes, Calif. Inst. Tech., Jet Propul- 
sion Lab., Prog. Rep. 30-10, April 20, 
1959, 9 pp. 

The Electric Conductivity of a Partially 
Ionized Gas, by A. C. Pipkin, Maryland 
Univ., Inst. for Fluid Dynamics and 
Appl. Math, TN BM-170, April 1959, 

0 pp. 

On Series Expansions in Magnetic 
——— Number, by Vernon J. Rossow, 
NASA TN D-10, Aug. 1959, 44 pp. 

A Contribution to the Theory of Meteor 
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Ablation, by George W. Sutton, General 
Electric, Co., Missile and Space Vehicle 
Devt., R59SD331, March 1959, 
23 pp. (Aerophys. Res. Memo. 28.) 

Transient Temperature With and With- 
out Ablation in the Skin of a Ballistic 
Missile Entering the Earth’s Atmosphere, 
by Shih-Yuan Chen, Boeing Airplane Co., 
Document D7-1138, June 1959, 40 pp. 

On Some Problems in Radiative Heat 
Transfer, by Saul S. Abarbanel, Mass. 
Inst. Tech., Fluid Dynamics Res. Group, 
Rep. no. 59-1, April 1959, 217 pp. (AF 
OSR TN 59-531; ASTIA AD 216, 538.) 

Heat Transfer to a Vaporizing, Ablating 
Surface, by Richard G. Fledderman, J. 
Aero/Space Sct., vol. 26, Sept. 1959, pp. 
604-605. 

How Far Are We on Re-entry Cooling? 
by Kurt R. Stehling, Space Aeron., vol. 32, 
no. 3, Sept. 1959, pp. 48-45. 

A Shallow-liquid Theory in Magneto- 
hydrodynamics, by L. E. Fraenkel, Calif. 
Inst. Tech., Guggenheim Aeron. Lab., June 
1959, 74 pp. (AFOSR TN-59-563; 
ASTIA AD-216,758.) 

On the Mechanism of Sub-cooled Nu- 
cleate Boiling, by S. G. Bankoff, Calif. 
Inst. Tech., Jet Prop. Lab., Memo., 30-8, 
Feb. 5, 1959, 50 pp. 

Experimental Studies on Heat Transfer 
and Fluid Flow ‘qo Hydrogen), by 
M. E. Ellion, E. L. Geery, L. Green Jr. 
and W. R. Thompson, Aerojet General 
Corp., Rep. 1520, vol. 11, Nov. 1958, 38 pp. 
(WADC TR-58-610, vol. IT.) 

Selected Topics from the Theory gf Gas 
Flow at High Temperatures, VI: (ome 
Aspects of Particle Interaction in Gases, 
by J. M. Burgers, Maryland Unir., Inst. 
Fluid Dynam. and Appl. Math., TN BN 
176, June 1959, 40 pp. 

Heat Transfer from Non-isothermal 
Surfaces, by Helmut Wolf, Princeton 
Univ., Proj. Squid Tech. Rep. PUR-38-0 
(ASTIA AD 215,140), April 1959, 25 pp. 
(Available only on microcards. ) 

A Study of the Mean Thickness of the 
Liquid Film and the Characteristics of the 
Interfacial Surface in Annular, Two-phase 
Flow in a Vertical Pipe, by D. A. Char- 
vonia, Proj. Squid Tech. Rep. PUR-39-T-R 
(ASTIA AD 216,823), May 1959, 253 pp. 
(Available only on microcards.) 

Entropy Production and Pressure 
Waves, by Gerald Rosen, Princeton Univ., 
Proj. Squid Tech. Rep. PR-90-P (ASTIA 
AD 217,238), June 1959, 11 pp. (Available 
only on microcards. ) 

Hydrodynamic Studies cn the Effect of 
Entrained Gases on Injection and Atom- 
ization, by R. W. Mascolo, North American 
Aviation, Inc., Rocketdyne Div., Rep. R- 
1466, June 5, 1959, 69 pp., 38 refs. 

Use of the Revolving Arm Method for 
the Determination of the Stagnation Tem- 
perature in an Ionized Gas, by F. Marcel 
Devienne, André F. Roustan, Henri L. 
Bosc, Robert G. Clapier and Raymond F. 
Meslet, Lab. Medilerranéen de Recherches 
Thermodynamiques, Feb. 1959, 21 pp. 
(ASTIA AD 211,473.) 

The Large High Pressure Arc Plasma 
Generator: A Facility for Simulating Mis- 
sile and Satellite Re-entry, by P. H. Rose, 
W. E. Powers and D. Hritzay, Avco-Everett 
Research Lab., Res. Rep. 56, June 1959, 30 
pp. 

Interaction of Two-dimensional Inviscid 
Incompressible Jet Facing a Hypersonic 
Stream, by S. H. Lam, Princeton Univ., 
Dept. Aeron. Engng., Rep. 447, March 
1959, 23 pp. (AFOSR- TN-B9-274: ASTIA 
AD 212,7 70: 

An At and Theoretical Study 
of Quartz Ablation at the Stagnation 


Point, by Mac C. Adams, William E. 
Powers and Steven Georgiev, Auvco- 
Everett Research Lab., Res. Rep. 57, June 
1959, 28 pp. 

Experimental Investigation of Coaxial 
Jet Mixing of Two Subsonic Streams at 
Various Temperatures, Mach Number, 
and Diameter Ratios for Three Configura- 
tions, by Richard R. Burley and Lively 
Bryant, NASA Memo 12-21-58E, Feb. 
1959, 29 pp. 

The Application of Constant Property 
Solutions to Mass Transfer Cooling Calcu- 
lations, by C. J. Scott, Minnesota Univ. 
Rosemount Aeron. Labs. Engng. Mem. 76, 
Dec. 1958, 31 pp. (AFOSR-TN-59-201; 
ASTIA AD 211,521.) 

Impact of High Speed Particles with 
Solids and Gases, by M. C. Kells and D. 
D. Keough, Stanford Res. Inst., Poulter 
Labs. Tech. Rep. 017-58, Dec. 15, 1958, 42 
pp. 

Burnout of Heating Surfaces in Water, 
by R. L. Menegus, AEC, DP-363, March 
1959, 52 pp. 

Nitrogen at High Temperatures (Pre- 
sented at ASME High Temperature 
Symposium, Purdue Univ., Feb. 1959), by 
Frank Martinek, General Electric Co., 
Flight Propulsion Lab., Cincinnati, Ohio, 
1959, 50 pp. (AFOSR-TN-59-264; 
ASTIA AD 212,467.) 

Radio-frequency Circuits for Plasma 
Physics, by R. L. Gamblin, Electronics, 
vol. 32, July 3, 1959, pp. 50-52. 

Entrance Effects and Prediction of Tur- 
bulence in Non-Newtonian Flow, by D. C. 
Bogue, Ind. Engng. Chem., vol. 51, July 
1959, pp. 874-878. 

Heat Transfer to Turbulent Non-New- 
tonian Fluids, by A. B. Metzner and P.S 
Friend, Ind. Engng. Chem., vol. 51, July 
1959, pp. 879-884. 

Laminar Heat Transfer Around Blunt 
Bodies in Dissociated Air, by Nelson H. 
Kemp, Peter H. Rose and Ralph W. Detra, 
J. Aero/Spac? Sci., vol. 26, July 1959, pp. 
421-430. 

Experimental Verification of Boundary 
Layer Corrections in Hypersonic Nozzles, 
by Donald L. Baradell, J. Aero/Space Sc7., 
vol. 26, July 1959, pp. 454-455. 

A Total-temperature Probe for High- 
temperature Boundary-layer Measure- 
ments, by M. Sibulkin, J. Aero/Space 
Sci., vol. 26, July 1959, pp. 458-459. 

Transport Phenomena in a Fluid Com- 
posed of Diatomic Molecules, by John S. 
Dahler, J. Chem. Phys., vol. 30, June 
1959, pp. 1447-1475. 

The Kinetic Radiation Equilibrium of 
Air, by P. Harteck, J. Phys. Chem., vol. 63, 
June 1959, pp. 956-961. 

Strong Hydromagnetic Disturbances in 
a Collision-free Plasma, by J. W. Dungey, 
Philosophical Magazine, vol. 4, May 1959, 
pp. 585-593. 

Stable Orbits of Charged Particles in an 
Oscillating Electromagnetic Field, by 
Erich S. Weibel, Phys. Rev., vol. 114, April 
1, 1959, pp. 18-23. 

Theory of Spatially Growing Plasma 
Waves, by Masao Sumi, Phys. Soc. 
Japan, J., vol. 14, May 1959, pp. 653- 
663. 

Flash X-ray Radiography for the Density 
Measurement in a Hypersonic Air Flow, 
by Fumio Tamaki and Chol-Soo Kim, 
Phys. Soc. Japan, J., vol. 14, May 1959, 
pp. 664-677. 

A Visualizing Method of Electron Tem- 
— and Plasma Density, by Hajime 

amagawa and Junji Fujita, Phys. Soc. 


Japan, J., vol. 14, May 1959, pp. 678- 
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Relativistic Hydrodynamics for a 
Charged Nonviscous Fluid, by Chau-Chin 
Wei, Phys. nt vol. 113, March 15, 1959, 


Trajectory Design Methods, Calif. Inst. 
Tech., Jet Prop. Lab., Pub. 30-18, Nov. 
1959, 28 pp. 

The Approximate Analysis of Zero-lift 
Trajectories, by C. P. Hoult, Air Force, 
Cambridge Res. Center, Geophys. Res. 
Directorate, GRD Res. Notes 17, AFC RC- 
TN-59-404, Aug. 1959, 19 pp. 

Effects of Small Changesin Position and 
Velocity at Burnout on the Orbit of a 
Ballistic Missile, by Edouard E. Perret 
and John C. Corbin Jr., Wright Air Dev. 
Center, TN 59-336, Oct. 1959, 25 pp. 

Trajectory Problems in Cislunar Space, 
Westinghouse Electric Corp., Air Arm Div., 
TN, Dec. 1959, 42 pp. (AF OSR TN 59- 
1284.) 

Orbits of Cosmic Rockets Toward the 
Moon, by L. I. Sedov, ARS JourNAL, vol. 
30, Jan. 1960, pp. 14-21. 

_ An Approximation to the Re-entry Tra- 
jectory, by Mildred M. Moe, ARS Jour- 
NAL, vol. 30, Jan. 1960, pp. 50-53. 

Analytical Estimates for Optimum 
Transfer Paths, by H. Munick and R. 
McGill, ARS JourNaL, vol. 30, Jan. 1960, 
pp. 120-121. 

Projected Orbits of 24-hour Earth 
Satellites, by Louis B. Wadel, ARS Jour- 
NAL, vol. 30, Jan. 1960, pp. 122-123. 

Satellite Librations of Large Amplitude, 
by W. B. Klemperer, ARS JourNAL, vol. 
30, Jan. 1960, pp. 123-124. 

Librations on a Slightly Eccentric Orbit, 
by Robert M. L. Baker Jr., ARS JourNaAL, 
vol. 30, Jan. 1960, pp. 124-126. 

Plane Libration of a Prolate Ellipsoidal 
Shell, by Robert M. L. Baker Jr., ARS 
JOURNAL, vol. 30, Jan. 1960, pp. 126-128. 


Rendezvous Compatible Orbits, by 
Norman V. Petersen and Robert S. Swan- 
son, Astron. Sci. Rev., vol. 1, Oct.-Dec. 
1959, pp. 13-14, 20. 

An Orbital Rendezvous Simulator, by 
E. Levin and J. Ward, Astron. Sci. Rev., 
vol. 1, Oct.-Dec. 1959, pp. 15-16, 21. 

Orbital Plane-change Maneuver, by 
Daniel B. DeBra, Astron. Sci. Rev., vol. 1, 
Oct.-Dec. 1959, pp. 19-20. 

On an Aspect of the Relations Between 
Mechanics and Artificial Earth Satellites, 
by Maurice Roy, La Recherche Aéro- 
nautique, no. 73, Nov.-Dec. 1959, pp. 
1-7. (In French.) 


Motion of a Satellite in the Earth’s 
Gravitational Field, by G. V. Groves, 
Proc. Roy. Soc., London, vol. 254, no. 1276, 
Jan. 19, 1960, pp. 48-65. 


A Method for Simulating the Atmos- 
heric Entry of Long-range Ballistic 

issiles, by A. J. Eggers Jr., NACA Rep. 
1378, 1958, 7 pp. 

Experimental Investigation of the Simu- 
lation of Atmospheric Entry of Ballistic 
Missiles, by Stanford E. Neice, James A. 
Carson and Bernard E. Gunningham, 
NACA Res. Mem. A57126, Dec. 1957, 32 
pp. (Declassified from Confidential by 
authority NASA TPA 8, 7-22-59.) 

Lagrange Multipliers and Quasi-steady 
Flight Mechanics, by Angelo Miele, J 
Aero/Space Sci., vol. 26, Sept. 1959, pp. 
592-598. 

Steady-state Circular Pitching and Yaw- 
ing Motion of Symmetric Missiles, by 
Charles H. Murphy Jr., J. Aero./Space 
Sci., vol. 26, Sept. 1959, pp. 600-601. 
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Similitude of Hypersonic Real-gas 
Flows over Slender Bodies with Blunted 
Noses, by Hsien K. Cheng, J. Aero/Space 
Sci., vol. 26, Sept. 1959, pp. 575-585. 

Stability of a Globular Cluster in the 
Oblate Earth’s Equatorial Plane, by 
Marvin H. Hewitt, ARS JournaAt, vol. 29, 
Aug. 1959, pp. 595-596. 

Ogive Having Minimum Drag in Two- 
dimensional Hypersonic Flow, by E. 
Mattioli, Aerotecnica, vol. 29, no. 1, Feb. 
1959, pp. 16-24. (In Italian.) 

Optimum Staging Techniques, by Leon 
Weisbord, ARS Journat, vol. 29, June 
1959, pp. 445-446. 

Some Design Aspects of an Interplane- 
tary Exploratory Vehicle, by Henry K. 
Hebeler and Richard D. White, Astro- 
nautical Sci. Rev., vol. 1, no. 2, April-June, 
1959, pp. 11-15. 

Graphic Determination of the Main In- 
fluences and Design Parameters of Mis- 
siles, by O. Scholze, Raketentechnik und 
Raumfahrtforschung, vol. 3, no. 2, April 
1959, pp. 36-44. (In German.) 

Design of Missile Bodies for Minimum 
Drag at Very High Speeds—Thickness 
Ratio, Lift, and Center of Pressure Given, 
by T. Strand, J. Aero/Space Sci., vol. 26, 
Sept. 1959, pp. 568-570. 

On a Class of Variational Problems in 
Rocket Flight, by G. Leitmann, J. Aero/- 
Space Sci., vol. 26, Sept. 1959, pp. 586- 
591. 

Theoretical Principles of Guided Mis- 
sile Systems, by E. G. C. Burt, J. Roy. 
A Soc., vol. 63, Aug. 1959, pp. 455- 
464. 


Vehicle Design, Testing and 
Performance 


Static Stability and Separation Charac- 
teristics of a Two-stage Rocket Configura- 
tion at Mach Numbers from 1.57 to 4.50, 
by Kenneth L. Turner, David S. Shaw 
and Lawrence W. Enderson Jr., NASA 
TN-D-188, Jan. 1960, 35 pp. 

Proposal for Standard Weight Defini- 
tions for Space Vehicles, by H. H. Koelle, 
Army Ballistic Missile Agency, Redstone 
Arsenal, Rep. DSP-TR-6-59, Dec. 1959, 
27 pp. 

Booster Design for a 20,000 Ib. Satellite, 
by S. B. Kramer and R. A. Byers, Space 
Aeron., vol. 31, June 1959, pp. 52-55. 

Atlas, Carrier Rocket for Important 
Space Flight Projects, by D. E. Kédlle, 
Raketentechnik und Raumfahrtforschung, 
vol. 3, no. 2, April 1959, pp. 50-52. (In 
German.) 

The First Artificial Planetoid Mechta, by 
D. E. KGlle, Raketentechnik und Raum- 
fahrtforschung, vol. 3, no. 2, April 1959, 
pp. 50-53. (In German.) 

Satellite Vanguard II with Infrared De- 
tectors, Raketentechnik und Raumfahrt- 
forschung, vol. 3, no. 2, April 1959, pp? 
55-56. (In German.) 

Space Sound Juno II, by D. E. Kédlle, 
Raketentechnik und Raumfahrtforschung, 
vol. 3, no. 2, April 1959, pp. 57-59. (In 
German.) 

Designing an ICBM Nose-cone Liner, 
by Bernard H. Gerberg, ASTRONAUTICS, 
vol. 4, Sept. 1959, pp. 38-39, 97. 

Launcher for Honest John, by Arnold 
A. Kester, Ordnance, vol. 43, May-June 
1959, p. 1008. 

Rocket Readiness, by J. W. Wiggins, 
Ordnance, vol. 43, May-June 1959, p. 1010. 

The Explorers, by W. von Braun, Astro- 
nautica Acta, vol. 5, no. 2, 1959, pp. 126- 
143. 


_ Lockheed X-17 Rocket Test Vehicle and 
Its Applications, by Ronald Smelt, ARS 
JOURNAL, vol. 29, Aug. 1959, pp. 565- 


Some Considerations Concerning Mag- 
netohydrodynamic Exhaust Control for 
Rocket Guidance, by Wulf B. Kunkel, 
Space Tech. Labs., Inc., Phys. Res. La). 
T R-59-0000-00608, Feb. 1959, 12 pp. 

Requirements of Trajectory Corrective 
Impulses During the Approach Phase of 
an Interplanetary Mission, by Alan |... 
Friedlander and David P. Harry III, 
NASA TN D-255, Jan. 1960, 53 pp. 

A Study of the Guidance of Space Ve- 
hicle Returning to a Braking Ellipse About 
the Earth, by Jack A. White, NASA 
TN-D-191, Jan. 1960, 65 pp. 

Midcourse Guidance Problem in Sate!- 
lite Interception, I, by Angelo J. Skala- 
furis and Donald H. Schiller, ARS Jovx- 
NAL, vol. 30, Jan. 1960, pp. 41-45. 


Missile Orientation Errors in Command 
Guidance Systems, by Kenneth ©. 
Mathews, ARS Journat, vol. 30, Jan. 
1960, pp. 46-49. 

Steady-state Behavior of Extended 
Surfaces in Space, by J. W. Tatom, ARS 
JOURNAL, vol. 30, Jan. 1960, pp. 118-119. 


Transients in Optimal Control Systems, 
by O. M. Kryzhanovski and V. M. Kuntse- 
vich, ARS Journal, vol. 30, Jan. 1960, pp. 
86-92. 

Contribution to the Theory of Inertial 
Damped Systems with Arbitrary Period, 
Invariant with Respect to Maneuver of the 
Object, by V. A. Bodner, V. P. Seleznov 
and V. E. Ovcharov, ARS Journat, vol. 
30, Jan. 1960, pp. 93-97. 

Effect of Thrust Misalignment on Gyro- 
stabilized Vehicles, by Russell P. Nagor- 
ski, ARS JourNAL, vol. 30, Jan. 1960, pp. 
128-129. 

Terminal Guidance for Rendezvous, by 
Raymond S. Wiltshire, Astron. Sci. Rev., 
vol. 1, Oct.-Dec. 1959, pp. 9-10. 

Multifrequency Pulse Devices for Tele- 
control, by F. A. Katkov, Automation and 
Remote Control (trans. of Avtomatika 1 
Telemekhanika), vol. 20, Jan. 1959, pp. 
49-55. 

Auto-oscillations in a Single-loop Auto- 
matic Control System Containing Two 
Symmetric Relays, by Tu Syui-Yan’ and 
Tei Lui-Vy, Automation and Remote Con- 
trol (trans. of Avtomatika i Telemekhanika), 
vol. 20, Jan. 1959, pp. 83-88. 

The Behavior of a Vibration Servo in an 
Electrical eg by I. N. Drutova, 
Automation and Remote Control (trans. of 
Avtomatika i Telemekhanika), vol. 20, 
Feb. 1959, pp. 113-121. 

The Stability of a Nonlinear Control 
System with a Neutral Object, by L. I. 
Kupriyanova, Automation and Remote 
Control (trans. of Avtomatika i Tele- 
mekhanika), vol. 20, Feb. 1959, pp. 122- 
129. 

Unified Twin-drive Servosystems, by 
P. F. Klubnikin, Automation and Remote 
Control (trans. of Avtomatika i Tele- 
mekhanika), vol. 20, Feb. 1959, pp. 154- 
168. 


Guidance and Control 


Choice of an Optimum Amplification 
Factor for a Self-aligning Control System, 
by I. I. Perel’man, Automation and Re- 
mote Control (trans. of Avtomatika i 
Telemekhanika), vol. 20, Feb. 1959, pp. 
177-184. 

Synthesis of the Control Unit of a Servo- 
system Optimum from the Point of View 
of Speed of Response, by Sung Chien, 
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Automation and Remote Control (trans. of 
Avtomatika i Telemkhanika), vol. 20, 
March 1959, pp. 261-275. 

Concerning the Realization of Optimum 
Weighting Functions of Pulse Servosys- 
tems, by V. P. Perov, Automation and 
Remote Control (trans. of Avtomatika i 
Telemekhanika), vol. 20, March 1959, pp. 
276-284. 

A Fast Transistorized Magnetic Push- 
pull Amplifier for a Servo Drive, by V. 8. 
Volodin, E. D. Larin, M. A. Rozenblat, 
and G. V. Subbotina, Automation and 
Remote Control (trans. Avtomatika i 
Telemekhanika), vol. 20, March 1959, pp. 
307-313. 

Choice of Characteristic for Nonlinear- 
rate Feedback in Position Servos, by 
B. N. Naumov, Automation and Remote 
Control (trans. of Avtomatika i Tele- 
mekhanika), vol. 20, March 1959, pp. 319 
329. 

Use of an Integral Square Estimate to 
Determine the Optimal Parameters of an 
Autopilot with Rate Feedback, by V. D. 
Matytsin and V. A. Ryapolov, Automa- 
tion and Remote Control (trans. of Avto- 
matika t Telemekhanika), vol. 20, April 
1959, pp. 403-408. 

The Dynamics of a Vibration Loop in an 
Electrical Servo During Free Oscillations. 
by I. N. Krutova, Automation and Remote 
Control (trans. of Avtomatika i Tele- 
mekhantka), vol. 20, April 1959, pp. 409- 
420. 


Compensation of a Continuous Auto- 
matic Control System by Means of a De- 
lay-Element Filter, by Wang Sin-Min, 
Automation and Remote Control (trans. of 
Avtomatika « Telemekhanika), vol. 20, 
April 1959, pp. 453-464. 


Stability of Nonlinear Automatic Con- 
trol Systems, by O. I. Komarnitskaia, 
PMM: Journal of Appl. Math. and 
Mech. (trans. of Prikladnaya Matematika i 
— vol. 23, no. 3, 1959, pp. 716- 

The Stability of Motion of a Gyroscope, 
by Chzhan Sy-In, PMM: J. Appl. 
Math. and Mech. (trans. of Prikladnaya 
Matematika i Mekhanika), vol. 23, no. 3, 
1959, p. 860. 

The Stability of Unperturbed Motion of 
a Gyroscope, by Chzhan Sy-In, PMM: 
J. Appl. Math. and Mech. (trans. of 
Prikladnaya Matematika i Mekhanika), 
vol. 23, no. 3, 1959, pp. 860-863. 

Reliability in Inertial Systems, by A. M. 
Voutsas and W. T. Chow, /nteravia, vol. 
14, July 1959, p. 823. 

LINS—Lightweight Inertial Navigation 
System, Jnteravia, vol. 14, July 1959, p. 
824, 

Comparison Study Between a Micro- 
wave Interferometer and an Electrical 
Contact System for Following the Motion 
of High-Speed Projectile, by Z. I. Slawsky, 
A. E. Siegel, L. T. Ho, H. J. Kopp and J. T. 
Vanderslice, Rev. Sci. Instr., vol. 30, Aug. 
1959, pp. 679-683. 

MIG Platform Extends Inertial System 
Use (Design Digest), by Bernard Kovit, 
Space/Aeron., vol. 32, July 1959, pp. 128- 
129. 

Precision Gyroscopics—New Demands 
on an Old Art, by Roland E. Barnaby Jr., 
Sperry Engineering Review, vol. 12, June 
1959, pp. 2-7. 

Determining Lunar Mission Guidance 
Requirements, by Benjamin M. Rosen, 
Sperry Engineering Review, vol. 12, June 


1959, pp. 27-38. 

Trade-Off Considerations in the Design 
of Guidance Equipment for Space Flight, by 
Charles J. Mundo Jr., Aero/Space Engng., 
vol. 18, June 1958, pp. 31-34. 


The Radar Appearance of a Rotating 
Corner Reflector, by R. J. Davis, Astro- 
nautica Acta, vol. 5, no. 1, 1959, pp. 1-14. 


Design Digest: Astro Compass Brings 
Space Navigation Closer, by Bernard 
Kovit, Space/Aeron., vol. 31, June 1959, 
pp. 127-129. 


The Performance of Inertial Compo- 
nents on an Unstabilized Base, by Henry F. 
Blazek, ARS Preprint, 775-59, April 30- 
May 1, 1959, 12 pp. 

Tracking and Recovery of the Jupiter 
Rocket Nose, Raketentechnik und Raum- 


fahriforschung, vol. 3, no. 1, Jan.-March, 


1959, pp. 25-26. (In German.) 


Instrumentation and 


Communications 


Theory of the Stagnation Point Lang- 
muir Probe, by L. Talbot, Univ. of Calif., 
Inst. Engng. Res., Tech. Rep. HE-150-168, 
March 30, 1959, 26 pp. (AFOSR-TN-59- 
474; ASTIA AD 215,717.) 


Commutators for Airbourne Multiplex 
Telemetering, by M. M. Kranzler, Elec- 
tronics, vol. 32, July 3, 1959, pp. 46-47. 


Six Rules for Evaluating CFVS Gen- 
erators, by F. T. DeWolf, E. A. Erdely, 
E. E. Kolatorowicz and W. R. Miller, 
Space/Aeron., vol. 32, July 1959, pp. 93- 
94, 96, 98, 102-103. 

The Payoff of Missile Design: Elec- 


Hypersonic Flow 
Rarefied Gas Phenomena 
Transport Phenomena 


tists in these fields. 


one or more of these areas. 


RESEARCH SCIENTISTS 


ARMOUR RESEARCH FOUNDATION has 
professionally challenging openings for high calibre 
research engineers and scientists to do theoretical 
and experimental research in the fields of: 


Shock and Detonation Wave Phenomena 


The Foundation encourages original research 
and publication of papers. These positions offer 
stimulating associations with other leading scien- 


Applicants must have an M.S. or Ph.D. degree 
in Engineering or Science, and have experience in 
In addition to the 
many cultural opportunities in Chicago, we offer 


SOUTHWEST 


SELF-ALIGNING BEARINGS 


PLAIN TYPES ROD END 
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PATENTED U.S. A. 
$ World Rights Reserved 
$ CHARACTERISTICS 
ANALYSIS RECOMMENDED USE 


and Race 
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Chrome Alloy Steel Ball 
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temperature (800-1200 degrees F.). 

For types operating under high radial 
ultimate loads (3000-893,000 Ibs.). 

For types operating under normal loads 


the advantage of being located on a university Steel Ba with minimum friction requirements. 
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tronic Fuzing, by William H. Doty, Space/- 
Aeron., vol. 32, July 1959, pp. 120-123. 

A Device for Dynamical Measurements 
of Pressure, by P. Bassi, R. Cano, §. 
Focardi, C. Rubbia, A. Micheline and F. 
Saporetti, Nuovo Cimento, vol. II, no. 4, 
Feb. 16, 1959, pp. 589-594. 

Reliable Optical Window for Cells 
Subjected to Widely Cycling Tempera- 
tures and Pressures, by W. C. Waggener 
and A. M. Tripp, Rev. Sci. Instr., vol. 30, 
Aug. 1959, pp. 677-680. 

Pressure Sensing Probes for Detecting 
Shock Waves, by Adolph B. Amster and 
Raymond L. Beauregard, Rev. Sci. Instr., 
vol. 39, Oct. 1959, pp. 942-944. 


High-sensitivity Pressure Gauge Coup- 
ling Two Spoon Gauges, by R. R. Mon- 
champ, E. Bannister and F. A. Cotton, 
Rev. Sci. Instr., vol. 39, Oct. 1959, pp. 
945-946. 

Basic Principles of Radar with Particular 
Reference to Aircraft and Missile Applica- 
tions, by F. D. Boardman, J. Roy. Aeron. 
Soc., vol. 63, Oct. 1959, pp. 581-596. 


Photographic Pyrometry of Rocket Ex- 
haust Jets, by Fredrick S. Simmons ry 
Arthur G. DeBell, Aircr. Engng., vol. 
May 1959, pp. 144-149. 


Low Frequency Vibrating Table for the 
Dynamic Calibration of Accelerometers, 
by L. Lemoine, La Recherche Aéronau- 
tique, no. 69, March-April 1959, pp. 27-37. 

Liquid Nitrogen Liquid Level Controller, 
by Rodger L. Gamblin, Ernest Goldberg 
and Dane T. Scag, Rev. of Sci. Insir., vol. 
30, May 1959, pp. 371-373. 

What the Designer Faces in Active 
Satellite Communications, Space/Aeron., 
vol. 31, June 1959, pp. 122-125. 

Representation and Analysis of Signals, 
Part I: The Use of Orthogonalized Ex- 
ponentials, by W. H. Huggins, Air Force 
Cambridge Res. Center, Tech. Rep. AFCRC 
TR 57-357, Sept. 30, 1957, 76 pp. 

Temperature Response of a Hot-wire 
Anemometer to Shock and Rarefaction 
Waves, by S. G. Datar, Inst. Aerophys., 
Univ. Toronto, UTIA TN 28, June 1959, 
20 pp. 

Study of a Phosphor Light Pulser, by 
M. R. Zatzick, Air Force Cambridge Res. 
Center, Geophys. Res. Directorate, AFC RC- 
T R-58-274, Sept. 1959, 19 pp. 

Study of a Pulsed Logarithmic Photome- 
ter, by M. R. Zatzick, Air Force Cam- 
bridge Res. Center, Geophys. Res. Direc- 
torate, Sept. 1959, 25 pp. 


General Model of a Semi-infinite 
Transmission Line with Mutual Induct- 
ance, by Richard H. Huddlestone, Space 
Tech. Labs., Inc., Phys. Res. Lab., T R-59- 
0000-00786, Oct. 1959, 16 pp. 


Survey of Analog Cross-spectral Ana- 
lyzers, by Lester J. Seltzer and Duane T. 
McRuer, Wright Air Dev. Center, Tech. 
Rep. 59-241, Oct.1959, 72 pp. 

The Optimal Parameters for Noise 
Stability in a Telemetry System, by V. A. 
Kashirin, Automation and Remote Control 
(trans. of Avtomatika i Telemekhanika), vol. 
20, Feb. 1959, pp. 216-225. 

Noise Stability in Pulse-width and 
Pulse-time Telemetry Subject to Strong 
Fluctuation Noise, by N. V. Pozin, Auto- 
mation and Remote Control (trans. of 
Avtomatika i Telemekhanika), vol. 20, 
Feb. 1959, pp. 227-235. 

Investigation of the Characteristics of 
Pneumatic Jet Elements, by L. A. Zal- 
monzon and A. I. Semikova, Automation 
and Remote Control (trans. of Avtomatika i 
— vol. 20, April 1959, pp. 
431-44 
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Strain Gauge with High Sensitivity 
Made of Semiconductor Material, Avia- 
tion Week, vol. 72, Feb. 1, 1960, pp. 74-75. 

Sensitive Method for the Measure- 
ment of Small Rotations, by J. King and 
J. W. Gates, J. Sci. Instr., vol. 36, Dec. 
1959, pp. 507-508. 


Simple Method for the Control of the 
Liquid Level in a Manometer, by J. G. 
Cunningham, J. Sci. Instr., vol. 36, Dec. 
1959, pp. 509. 


Effect of Antenna Size on Gain, Band- 
width, and Efficiency, by Roger F. Har- 
rington, J. Res. Nat. Bur. Standards, D 
Radio Propagation, vol. 64D, no. 1, Jan.- 
Feb. 1960, pp. 1-12. 


A Very-low-frequency Antenna for In- 
vestigating the Ionosphere with Horizon- 
tally Polarized Radio Waves, by R. S. 
Maemillan, W. V. T. Rusch and R. M. 
Golden, J. Res., Nat. Bur. Standards, D. 
Radio Propagation, vol. 64D, no. 1, Jan.- 
Feb. 1960, pp. 27-36. 

Sampling Oscilloscope for Millimicro- 
second Pulses at a 30-Mc Repetition Rate, 
by Arnold S. Farber, Rev. Sci. Instr., vol. 
31, Jan. 1960, pp. 15-17. 

Precise Measurement of Density and 
Surface Tension at Temperatures up to 
1000°C in One Apparatus, by George J. 
Janz and Max R. Lorenz, Rev. Sci. Instr., 
vol. 31, Jan. 1960, pp. 18-22 


Fluctuation Noise in Frequency-modu- 
lated Radio Lines Transmitting a Multi- 
channel Signal with Pulse-phase Modula- 
tion, by G. A. Malolepshii, Radio Engng. 
(trans. of Radiotekhnika), vol. 13, no. 7, 
Nov. 1959, pp. 12-30. 

Some Problems of Magnetic Pulse Re- 
cording, by R. G. Ofengenden, Radio 
Engng. (trans. of Radiotekhnika), vol. 13, 
no. 7, Nov. 1959, pp. 95-102. 


Atmospheric and Space Physics 


Atmospheric Refraction of Infrared Ra- 
diation, by T. P. Condron, Air Force 
Cambridge Research Center, Geophys. Res. 
Directorate, GRD Res. Notes 23, AFCRC- 
TN-59-617, Oct. 1959, 15 pp. 

Results of Scientific Investigations 
Made by Soviet Sputniks and Cosmic 
Rockets, by V. I. Krasovskii, ARS Jour- 
NAL, vol. 30, Jan. 1960, pp. 27-33. 


Summary of Rawinsonde Measure- 
ments of Temperatures, Pressure Heights, 
and Winds Above 50,000 Feet Along a 
Flight-test Range in the Southwestern 
United States, by Terry J. Larson and 
Harold P. Washington, "NV ASA TN D-192, 
Jan. 1960, 80 pp. 

Interplanetary Gas, 1: Hydrogen Radia- 
tion in the Night Sky, by John C. Brandt 
and Joseph W. Chamberlain, Astrophys. J., 
vol. 130, no. 2, Sept. 1959, pp. 670-676. 


Carbon Thermonuclear Reactions and 
the Formation of Heavy Elements, by 
A. G. W. Cameron, Astrophys. J., vol. 130, 
no. 2, Sept. 1959, pp. 429-434. 

Advances in Astron. Sci., vol. 4 (Proc., 
5th Annual Meet., American Astron. Soc., 
Nov. 1958), Plenum Press, N. Y., 1959: 

Some Aspects of the Electrical Proper- 

ties of the Upper Atmosphere, by P. A. 

Goldberg, pp. 37-50. 

On Gravitation, by Vjekoslav Gradecak, 

pp. 98-113. 

Corpuscular Radiation Around the 

Earth, by James A. Van Allen, pp. 329- 

334. 

Artificial Modification of the Earth’s 

a Belt, by 8. Fred Singer, pp. 


Interplanetary Dust Distribution and 
Erosion Effects, by David B. Beard, 
Advances in Astron. Sci., vol. 3 (Proc. 
Western Regional Meet., American Astron. 
Soc., Aug. 1958), Plenum Press, N. 
1958, Paper 23, 6 pp. 

Advances in Astron. Sct., vol. 2 (Proc., 
4th Annual Meet., American Astron. 
aaa Dec. 1957), Plenum Press, N. Y., 
1958: 


A Recoverable Emulsion Package for 

Extra-atmospheric Study of Cosmic 

Rays, by Robert A. Webster, Paper 32, 

12 pp. 

Plastic Balloons for Planetary Research, 

by Malcolm D. Ross, Paper 42, 13 pp. 

A Preliminary Determination of the 

Positions of Some First, Second, and 

Third Magnitude Stars on Mars’ 

Celestial Sphere, by Frederick West J»., 

Paper 44, 31 pp. 

Meteor Trains, by Robert F. Hugés, 
Smithsonian Contribs. to Astrophys., vol. 3, 
no. 8, 1959, pp. 79-94. 

Exploring the Depth of the Surface 
Layer of the Moon from a Radar Space 
Observatory, by W. E. Fensler and K. MI. 
Siegel, Aero/Space Engng., vol. 18, Nov. 
1959, pp. 38-41. 

Should the Laws of Gravitation be Re- 
considered? Part II: Experiments in 
Connection with the Abnormalities Noted 
in the Motion of the Paraconical Pendulum 
with an Anisotropic Support, by Maurice 
F. C. Allais, Aero/Space Engng., vol. 18, 
Oct. 1959, pp. 51-55. 

The Study of the Terrestrial Corpsucular 
Radiation and Cosmic Rays During the 
Flight of a Cosmic Rocket, by S. N. 
Vernov, A. E. Chudakov, P. V. Vakulov 
and Yu. I. Logac hey, Soviet Physics- 
es vol. 4, no. 2, Oct. 1959, pp. 338- 


= Nonlinear Vacuum Effect in Gravita- 
tion Theory, by N. V. Mitskevich, Soviet 
Physics-J ETP, 36(9), no. 4, Oct. 
1959, pp. 859-86 

Atomospheric by A.S. Monin, 
Soviet Physics-Uspekhi, vol. 2(67), no. 1, 
Jan.-Feb. 1959, pp. 50-58. 

Cosmic Ray Detection by Visual Scintilla- 
tions, by Herman Yagoda, Jnst. Aeron. 
Sci., Rep. 59-45, Jan. 1959, 8 pp. 

Infra-red Spectrescopy from Balloons 
and the Possibility of Some Observations 
on the Biosphere, by David M. Gate, 
Inst. Aeron. Sci., Rep. 59-47, Jan. 1959 5 
pp. 

Vistas in Astronautics, vol. 2 (2nd 
Annual Astron. Symposium), Pergamon 
Press, N. Y., 1959, 318 pp: 

How High the Moon, by C. C. Furnas, 

pp. 9-17. 

Cosmic Debris of Interplanetary Space, 

by Maurice Dubin, pp. 39-45. 

The Scientific and Engineering Ex- 

ploration of the Moon, by John L. 

Barnes, pp. 255-260. 

Dust on the Moon, by Thomas Gold, 

pp. 261-266. 

On the Lunar Dust Layer, by Fred L. 

Whipple, pp. 267-272. 

The Exploration of the Moon, by Gerald 

P. Kuiper, pp. 273-313. 

High-energy Neutrons in Cosmic Rays, 
by G. N. Flerov, V. I. Kalashnikova, A. V. 
Podgurskaya, E. D. Vorob’ev and G. A. 
Stolyarov, Soviet Physics-JETP, vol. 
36(9), no. 3, Sept. 1959, pp. 511-515. 

Should the Laws of Gravitation Be Re- 
considered? PartI: Abnormalities in the 
Motion of a Paraconical Pendulum on an 
Anisotropic Support, by Maurice F. C. 
Allais, Aero/Space Engng., Rev., vol. 18, 


1959, Pp. 
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The Hund Gravitational Equations and 
Geons, by L1. G. Chambers, Canadian J. 
Phys., vol. 37, Sept. 1959, pp. 1008-1016. 

Sidereal Anisotropy of High Energy Cos- 
mic Rays, by I. Escobar, V: N. Nerurkar 
and R. Weil, Planetary and Space Sci., vol. 
1, no. 3, Aug. 1959, pp. 155-160. 

Nitric Oxide and Molecular Oxygen in 
the Earth’s Upper Atmosphere, by A. 8. 
Jursa, Y. Tanaka and F. LeBlanc, Plane- 
tary and Space Sci., vol. 1, no. 3, Aug. 
1959, pp. 161-172. 

The Relation Between Lightning Dis- 
charges and Whistlers, by Harald Norin- 
der and Edgar Knudsen, Planetary and 
Spuce Sci., vol. 1, no. 3, Aug. 1959, pp. 
173-183. 

Strongly Absorbing Layers Below 50km, 
by Bengt Hultquist and Johannes Ortner, 
Planetary and Space Sci., vol. 1, no. 3, 
Aug. 1959, pp. 193-204. 

On Some Characteristics of the E,- 

layer, by M. Bossolasco and A. Elena, 
Planetary and Space Sc7., vol. 1, no. 3, 
Aug. 1959, pp. 205-212. 
The Thermal Radiation of the Moon at 
1420 Mc/s, by P. G. Mezger and H. 
Strassl, Planetary and Space Sci., vol. 1, 
no. 3, Aug. 1959, pp. 213-226. 

Artificial Electron Clouds, I: Summary 
Report on the Creation of Artificial Elec- 
tron Clouds in the Upper Atmosphere, by 
F. F. Marmo, L. M. Aschenbrand, and J 
Pressman, Planetary and Space Sci., vol. 
1, no. 3, Aug. 1959, pp. 227-237. 

Determination of Upper-atmosphere 
Air Density and Scale Height from Satel- 
lite Observations, by G. V. Groves, Proc. 
Roy. Soc., London, vol. A 252, Aug.. 25, 
1959, pp. 16-27. 

Determination of Upper-atmosphere 
Air-density Profile from Satellite Observa- 
tions, by G. V. Groves, Proc. Roy. Soc., 
London, vol. A 252, Aug. 25, 1959, pp. 
28-36. 

The Earth’s Simple Shadow-Effect on 
Cosmic Radiation, by J. E. Kasper, Nuovo 
Supplemenio, no. 1, 1959, pp. 


‘een and Simple Shadow Cone of 
Cosmic Radiation, by M. Schwartz, 
Nuovo Cimento Supplemento, no. 1, 1959, 
pp. 27-35. 

Propagation of Radio Waves along an 
Inhomogeneous Surface, by E. L. Fein- 
berg, Nuovo Cimento Supplemento, no. 1, 
1959, pp. 60-66. 

Polarization of the Red-oxygen Auroral 
Line, by R. A. Duncan, Planetary and 
Space Sci., vol. 1, no. 2, April 1959, pp. 
112-114. 

Observations of the Russian Satellites 
and the Structure of the Outer Terrestrial 
Atmosphere, by H. K. Paetzold, Planetary 
and Space Sci., vol. 1, no. 2, April 1959, pp. 
115-124. 

Ambipolar Diffusion of a Meteor Train 
and its Relation with Height, by E. L. 
Murray, Planetary and Space Sci., vol. 1, 
no. 2, April 1959, pp. 125-129. 

Absorption and Electron Distribution in 
the F, Layer Determined from Measure- 
ments of Transmitted Radio Signals from 
Earth Satellites, by A. N. Kazantsev, 
Planetary and Space Sci., vol. 1, no. 2, 
April 1959, pp. 130-135. 

Derivation and Analysis of Atmospheric 
Density from Observations of Satellite 
1958 Epsilon, by G. F. Schilling and C. A. 
Whitney, Planetary and Space Sci., vol. 1, 
no. 2, April 1959, pp. 136-145. 

Some Remarks on the Nature and 
Origin of Noctilucent Cloud Particles, by 
Drian Deirmendjian and E. H. Vestine, 
Planetary and Space Sci., vol. 1, no. 2, 
April 1959, pp. 146-153. 


JuNnE 1960 


Trends in Upper Atmosphere Research, 
by R. D. Cadle, SRI J., vol. 3, 2nd Quar- 
ter, 1959, pp. 41-45. 

Rocketsondes for High-altitude Sound- 
ings, by I. G. Poppoff, SRI J., vol. 3, 
2nd Quarter, 1959, pp. 56-61. 

Highly Directional Detector for Cosmic 
Ray Particles, by G. W. Hutchinson, 
Nuovo Cimento, vol. II, no. 3, Feb. 1, 
1959, pp. 377-382. 

An Experimental Analysis of Two Large 
Cosmic Ray Jets, by F. A. Brisbout and 
C. B. A. McCusker, Nuovo Cimento, vol. 
11, no. 4, Feb. 16, 1959, pp. 484—490. 

Gravitational Motion and Radiation I, 
by A. Peres, Nuovo Cimento, vol. IT, no. 5, 
March 1, 1959, pp. 617-627. 

Gravitational Motion and Radiation II, 
by A. Peres, Nuovo Cimento, vol. II, no. 
5, March 1, 1959, pp. 644-655. 


An Astronomical Telescope in Space, by 
R. J. Davis, F. L. Whipple and C. A. 
Whitney, Astron. Sci. Rev., vol. 1, no. 1, 
Jan.—March, 1959, pp. 9-12, 14. 


On the Existence of a Strong Magneto- 
Ionic Effect Topside of the F Maximum of 
the Kenelly-heaviside Layer, by P. R. 
Arendt, J. Appl. Phys. vol. 30, no. 5, 
May 1959, pp. 793-794. 


Near Infrared Solar Radiation Measure- 
ments by Balloon to an Altitude of 100,000 
Feet, by D. M. Gates, D. G. Murcray, 
C. C. Shaw and R. Herbold, J. Optical 
Soc. Amer., vol. 48, no. 12, Dec. 1958, 
pp. 1010-1017. 


Scientific Problems in Cislunar Space 


and Their Exploration with Rocket Ve- 
hicles, by S. F. Singer, Astronautica Acta, 
vol. 5, no. 2, 1959, pp. 116-125. 


Techniques of Analyzing Terrestrial 
Radio and Optical Observations of Earth 
Satellites, by R. H. Merson, Astronautica 
Acta, vol. 5, no. 1, 1959, pp. 26-39. 

The Physical Nature of the Surface of 
the Moon, by Gilbert Fielder, Brit. Inter- 
planet. Soc., J., vol. 17, March-April 1959, 
pp. 57-58. 

The Detection of Planets at Interstellar 
Distances, by M. H. Briggs, Brit. Inter- 
planet. Soc., J., vol. 17, March-April 
1959, pp. 59-60. 

Size of Irregularities in the E region of 
the Ionosphere, by M. Srirama Rao, 
Canadian J. Phys., vol. 37, May 1959, pp. 
557-568. 

Transient Decreases in Cosmic Ray In- 
tensity During the Period October 1956 to 
January 1958, by A. G. Fenton, K. G. 
McCracken, D. C. Rose and B. G. Wilson, 
Canadian J. Phys., vol. 27, May 1959, pp. 
569-578. 

Further Observations from the Russian 
Satellites, by H. K. Paetzold, Raketen- 
technik und Raumfahrtforschung, vol. 3, 
no. 2, April 1959, pp. 45-49. (In Ger- 
man.) 

Balloon Astronomy, by Martin and 
Barbara Schwarzchild, Scientific American, 
vol. 200, May 1959, pp. 52-59. 

Artificial Satellites and Relativity, by 
V. L. Ginzburg, Scientific American, vol. 
200, May 1959, pp. 149-160. 

Rocket Astronomy, by Herbert Fried- 
man, Scientific American, vol. 200, June 
1959, pp. 52-59. 

Artificial Modification of the Earth’s 
Radiation Belt, by S. F. Singer, J. Astron. 
Sci., vol. 6, no. 1, Spring 1959, pp. 1-10. 

The Sun and the Upper Atmosphere, by 
Herbert Friedman, ASTRONAUTICS, vol. 4, 
July 1959, pp. 20-22, 55-56. 

Density and rs a erature of the Upper 
Atmosphere, by Robert Jastrow, AsTRo- 
NAuTIcs, vol. 4, July 1959, pp. 24-25. 


13 Noted Scientists Present 


A SPACE AGE SYMPOSIUM... 


Robert Jastrow 
Fred L. Whipple 
Thomas Gold 
Eugene Parker 
James A. Van Allen 
N. C. Christofilos 
Homer E. Newell 
Gerard P. Kuiper 
Harold C. Urey 

G. de Vaucouleurs 
Herbert Friedman 
Leo Goldberg 

J. W. Townsend, Jr. 


THE EXPLORATION OF SPACE 


edited by Robert Jastrow 
illustrated, $5.50 

Van Allen discusses the geomagneti- 
cally trapped corpuscular radiation; 
Christofilos reports on “The Argus 
Experiment;” de Vaucouleurs, on Mars 
and Venus; Goldberg, on astronomy 
from satellites and space vehicles; 
Whipple, on solid particles in the 
solar system. Other vital topics include 
rocket astronomy, plasma and magnet- 
ic fields in the solar system, the U.S. 
space-exploration program, and the 
U.S.S.R. report of an active volcano 
on the moon. 
Order Today From Your Nearest 

Bookstore Or Write Dept. ARS-1. 


The Macmillan Company 


60 FIFTH AVENUE, NEW YORK 11, N.Y. 


Aeronautical 
Paperbacks 


New paperbacks—published in 
response to a growing demand 
for the most essential sections of 
the HIGH SPEED AERODYNAMIC 
AND JET PROPULSION series— 
bring to the student and re- 
search engineer directly useful 
information, at prices they can 
easily afford. 


1. LIQUID PROPELLANT ROCKETS 
David Altman, James M. Carter, 
S. S. Penner, Martin Summer- 


field. 196 pp. $2.95 


2. SOLID PROPELLANT ROCKETS 
Clayton Huggett, C. E. Bartley 
and Mark M. Mills, 176 pp. $2.45 


3. GASDYNAMIC DISCONTINUITIES 
Wallace D. Hayes. 76 pp. $1.45 


4. SMALL PERTURBATION THEORY 
W. R. Sears. 72 pp. $1.45 


5. HIGHER APPROXIMATIONS IN 
AERODYNAMIC THEORY 
M. J. Lighthill. 156 pp. $1.95 


6. HIGH SPEED WING THEORY 
Robert T. Jones and Doris 
Cohen. 248 pp. $2.95 
At your bookstore 
Princeton University Press 
New Jersey 
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